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HAMILTON HARBOUR REPORT 


SUMMARY 


The studies in Hamilton Harbour in 1977 were a continuation of 
studies from earlier years. (MOE 1975, 1977, 1978.) In addition to 
the work carried out by the Water Resources Branch, investigations 
of plankton densities and speciation and their correlation with 
nutrient concentration and light regimes were done under a contract 
by Dr. G. Harris of McMaster University. An extensive study of 
loadings to the Harbour from industry, municipal sewage treatment 
plants, tributaries and storm sewers was compiled under a contract 
by Dr. W. Snodgrass of McMaster University. This is included as 
Volume II of this report. The main findings are summarized below. 


VOLUME I 
A. ARTIFICIAL MIXING 


Operationally, this was the best year to date, although the airflow 
was below the estimated critical level required for successful 
turnover of the water column and reaeration. The results expected 
were: an increase in dissolved oxygen levels; reduction in the 
levels of dissolved heavy metals, nutrients and dissolved organics; 


and improved light penetration. 


The results of the bi-weekly sampling for water quality in the 

Harbour did not show any evidence that the temperature and oxygen 
profiles had been affected by the mixing. The hypolimnion oxygen 
concentration increased at the end of July as it had in previous 


years. 


No relationships between mixing and the concentrations of other 
water quality parameters such as nutrients, suspended solids, 
phytoplankton and bacteria could be observed. The one possible 
exception was the increase in the number of nitrogen cycle bacteria. 


B. OXYGEN BUDGET 


The average oxygen stocks in the Harbour were the same as in 
previous years, reaching a minimum in early July when the 
epiliminion was also undersaturated. The lenses of oxygen-rich 
water in the hypolimnion during June and July are, in all 
likelihood, lake water entering via the Burlington Canal. 


The vertical turbulent fluxes of oxygen were calculated for 
different seasons using a model which divides the Harbour vertically 
into layers of equal thickness. The layers are considered as 
horizontally homogeneous. The vertical diffusion coefficient, D, 
was calculated from the fluxes and the known oxygen gradients. A 
very pronounced minimum in D was found at the thermocline, but, 
because of the high oxygen gradients, the oxygen flux was still 
considerable and at no time during the sunmer was the hypolimnion 
effectively isolated from the epilimnion. 


A simple single box predictive model for oxygen concentrations was 
developed, using sediment oxygen demand and water column oxygen 
demand as sinks and reaeration, photosynthesis and lake-harbour 
exchange as sources. Quite good agreement between predicted and 
observed mean oxygen concentrations was achieved. The model 
Suggests that to reach an oxygen concentration of 5 mg/L, which is 
the Ontario objective for the protection of warm-water biota at 
temperatures between 10°C and 20°C the water column demand 

Should be reduced by 20%. 


C. WATER CHEMISTRY 


In 1976 and 1977, four stations were sampled approximately bi-weekly 
to determine trends in water quality, the oxygen budget and the 
effects of abatement procedures. In addition, depth profiles of 
temperature, oxygen, Eh, pH and conductivity were measured weekly. 
In 1977 a survey of stormwater runoff at 12 locations in the harbour 
was done. 


Oxygen has already been dealt with in Section B. Of the nutrients, 
ammonia was much higher in the spring of 1977 than in the previous 
two years. The source is apparently the Hamilton Sewage Treatment 
Plant. The year-to-year or seasonal differences in total phosphorus 
concentrations were not statistically significant. Silica 
concentrations were again higher in the hypolimnion than in the 
epilimnion. 


The heavy metals iron and zinc show no consistent spatial and 
temporal variations. Zinc is present at concentrations high enough 
to be potentially toxic to algae. 


Organic matter in the harbour was determined as chemical and 
biochemical oxygen demand and total organic carbon. These 
parameters are not correlated to one another since different tests 
measure significantly different substances. Each parameter varies 
Significantly with sampling date. The maximum BOD values were found 
in June at the time corresponding to the maximum nitrate 
concentrations. Dissolved organic carbon , however, correlated with 
BOD, suggesting that this form of organic carbon plays an important 
role in the oxygen demand of the water column. Most of the 
industrial, municipal and stormwater outfalls are sources of 
dissolved organic carbon. 


D. SEDIMENT CHEMISTRY 


Sediment samples were collected at 25 stations in May and October of 
1977 and at 4 stations monthly. There were no significant 
differences, except for ammonia, seasonally or from 1976 to 1977. 
The sediments are highly polluted with zinc, lead and chromium and 
to a lesser extent cadmium and mercury. The highest concentrations 
are found in deep waters and adjacent to the outfalls. Factor 
analysis classified the sediments into four main groups: the deep 
water group, those close to the north shore, those near outfalls and 
those in the canal region which is influenced by lake-harbour 


exchange. 


E. PHYTOPLANKTON AND BACTERIA 


The highest phytoplankton cell densities were observed in the spring 
with a secondary maximum during the summer. The maximum cell 
concentrations at any period were found in the upper three to four 
metres of the water column, which was also the euphotic depth. The 
community structure was highly dynamic, influenced as it is by the 
vertically and horizontally unstable physical structure of the water 


column. 


The dominant species throughout most of the year were the 
Chlorophyta (green algae), the Cryptophyta were the next important 
group, with Chrysophyta (diatoms) present in small numbers. The 
bacterial populations were dominated numerically by heterotrophs. 
These were followed in descending numbers by sulphur oxidizers, 
nitrifiers, and sulphur reducers. Although processes and reaction 
rates cannot be accurately assessed from bacterial numbers, the 
results suggest that organic carbon breakdown and sulphur oxidation 
are the major sources of the large and sustained water column oxygen 
demand in the Harbour. 


F. 2-DIMENSIONAL NUMERICAL MODEL 


The two-dimensional numerical model described in earlier Ministry of 
Environment publications has been modified to include real-time data 
for industrial inputs and lake-harbour exchange. Nevertheless, the 
model simulation is as yet not satisfactory as it does not describe 
short-term variations at all, and the agreement between observed and 
models results is better at grid squares removed from meter 
locations. The major deficiency in this modelling is the use of a 
two-dimensional model to describe a three-dimensional stratified 
Situation. 


The preliminary simulation of total dissolved solids distribution in 
the harbour was done, and qualitative agreement was obtained, 
showing high concentrations in the southeastern region and lower 
concentrations due to lake water intrusions at the canal. 


G. 3-DIMENSIONAL NUMERICAL MODEL 


The three-dimensional numerical model of hydrodynamics, temperature 
and water quality as developed by Leendertse et al (1973) has been 
modified for use in lake systems. The theory of the model is 
described in detail, with particular reference to the physical 
parameters for which calibration values must be specified. 


Initial calibration runs of the hydrodynamic portion of the model 
have been performed for Hamilton Harbour using real-time input 

data. Model performance was assessed by comparing modelled currents 
to those measured at a harbour location not used as model input. 
Attempts were made to improve the fit to observed currents by 
varying the values of the input physical parameters. 


Poor model performance was likely the result of vertical instability 
in the harbour. Water flowing past a current meter which has a 
fixed vertical position may be representative of different model 
layers at different times. 


The results suggest that any future modelling should involve a 
larger number of model horizontal layers, as well as vertical 
strings of current meters for the calibration process. 


H. WATER MOVEMENTS 


Currents were measured in both 1976 and 1977 in mid-harbour and in 


the Burlington and Desjardins Canals. 


In summer, .the mid-harbour currents were to the southerly quadrant 
for both years. Spectral analyses show diurnal and semi-diurnal 
periods and the effects of Lake Ontario. 


The Burlington Canal flows were uni-directional and towards the Lake 
in the sunmer of 1976. During the same period in 1977, the flow was 
outward in the upper layers but inward below 3.5 metres. 


The net exchange for the period from June to November 1976, 
calculated using the excursion-episode method, varied from 6 to 45 
m?/s representing 0.2% to 1.4% of the harbour volume per day. Net 
flows calculated from mass balances of conservative substances, 
assuming a steady-state for a year, are about 35% to 60% of the 
flows in this report. This discrepancy between the two methods 
still has to be resolved. 


The net flow through the Desjardins Canal varied from 0.6 to 1.5 
m?/s representing 0.02% to 0.05% of the harbour volume per day. 


VOLUME II 
INPUTS OF MATERIALS TO THE HARBOUR 


The inputs of pollutants from tributaries, storm sewers, industrial 
outfalls and municipal sewage treatment plants to the harbour were 
obtained under contract by W. J. Snodgrass, McMaster University. 

The results are given in Volume II of this report. These include 
estimates of water input from all major watersheds, sewage treatment 
plants, lake-harbour exchange, rainfall and storm sewers from below 
the escarpment. The author also points out that exchange flows 
calculated on the excursion-episode method can lead to unacceptably 
large theoretical lowerings of the harbour level, as the inflows are 
not sufficient to balance the outflows. 


The estimates of industrial loadings from the two integrated steel 
mills - Stelco and Dofasco suffer from scarcity of flow and 
concentration measurements. The major net loadings from Stelco are 
suspended solids, BOD and iron; from Dofasco, suspended solids, COD, 
and sulphates. In addition, ammonia, phenols and cyanides are added 
from both. 


The Hamilton sewage treatment plant removes about 80% to 90% of 
incoming BOD, COD, suspended solids, total phosphorus and inon,+ but 
only about 5% to 20% of the nitrogen. The Burlington sewage 
treatment plant removes 90% of suspended solids, BOD and total 
phosphorus through chemical precipitation and 70% of total nitrogen 
(biological removal). The Dundas sewage treatment plant is a large 
source of phosphorus and ammonia, much of which is removed in Cootes 


Paradise. 


Hamilton has a combined sewer system. The overflows from above the 
escarpment drain into Redhill and Chedoke Creeks; from below, 
directly into the harbour. The studies of stormwater inputs are 
insufficient to give a good estimate of loadings, but those draining 
directly into the harbour contribute, at the most, 5% of total 
pollutant load. The author estimates that the municipal sewage 
treatment plants are the major contributors (greater than 50% of 
total load) of total phosphorus, soluble phosphorus, total nitrogen, 
nitrate, total dissolved solids and chloride. Industrial and 
municipal sources each contribute 30% to 40% of total organic carbon 
and BOD, whereas industry is a major source of COD, suspended solids 
and phenols. Streams are significant sources (20% to 25%) only for 
nitrates, total dissolved solids and chlorides. 


The author concludes that improved control of municipal and 
industrial point sources would produce the most significant 
improvements in water quality, whereas control on combined sewers 
would effect only minor improvements. 
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HAMILTON HARBOUR STUDY 
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VOLUME I 


SECTION A 


ARTIFICIAL MIXING 


dite. MUOSRAH WOTSIIMAH 


SUMMAR Y 


The primary artificial mixing system employed was identical to that 
of 1976. The system was operational at the beginning of June. At 
this time, water column dissolved oxygen concentrations were already 
severely reduced. The theoretically required air flow for maximum 
column destabilization was not obtained until late July. This flow 
was sustained for only a few weeks due to subsequent system 
failures. This year's activity represents the best operational 
record to date, but it was inadequate from both a theoretical and an 
observational viewpoint. 


Neither vertical thermal structures nor dissolved oxygen 
concentration could be directly related to aerator activity due to 
the harbour's naturally vertically diffuse characteristics. The 
customary recovery in the severely depleted bottom oxygen stock was 
again seen at the end of July. This recovery may have been due to 
aerator activity, however a clear cause and effect relationship 
cannot be established at this time. 


Bi-weekly samples of some water column chemical parameters, 
phytoplankton and bacteria were taken. Despite significant changes 
in some chemical parameters, such as the elevation of ammonia 
levels, no direct relationships between water column chemistry and 
artificial mixing could be established. The same was true for both 
the phytoplankton and the bacteria measurements. There did, 
however, appear to be more subtle or indirect effects of artificial 
mixing on the harbour. Despite higher in situ concentrations of 
chemical parameters such as ammonia, the assimilative capacity of 
the harbour did not diminish appreciably. Also, nutrient cycle 
bacteria appeared to respond to the aerator operation by maintaining 
and even augmenting their standing stocks. 
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INTRODUCTION 


Water quality studies have been conducted in Hamilton Harbour since 
1972. On the basis of the 1972 to 1974 information it was 
recommended that artificial mixing be instituted as a means for 
improving the poor water quality which existed. It was anticipated 
that the aerator system would increase dissolved oxygen 
concentrations, reduce the levels of heavy metals, nutrients, and 
dissolved organics, and improve the light penetration. Extensive 
monitoring and process oriented studies were conducted to examine 
the impact of aeration on the harbour. These programs have also 
disclosed in greater detail the physical, chemical and biological 
processes of the harbour. The monitoring program consisted of a 
bi-weekly vertical sampling at the four stations (Figure 1) for 
nutrients, metals, bacteria, phytoplankton, alkalinity, pH, oxygen, 
chlorophyll and temperature. 


In addition to the artificial mixing device operated by the 
Ministry, a secondary mixing technique developed by Professor M. 
Baird of McMaster University was examined. 


THE MIXING SYTSTEMS 


The primary mixing system employed was identical in design to that 
used in 1975 and 1976, and has been described in detail in the 
previous reports (MOE, 1977, 1978). The compressor and lines were 
located off the north face of the Stelco ore dock as in 1976 
(Figure 1) however some of the existing lines were replaced. 


The operational record is displayed in Figure 2. The 1976 program 
represents the best performance record to date, with the system 
functioning from the beginning of June until early September when 
operations were terminated. The system was completely 
non-functional for brief periods at the beginning of July and the 
beginning and end of August. Compressed air flow through the 
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various lines was briefly interrupted over the sunmer and line #1 
was inoperative in the middle of August. As a result of these 
problems, the average daily flow from the system (Figure 2) was 
below the critical, theoretical requirement (300 cfm) throughout 
most of the summer (MOE, 1978). Therefore, it is not surprising 
that the system could neither promote an oxic environment in the 
hypolimnion nor a large increase in the dissolved oxygen stocks of 


the harbour. 


The secondary mixing technique employed was that of propelling a 
series of tori (ring-shaped water masses) downward from a special 
generator mounted on a raft. The vortex generator consisted of a 
compressor, a timer, a surge tank and the vortex tank itself which 
was an open ended cylinder inverted into the water. Pulses of air 
were released from the surge tank at timed intervals and fed into 
the vortex chamber which generated ring vortices of water. The 
vortex of water so generated was 25.4 cm in diameter and displaced 
approximately 9 m?/hr of water downwards. Entrainment effects 

must be added to this estimate and an order of magnitude estimate of 
transport of aerated water downwards would be 20 m>/hr. The 

device was operated for a few hours on several days during July and 
August in a water column depth of 18 m. Hourly recordings of 
dissolved oxygen concentration and temperature were taken at depths 
of 0, 5, 10 and 15 metres. 


TEMPERATURE AND DISSOLVED OXYGEN 


A definite, harbour-wide surface cooling was noted at the end of 
July and the beginning of August (Figure 3). This surface cooling 
was accompanied by a deepening of the thermocline and a slight 
warming of bottom waters. The phenomena corresponded to that period 
of peak aerator output (Figure 2) where delivery rates exceeded the 
theoretical critical flow requirements. 


Local surface coolings were observed at all stations periodically 
throughout the summer and are unlikely to be attributable to aerator 
Operations. Since the harbour exhibits naturally diffuse vertical 
thermal gradients, periodic changes in the temperature profiles 


cannot be conclusively related to aerator activity. Large 
differences in both the absolute temperatures and vertical 
temperature gradients often existed between stations which further 
prevented a clear interpretation of the aerator effects. Mixing did 
not stimulate vertical exchange sufficiently to achieve thermal 
homogeneity. 


The oxygen profiles were likewise unrelated to artificial mixing. 
The oxygen demand rates were large and exceeded the natural and 
artifically induced supply rates (refer to Section 'B'). The 
dissolved oxygen concentrations of the harbour declined rapidly from 
Saturation levels during the spring to bottom anoxia and surface 
sub-saturation levels during the summer. A slight recovery in both 
the bottom and surface concentrations was evident near the end of 
July and the beginning of August. After autumn turnover, which 
occurred in late September, the vertical oxygen gradient was 
eliminated. The oxygen concentrations did, however, remain 
sub-saturated until sometime after survey activity ceased at the end 
of October (refer to Figures 4 and 5). 


The total oxygen stock of the harbour declined rapidly from 3.1 x 
10° kg of oxygen in April to 0.9 x 10° kg of oxygen in early 

July (Figure 6a, Table 1). The original saturated stock levels were 
not again attained during the survey period (Figure 6b). Total 
oxygen stock did, however, show an increase in late summer which 
corresponded to that period of maximum aerator output. The overall 
effect of the aerator system on the oxygen stocks was far less than 
was originally anticipated. This was partly due to late initiation 
of operation and partly due to the systems inability to deliver 
compressed air at the required flow rates. 


The results of the vortex mixing experiment showed that there was a 
Slow but definite increase in the dissolved oxygen concentrations 
directly beneath the raft. Vertical mixing was also demonstrated by 
the temperature profile data.. A complete analysis of this 
experiment, along with a discussion of laboratory experiments and 
vortex mixing theory can be found in Rohatgi, Baird, and Wairegi 
(1979). 


NUTRIENTS 


Filtered reactive phosphorus (FRP) concentrations (as P) ranged from 
less than 1 ug/L to 96 ug/L. No definite seasonal trend was evident 
from the survey data. Also, no spatial trends could be found in 
either the horizontal or vertical dimensions. In general, the FRP 
concentrations were greater than those of 1976, particularly at 
stations 4 and 252. 


The total phosphorus (TP) concentrations (as P) ranged from 20 to 
128 ug/L and, like the FRP, showed no definite trends in space or 
time. Unlike the 1976 results, no pattern of higher concentrations 
associated with sediment disturbance was seen. Surface mean 
concentrations were in general as high or higher than those of the 


whole water column. 


Neither total nor dissolved phosphorus concentrations appeared to be 
related to aerator activity. Concentrations were not reduced 
Significantly when compared to previous years, despite the 
significant improvement in aerator operations. This may have been 
due to the large and continuous loadings of phosphorus to the 
harbour (refer to Section 'H'). 


Anmonia concentrations which were approximately uniform with depth, 
ranged from 0.012 to 6.30 mg/L throughout the year. Concentrations 
were highest in the spring and rapidly declined until the end of 
June. Thereafter, concentrations decreased much more slowly. While 
the seasonal trend resembled those of previous years, it can be seen 
that the spring concentrations were much higher in 1977 (refer to 
Section 'C'). The greatest concentrations were found in the eastern 
end of the harbour, particularly at station oe 


The nitrite concentrations in the harbour ranged from 0.011 - 0.370 
mg/L during the sampling period. The lowest concentrations were 
found in the hypolimnion during the summer. Concentrations changed 
inversely as the ammonia concentrations during the first part of the 
year. After the ammonia levels became greatly reduced, no direct 
relationship could be ascertained between the two parameters. 
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Nitrate concentrations displayed a great deal more vertical tha.. 
horizontal variation, in contrast to other nitrogen parameters. The 
concentrations ranged from 0.20 - 3.37 mg/L over the year with the 
lowest concentrations being found near the bottom of the harbour. 
The mean harbour concentrations dropped slightly from early to late 
spring and then increased to a mid-summer maximum. Although the 
values remained high throughout the late summer and autumn, they 
oscillated from one sampling date to the next. 


A total Kjeldahl nitrogen, whose concentrations ranged from 0.21 - 
7.00 mg/L, showed a seasonal trend similar to ammonia. The 
concentrations did, however, show a great deal of variation in both 
the vertical and horizontal dimensions. 


While the total nitrogen concentrations were in general greater in 
1977 than in 1976, no relationship between either the concentration 
changes or speciation and artificial mixing could readily be seen. 


Silica concentrations, which ranged from 1.10 to less than the 
detection limit of 0.05 mg/L, were highest in the spring and lowest 
during the mid-summer period. The dramatic decline in silica, which 
was observed in late June to early July, coincided with the diatom 
bloom (refer to Section 'E'). During this period the highest 
concentrations were observed near the bottom of the harbour. These 
conditions persisted until Autumn overturn. Therefore, no 
artificial mixing effects on the concentration or distribution of 


silica were evident. 
ORGANIC CARBON, BOD AND COD 


Total organic carbon (TOC) concentrations ranged from 1-16 mg/L. 
There were no consistent trends in the concentrations with time, 
station or depth. The lowest values observed occurred on August 10 
and these were similar throughout the whole of the harbour. The 
concentrations, however, had reattained previous levels within two 


weeks. 
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Analysis of the total filtered organic carbon (i.e. organic carbon 
compounds below the 1-2 micron size) showed that it comprised 76-96% 
on average, of the total organic carbon. Therefore, the changes 
observed in the total organic carbon concentrations were also seen 
in the filtered organic carbon concentrations. The mean harbour 
concentrations observed were similar to those in 1976. None of the 
changes seen in the organic carbon parameters could be related to 


aerator activity. 


Total inorganic carbon concentrations ranged from 21-29 mg/L, being 
slightly lower through July and August than in the spring and 
autumn. The concentrations showed some depth variation throughout 


the year. 


The biochemical oxygen demand (BOD, ) showed no consistent seasonal 
trend. The BOD's were both horizontally and vertically variable. 
The values ranged from 1.0 to greater than 8.0 mg/L and in general, 
the highest values appeared in the epilimnion. The mean values were 
somewhat higher than those of 1976. BOD values were unrelated to 
the total organic carbon concentrations. 


The chemical oxygen demand (COD) also showed no definite seasonal 
trend. Its values ranged from 10-150 mg/L and were equally variable 
horizontally and vertically. In general, the highest values were 
observed in the epilimnion. The mean harbour values were similar to 
those observed in 1976 during the spring and early part of the 
summer. The late summer and autumn values, however, were 
consistently higher in 1977. The COD showed no direct relationship 
with either the TOC or the BOD... 
No direct relationships between artificial mixing and TOC, BOD, or 
COD could be found. ) 
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ALKALINITY, PH AND QONDUCTIVITY 


The carbonate alkalinity of the harbour decreased gradually from a 
high of 123 mg/L in the spring to 87 mg/L in the late summer. The 
harbour waters are therefore well buffered and the pH changes due to 
primary production or artificial mixing would be minimized. The pH 
ranged fran 7.15-8.60 and was highest in the epilimnetic waters. 
Values in general, declined fran spring to autumn, however, the 
highest levels were observed in the surface waters in late July and 
early August which coincided with large increases in the 

phytop! ankton commmity. 


Conductivity ranged from 410-650 u Siemens cm. The mean column 
values declined throughout the year and the highest values were 
found in the surface waters. 


The values and trends observed for all of the above parameters were 
consistent with observations made in previous years. No aerator 
effects were apparent despite improved operations. 


TURBIDITY AND SUSPENDED SOLIDS 


The turbidity levels in the harbour ranged fron 1.5 to 7.5 Formazin 
units which were sanewhat greater than the 1.5-2.3 (Formazin units) 
range of 1976. Turbidity levels were in general, greatest in the 
epilimnion. There was a good deal of variation in turbidity between 
stations and with depth. The stations at the extreme end of the 
harbour, 270 and 4 (refer to Figure 1), had the highest turbidity 
levels. Values were greatest at all stations in the late spring and 
early summer, however no definite relationship between artificial 
mixing and turbidity was found. 


The suSpended solids ranged fran 2-14 mg/L over the year. The 
highest values occurred in the epilimnion but in general, showed 
both spatial and temporal variability. The average harbour values 
were bimodally distributed, with peaks in early May and late July. 


Turbidity and suspended solids were linearly correlated, 


Turbidity = 0.34 (suspended solids) - 1.10; 


At the 0.001 level with r = 0.80. This relationship shows that when 
turbidity is 0, the suspended solids content of the harbour is 1.10 
mg/L. This value was similar to that found in 1976 (1.78 mg/L) and 
suggests that particles less than 2 microns were abundant. These 
particles have a forward scattering effect on penetrating light and 


are therefore not measured as turbidity. 


The slope of the regression line (0.34) was significantly different 
from that obtained in 1976 (0.97). The reason for this change in 
slope, which implies that turbidity values increased much more 
slowly with increased suspended solids concentrations in 1977 than 
in 1976, is not known. 


HEAVY METALS 


Iron and zinc, which had been identified as major environmental 
concerns in previous years, were monitored bi-weekly at the four 
Stations. 


Iron concentrations ranged from 0.04-1.3 mg/L. The values were in 
general greater at station 4 than other stations and were in excess 
of the IJC guidelines of 0.30 mg/L throughout the harbour. No 
definite seasonal or vertical trend was apparent in the data. 


Zinc was again present in high concentrations with a range of 
0.02-0.46 mg/L. Although the values were spatially and temporally 
variable, the average column concentrations normally exceeded the 
IJC objective of 0.03 mg/L. No effects of artificial mixing on 
either the zinc or iron concentrations were apparent. 


Other metals such as manganese, copper and mercury were sampled too 
infrequently to relate them to aerator activity. 


A more detailed analysis and discussion of the chemical water 
quality can be found in Section C of this report. 
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BIOLOGICAL RESPONSES 
BACTERIAL POPULATIONS 


Bi-weekly samples were taken for fecal coliforms, heterotrophs, 
Sulphur cycle and nitrogen cycle bacteria from stations 4, 258 and 
270 at the surface, mid-depth and bottom. Similar samples were also 
taken from the surface and bottom at station 252. 


The fecal coliform population densities, which indicate domestic 
Sewage contamination, showed a slight increase over the 1976 

levels. The highest concentrations were observed at station 4. The 
population densities were highly variable in space and time which 
precluded a clear linkage of population size and distribution with 
aerator activity. Fecal coliform levels in the surface water 
frequently exceeded the 100 counts per 100 ml level which is the 
upper limit permissible for body contact recreation. 


Heterotroph population densities showed a significant increase over 
the previous year at all stations. This trend of increase in 
population densities, which has been observed since 1975, may be due 
to aerator activity, however, changes in the population densities 
observed within a single year cannot be directly linked to aerator 
activity. The highest levels of heterotrophs were generally found 
at station 4 (refer to Table 1). 


Both the sulphur oxidizers and sulphur reducers showed notable 
declines in population size. This may have been due to the decline 
in sulphur, particularly reduced sulphur, concentrations observed in 
the sediments in 1976. Sulphur oxidizer levels were greatest at 
station 4, while the sulphur reducer levels were greatest at station 
252 (refer to Table 1). A good deal of spatial and temporal 
variability was seen in both populations and no direct effect of 


aerator activity could be shown. 
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The nitrogen cycle bacteria had shown a dramatic increase in 
population densities since 1976. This large increase was probably 
due to aerator activity. The greatest increase during the year 
centred around peak aerator operation. This increase was most 
likely a result of artificial mixing since the ammonia levels were 
coincidentally at their lowest levels for the year. 


PHYTOPLANKTON 


Phytoplankton samples were collected weekly at station 258 from the 
beginning of May until the middle of October. Samples were also 
collected for chlorophyll a analysis, which is used as an indicator 
of algal abundance, however, due to analytical difficulties the 
results cannot be used. The phytoplankton numbers demonstrated a 
large amount of vertical variability from week to week, but neither 
total algal abundance nor vertical arrangement, nor species 
composition could be directly related to aerator activity. 


ZOOPLANKTON 


Zooplankton were monitored, as well as surface phytoplankton 
throughout the harbour and primary productivity was measured by 
McMaster University. The results of this program are discussed in 
Harris, Piccinin, Haffner, Snodgrass and Polak (1979). 


OTHERS 


Bottom fauna were collected at various times during the year by the 
West Central Region (MOE), however, analysis is aS yet incomplete. 


Various fish species were examined for organochlorine contaminants 
during the year and the results of these findings are available from 
the Ontario Ministry of the Environment, Water Resources Branch. 
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Table 1: 


Bacterial Populations of Hamilton Harbour - 1977 


(Annual Geometric Means), Counts/100 mL 


Fecal Coliforms 
Heterotrophs 
Sulphur Oxidizers 
Sulphur Reducers 
Nitrosomonas 
Nitrobacter 
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Sulphur Reducers 
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Nitrobacter 
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Nitrobacter 
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HAMILTON HARBOUR STUDY 


1977 


VOLUME I 
SECTION B 


OXYGEN BUDGET 


SUMMARY 


The average oxygen stocks in the Harbour were the same as in 
previous years, reaching a minimum in early July when the epilimnion 
Was also undersaturated. The lenses of oxygen-rich water in the 
hypolimnion during June and July are, in all likelihood, lake water 
entering via the Burlington Canal. 


The vertical turbulent fluxes of oxygen were calculated for 
different seasons using a model which divides the Harbour vertically 
into layers of equal thickness. The layers are considered as 
horizontally homogeneous. The vertical diffusion coefficient, D, 
was calculated from the fluxes and the known oxygen gradients. A 
very pronounced minimum in D was found at the thermocline, but, 
because of the high oxygen gradients, the oxygen flux was still 
considerable and at no time during the summer was the hypolimnion 
effectively isolated from the epilimnion. 


A simple single box predictive model for oxygen concentrations was 
developed, using sediment oxygen demand and water column oxygen 
demand as sinks and reaeration, photosynthesis and lake-harbour 
exchange as sources. Quite good agreement between predicted and 
observed mean oxygen concentrations was achieved. The model 
suggests that to reach an oxygen concentration of 5 mg/L, which is 
the Ontario objective for the protection of warm-water biota, the 
water column demand should be reduced by 20%. 
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INTRODUCTION 


A rather extensive study of oxygen conditions in Hamilton Harbour 
was conducted in 1976 and was published in the 1976 Hamilton Harbour 
Report(MOE, 1978). In 1977, weekly oxygen vertical profiles were 
collected at the four main sampling stations (Figure 1). This 
section summarizes some of the findings of the oxygen budget of the 
harbour in 1977 and shows some of the new techniques applied to the 
data. 


HAMILTON HARBOUR DISSOLVED OXYGEN STOCKS 


The oxygen profiles measured weekly in 1977 were used to calculate 
oxygen contents of the harbour by the same method as described in 
the 1976 report. The total amount of oxygen dissolved in the 
harbour waters is given in Figure 2. The thick line in this figure 
indicates the 1977 dissolved oxygen stocks. As can been seen, the 
seasonal changes of the 1977 stocks are similar to those of the 
previous years. The minimum in 1977 was reached in early July, 
which is later than in 1976 but earlier than in 1975. The minimum 
oxygen content of the harbour is somewhat higher than in the 
previous years but the difference is not substantial. Following the 
July minimum, the DO stock increased slowly until the end of August, 
and again decreased to a secondary minimum in mid-September. This 
secondary minimum could be related to oxygen depletion occurring 
after the cessation of artificial mixing, but before the 
commencement of natural overturn. The late summer decrease seems to 
be a regular pattern as it was also observed in 1975 and 1976. 


Figure 3 shows the average 1977 oxygen concentrations in the 
epilimnion (0 to 3 m from surface) and hypolimnion layers (7 m - 
bottom). The average hypolimnion concentration was below 5 mg L7 
between the middle of June and the end of September. The Ontario 
objective for the protection of warm-water biota is 5 mg L71 when 
the temperature is between 10 and 15°C. The average epilimnion 
concentration was always greater than 5 mg L7l, reaching a minimum 
slightly below 6 mg L~+ near the end of September. The epilimnion 
is, however, below saturation from the end of May to the middle of 
October. 
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VERTICAL EXCHANGE 


A simple model of vertical oxygen exchange can be built using 
previously obtained information on sediment and water column oxygen 
demands. The vertical turbulent exchanges can then be calculated 
from the model results. 


The water reservoir is divided into a finite number of layers of 
equal thickness (Figure 4). Any particular layer i has a volume 
Vj and upper and lower boundaries of area Aj4y] and Aj, 
respectively, because of the slope of the basin. The area of 
sediment exposed in the layer is Sj and Ajy] = Aj+Sj, except 

for the bottom layer where A]=Sj=A2. The oxygen concentration 
is uniform in each layer and equal to (Cj. 


In this model all oxygen in the layer i is assumed to enter by 
vertical exchange from layer i+] (flux rj4,). Oxygen is being 
lost by a) vertical transfer to the lower layer i-1l CELO er ern 

b) by water column oxygen demand (WCOD) and c) by sediment oxygen 
demand (SOD). 


Thus, the change in the dissolved oxygen stock in layer 1 (AM; ) between 
two dates t; and tg is 


AMs= Ri41] - Ry - SOD; - WCOD; (1) 
where Rj, ) = rate of oxygen input from layer itl 
Ri = rate of oxygen loss to layer i-l 


For input from layer itl 


where r;,] 1s the flux 
Assuming that rj,, is constant between the two dates. 
Rit = Aisi Tita (t2 - t)) (2) 
Loss of oxygen to the lower layer can be expressed similarly as 
Rao Ap i (ty 57) (3) 


Amount of oxygen used by sediments 


t 
SOD; = ssf dt (4) 


oI 
Sediment oxygen demand k, in Hamilton Harbour can be approximated 
(Polak and Haffner, 1978) by equation 
k, = 0.72 +0.26c (g 09 m-2a7!) ©) 
where c is the oxygen concentration at time t. 


Assuming linear change of oxygen concentration between the dates 
t,; and t9 


c = COL = Ci7 Cty = tips (co = Gu) aH Cy (6) 


Substituting (5) and (6) into (4) 


c 
SOD = ssf fo.n > Ondo £7 (co - c,)/(t - t)) x 
1 


0126) cytbOs 26site lence eq diiGtoee 1) dt 
or after integration and rearrangement 


SODy =9S, Cty o= ty) C0728 0013. (ch 548, 7) 


(7) 
Similarly, using the equation for water column oxygen demand in 
Hamilton Harbour (Polak and Haffner, 1978) 
ky = 0.24 + 0.07¢ (g 05 m3 a7!) (8) 
the total oxygen used in the water of layer i is 
WCOD; = Vy(to-t)) (0.24+0.035(cj+c9)) (9) 
Combining’ (1), °¢€2);°(3),; (7) and’ (9) 
iMpnow (Ga sth) (Ava aGas) cee ko 
- Sj (ty - ty) (0.72 + 0.13(cy + c9)) 
=Ve "(tp = 't)) (0.24 + 0.035(e] + c5)) (10) 


From this equation, the vertical turbulent fluxes rj and rj4] 
can be calculated since they are the only unknowns, using the 
following procedure. 


For the bottom layer, there is no turbulent exchange downwards and 
rj, and R, therefore equal zero. From eq. (1), 


Ry = 4M)+SOD,+WCOD, Cet) 
or, expressed in terms of eq. (10), 
ry = {amy +8; (e2-t])(0.7240.13(ey+e9)) 
+V} (to-t))(0.244+0.035 (6140) b/g (tg-ty) (12) 
Knowing now ro and Rg, R3 and r3 can be calculated from 
R3 = A Mg+R2+SOD2+WCOD2 Cia) 


By iterating this process for the remaining layers, all the other r 
terms can be calculated. 


The vertical oxygen fluxes for different periods in 1976 were 
calculated (Fig. 5). In April, the fluxes were as high as log 09 
m~2d~! at the surface and 6.5 at 18 m. The fluxes decreased 
considerably to between 6.3 and 1.1 by June 20. The averaged fluxes 
for the whole field season (April 18-October 27) and the stratified 
period (June 11-September 16) are also shown. 


The flux is a function of both the eee radients and the 
coefficient of eddy diffusion D (in cm“ s~*) and can be 
expressed by 


F = 8.64 p 2 (14) 
dx 
where F is the oxygen flux through a unit area in gm72q71, Cue 
the oxygen concentration in mg Uae x is the distance in metres and 
8.64 a conversion factor. 


Coefficients D calculated from average oxygen fluxes and oxygen 
gradients for stratified periods of 1975 (June 17 to September 16), 
1976 (June 11 to September 16) and 1977 (June 9 to September 22) are 
shown in Figure 6 along with average depth profiles of temperature 
and dissolved oxygen concentrations for the stratified period of 
1976. The eddy diffusion coefficients displayed a marked minimum 
near the thermocline. This indicates that the turbulent exchange is 
at a minimum in this layer. However, oxygen flow through the 
thermocline was still considerable because of the steep oxygen 
gradients (Figure 5) and at no time during the year was the 
hypolimnion effectively isolated from the epilimnion. The vertical 
exchange along with the water column oxygen demand explains why even 
the epilimnetic waters were undersaturated for prolonged periods of 
time in the summer. The only occasions in which the mean dissolved 
oxygen concentrations exceeded saturation were at times when rapid 
heating was taking place. This effectively formed a more stable 
water column during this period, cutting off exchange between the 
epilimnion and hypolimnion. The occasions on which the dissolved 
oxygen depletion in the surface layer was most severe coincided with 
periods of higher vertical exchange, early in the season during the 
partial overturn in June and late in the season during fall overturn. 


During all three years, lenses of oxygen rich water were observed in 
the otherwise anoxic hypolimnetic waters of the harbour. Figure 6 
shows a marked local maximum at depths of 16 to 17 m. The lens of 
oxygen-rich water was most evident at station 4 during late June and 
July 1977. It is believed that these lenses represent a portion of 
oxygen-rich waters from Lake Ontario entering the harbour through the 
Burlington ship canal. 


PREDICTIVE MODEL 


Some of the results and equations developed previously for the SOD 
and WCOD in Hamilton Harbour (MOE, 1978; Polak and Haffner 1978) can 
be combined into a simple one layer box model capable of predicting 
how different water quality abatement measures could affect the 
oxygen situation of the harbour. 


If on one date the oxygen stock in the harbour is OS], then at some 
date n days later the stock OS2 will be given by 


0S2=0S1+R+LH+PHOT-WCOD-SOD Cis) 
with the individual terms described and estimated below. 


Assuming that the harbour volume V does not change between those two 
dates we can write 


0S1 


Ty (g 0) (16) 


OS2 = coV (g 09) (17) 


where ay and C9 are the average oxygen concentrations in mg Lot 
at the two dates. 


The amount of oxygen used by the water column between the two dates 
(WCOD) can be estimated using the equation 


WCOD = nV(0.24 + 0.07 (c) + co)/2) (g 02) (18) 
The sediment will consume 
SOD = nA (0.72 + 0.26 (cy; + c9)/2) (g 09) tg) 


Here A is the area of the harbour in m2 and V is the volume of the 
harbour in m”, 


The aeration rate for the early season was estimated (Polak and 
Haffner 1978) at 9.0 g 09 m2 d-! and therefore the oxygen 
supplied by aeration (R) is 


R= 9nA (g 0) (20) 


The lake-harbour net oxygen exchange (LH) at 1% of the harbour volume 
per day can be estimated by the equation 


LH = 0.01nV(12-(cy+c7)/2) (g 09) (21) 


assuming that the average lake water oxygen concentration for the 
period was 12 mg L7!. Photosynthesis (PHOT) was estimated (Polak 
and Haffner, 1978) to be 1.0 g 09 m~2 q-!, Its contribution is 


PHOT=nA_ = (g 09) (22) 


The extent of reduction in the WCOD can be expressed by (100-x)/100 
where x is the percent reduction of the present WCOD. Thus x = 0 
fixes the WCOD at the present level and x = 100 means complete 
elimination of the WCOD. Introducing equations (16) to (22) into 
(15) and rearranging results in 


Co =[afp-2aaevP 12-0.24(100-)/10e ¥en 0.00540,035(100-x)/100)-0. 1303] 
{¥+n(0.005V+0.13A+0.035V(100-x)/100)} (23) 


As an example of the use of this equation, concentrations for April 
23rd to June 18, 1976 were calculated using different values of x 

(Fig. 7). With no reduction in WCOD (x=0), the calculated average 
harbour dissolved oxygen concentration on June 19, 1976 was 2.6 mg LE. 
This is in close agreement with the average of 2.4 mg L* based on 
survey data. By reducing the WCOD by 20 percent, the calculated 
average oxygen concentration was 5.2 mg L ~. To reach an average 
concentration of 9.6 mg L7! of oxygen, which is the saturated 
dissolved oxygen concentration at the harbour average temperature, 

the WCOD must be decreased by 46 percent of the 1976 value. 


In interpreting Figure 7 it has to be borne in mind that the 
calculated values are only estimates because of all of the 
approximations used during development of equation (23). One of the 
most serious approximations is that the atmospheric reaeration, 
which is the main source of the oxygen to the harbour, remains 
constant. In reality, the extent of reaeration decreases with 
increasing surface water DO, and is also dependent upon the wind 
velocity. Thus the actual decrease would probably be represented 
more closely by a curve located below the one shown in Figure 7, and 
higher reductions in WCOD would be required to retain complete 
Saturation. The model assumes a one layer system while during the 
stratification period two layers are formed and oxygen depletion in 
the hypolimnion is speeded up by its partial separation from the 
oxygen-rich epilimnion. 


RELATIONSHIP OF SEDIMENT OXYGEN DEMAND TO OXYGEN CONCENTRATION AND 
TEMPERATURE 


It was shown (MOE, 1978) that sediment oxygen demand as measured in 
the laboratory depends on the temperature and on the oxygen 
concentrations in the overlying water. The dependence was given in 
a graphical form. It would be more convenient, however, to have a 
general mathematical function which will express this relationship. 
This would allow the use of a single measurement of the sediment 
oxygen demand (preferably measured in situ) to determine values for 
a wider concentration and temperature range. A simple mathematical 
expression for this dependence can be derived as follows. 


Assuming that the sediment oxygen demand is only a function of 
oxygen concentration c and temperature t we can write 


ks = kg(c,t) (24) 


by differentiation 


dk,= Oks dc + dkg\ dt 
Ocese dt fe (25) 


B=6 


The problem is to find the partial differentials (Qk,/dc)t and 
(dkg/dt),.. It was shown in the 1976 report that the dependence 
of the sediment oxygen demand on oxygen concentration in the 
over lying water at constant temperature can be expressed by the 
equation 


k, = acD (26) 


This equation was shown to be valid for diffusion accompanied by a 
chemical reaction (Edwards and Rolley, 1965). 


Equations expressing the dependence of the sediment oxygen demand on 
temperature at constant oxygen concentration have been suggested 
previously. Hargrave (1969) used an equation similar to equation 
(26). Generally, a classical Arrhenius type relationship can be 
used for the dependence of the reaction rates of both chemical and 
biological reactions on temperature. This can be written for k, as 


where R, S are constants, and T is the temperature in degrees K. 


Equations (22) and (23) can be used to express the partial 
differentials in equation (21). 


Assuming that a and b are independent of oxygen concentration, then 
from equation (26) 


(dk,/dc)_ = b(k,/c) (28) 


Similarly from equation (27) assuming R and S are independent 
of temperature 


(dk,/dt), = Sk,/T? (29) 
Introducing (28) and (29) into (25) we get 


dk, = b(k,/c)dc+S(k,/T*)dt 


s 
which gives after integration 


keg cb,-5/T (30) 


Equation (30) was used for correlation of the 1976 experimental 
laboratory data. Since it is nonlinear for parameters Q, b and S, 
the Newton-Raphson method with the least-square fit of the 

increments was used to evaluate their numerical values. For the 

1976 laboratory data the final equation is 

= 3.53 x 10/7 <0-346,-4989/T (31) 


Ks 


where kg is in g 09 m-2d-li¢ is in mg L7! and T is in degrees K. 


Equation (30) can be used to calculate the sediment oxygen demand at 
different temperatures and concentrations. Thus, k,(t2,¢9) 
can be calculated from k, (ty cy) using: 


a ee k,(ty,¢])(e9/e,)® exp {s(1/T,-1/T2)} 
if b and S are known. C325) 


During derivation of equation (30) some approximations were used. 
These include independence of coefficient b on concentration and 
coefficient S on temperature. As this is strictly not true, 
equation (30) is valid only for limited temperature and 
concentration intervals. Also, coefficients given in equation (31) 
were calculated for one type of sediment only and specific 
laboratory conditions. For different sediment and/or experimental 
conditions the coefficients will have different values. 


It was not possible to calculate the coefficients of equation (30) 
from the 1976 Hamilton Harbour field data, as the temperature varied 
very little from survey to survey. Data of Edwards and Rolley 
(1965) were fitted successfully. For the two different sediment 
types given in their paper, values of b were 0.432 and 0.383 while §S 
had values of 6488 and 7288, respectively. There was much larger 
variation in the value of Q. Variations in b and S were not too 
large considering that sediments were from different locations and 
origins (river mud affected by pollution), and agree well with the 
constants of equation (31). Thus, even though the absolute value of 
the sediment oxygen demand is different, the temperature and 
concentration dependence is similar and equation (32) can be used 
for interpolation with average b of 0.39 and S of 6250. Equation 
(32) can be rewritten in the form 


Kettoycon = kg (ty ,¢])(eg/c,)9-39 exp {6250(1/T)-1/T2)} (23) 
This general equation is believed to be suitable for temperature and 


concentration transformation of SOD data for sediments affected by 
pollution. 


CONCLUS IONS 


As 


The average oxygen concentration in the hypolimnion of Hamilton 
Harbour was below the provincial objective of 5 mg L-! between 
the middle of June and the end of September, 1977. This 
Objectives is designed for the protection of aquatic biota when 
the temperature is between 10° and 15°C. The average 
epilimnion concentration was always greater than 5 mg ish, 
reaching a minimum slightly below 6 mg L7l near the end of 
September. The epilimnion is, however, below saturation from 
the end of May to the middle of October. 


A greater simplified one-box model of the oxygen concentrations 
in the Harbour predicts that, to achieve average oxygen 
concentration greater than 5 mg L7!, the total water column 
oxygen demand should be reduced by around 20%. The 
relationships between the water column oxygen demand and the 
biochemical processes in the Harbor and between the demand and 
the loadings to the Harbor from industry, municipal sewage 
treatment plants and stormwater runoff have to be determined 
before definite recommendations for water quality management can 
be made. 


A very pronounced minimum in the calculated D, the vertical 
diffusion coefficient, was found at the thermocline. Because of 
the high oxygen gradients, the oxygen flux was still 
considerable and at no time during the summer was the 
hypolimnion effectively isolated from the epilimnion. 
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FIGURE 7 ESTIMATED FINAL CONCENTRATION ON JUNE 18, 1976 FOR DIFFER 
—ENT REDUCTION OF WATER COLUMN OXYGEN DEMANDS. 
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HAMILTON HARBOUR STUDY 


1977 


Volume I 


SECTION C 


WATER CHEMISTRY 


SUMMAR Y 


The water quality survey initiated in 1975 was continued during 1976 
and 1977 to assess the effectiveness of the artificial mixing 
program and extend the knowledge of harbour water quality trends. 

In addition, a survey of nutrient and organic chemical pollution at 
storm water outfalls was conducted in 1977. 


Dissolved oxygen results in both years showed the effect of water 
column instability, lake water intrusion and artificial mixing. 
Midsummer DO values next to the sediment were as high as 6 mg/L at 
station 4 and 3 mg/L at station 258. 


Ammonia values in spring 1977 were considerably higher than those 
observed in recent years with the mean value being 5 mg/L. 
According to outfall survey samples, the source of the higher 
anmonia appeared to be the Hamilton WPCP. The seasonal nitrate 
pattern was dominated by a strong nitrification peak at the end of 
June which coincided with the time of maximum harbour BOD. Maximum 
nitrate values (in July) averaged 2.6 mg/L. 


Total phosphorus values were slightly lower at 0.06 mg/L in 1976 
than in 1975. In 1977, dissolved phosphorus results were slightly 
higher at the east end of the harbour. Otherwise, there were no 
Significant changes in this parameter. 


Silica concentrations were again higher in the hypolimnion than in 
the epilimnion in both years. In August 1976, however, an unusually 
high hypolimnetic silica concentration occurred throughout the 
harbour. 


No significant seasonal or year-to-year trends in iron and zinc 
concentrations were found in the period 1975-77. Mean iron and zinc 
concentrations were 0.3 and 0.08 mg/L, respectively. This zinc 
concentration is high enough to cause a potential algal toxicity, 
but the extent of toxicity is dependent upon such factors as 
chemical speciation and interaction with other metals. 
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Determinations of organic matter in the harbour by chemical and 
biochemical oxygen demand and total organic carbon indicated that 
these parameters are not correlated with one another. This 
situation has been frequently observed in natural waters and occurs 
because the different tests measure significantly different 
substances. Each parameter varies significantly with sampling date; 
maximum BOD values were found in June at a time corresponding to the 
maximum nitrification rate. The observed correlation of dissolved 
organic carbon with BOD suggests that organic carbon plays an 
important role in the water column oxygen demand. Most of the 
industrial, municipal and storm water outfalls are sources of 
dissolved organic carbon. 


INTRODUCTION 


As part of the multi-year water quality study in Hamilton Harbour, 
an intensive survey program established in 1975 was continued in 
modified form during 1976 and 1977. The same four stations (4, 252, 
258, and 270; Figure 1) were used and three to five depths per 
station were sampled according to the total water depth. In 1976, 
the harbour was surveyed once every two weeks fron late April until 
the commencement of artificial mixing in June, and then 
approximately weekly until late August, and irregularly thereafter. 
In 1977, surveys were done every two weeks from late April to 
mid-October. Some monthly sampling runs for nutrients (8 stations 
in 1976 and 10 stations in 1977) were also done. In addition, 
detailed depth profiles of temperature, dissolved oxygen, redox 
potential, pH, and conductivity were obtained by NERA water quality 
monitoring, approximately weekly during both field seasons. During 
1977, a survey of stormwater runoff flow and chemistry was also 
conducted. This consisted of continuous conductivity, half-hourly 
temperature and dissolved oxygen, and hourly biological oxygen 
demand (BOD), chemical oxygen demand (COD), total organic carbon 
(TOC), and filtered organic carbon (FOC) measurements for six hours 
at the locations shown in Figure la, approximately three times per 
location. One sample for nutrients (N and P) was also taken during 
each survey. The purposes of the chemical survey were to extend the 
knowledge of water quality trends in Hamilton Harbour, to assess the 
dissolved oxygen deficit and measure the effects of abatement. 
procedures in operation. The purpose of the stormwater survey was 
to obtain accurate real-time histories of flows and water quality at 
the storm sewer outfalls, suitable for input to the numerical 
models, as well as improving our knowledge of the time variation of 


the harbour processes. 


Details on the artificial mixing procedure and the results obtained 
are covered in Section A, Artificial Mixing. Measurements of 
sediment and water column oxygen demand and their contributions to 
harbour dissolved oxygen stock has been covered in detail in section 
B, Oxygen Budget. The dependence of water nutrient and heavy metal 
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concentrations upon events at the sediment-water interface is 
discussed with respect to a series of proposed experiments on 
sediment extraction and sediment-water exchange in Section D, 
Sediment Chemistry. 


This section consists mainly of a water quality update in which the 
trends in previously measured parameters are extended to 1977 and 
comparisons of water quality data for the interval 1972-77 are 
discussed with reference to the above-mentioned sections, where 
appropriate. As in the 1975 report, it has been broken down into 
several subsections to facilitate discussion of the harbour quality 
with respect to several chemical parameters. In this section the 
effect of the control procedures has been included where 
appropriate. 


REDOX POTENTIAL AND DISSOLVED OXYGEN 


In Figure 2, the depth-time distributions of redox potential are 
plotted for 1976. As in 1975, all results indicated potentials much 
lower than expected for oxygenated surface waters (Hutchinson, 1957) 
with average values being about 200 mv and considerable fluctuations 
occurring during the period of artificial mixing. 


As discussed in the 1975 report (MOE, 1977), it is difficult to 
explain the redox potential results. Various explanations suggested 
include exposure of the redox probe to sulphide - containing 
sediments (which produces a platinum sulphide film resulting in low, 
drifting potentials), or various iron redox couples. Equipment 
malfunctions during 1977 resulted in no reliable redox data being 
obtained for that year. 


Isopleths of dissolved oxygen results obtained at the four major 
sampling stations are given in Figure 3 for 1976 and Figure 4 for 
1977. Periods of anoxic or near-anoxic conditions (DO <1 mg/L) are 
indicated by shaded areas. The data for both years show the effect 
of harbour instability as a series of oscillations in the depth of 
constant DO values during the summer months at a depth of 6 to 10 

m. These apparent oscillations are due to the internal wave 
Structure within the harbour caused by the effect of Lake Ontario as 
well as wind-induced oscillations. 
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In 1976, the aerator supplied flows close to the critical value of 
0.14 ms74 (300 cfm) for most of the period from mid-July to the 
end of August (section A, artificial mixing). In 1977, the aerator 
produced flows above 0.14 ms72 only for the last three weeks of 
July; average flows were 0.11 m>s~+ in August and 0.09 m3s~1 

in the period June 13-July 10. By comparison of these dates with 
the DO-depth data given for stations 4 and 258 in both years, it is 
apparent that oxygenation of the near-bottom waters, which reached a 
peak in July of both years, and was stronger at Station 4 than 258, 
is a combined effect of the artificial mixing with the harbour-lake 
exchange. Colder lake water entering the harbour sinks to the 
appropriate density level and mixes with harbour water, but produces 
a stronger effect near the ship canal. However, this phenomenon 
should occur for a greater portion of the summer and not be 
concentrated within a period of a few weeks. In addition, the 
efficiency of the mixing in August 1976 was as good as that in the 
last half of July, but DO levels remained below 1 mg/L throughout 
August at both stations 4 and 258. 


In 1977, hypolimnetic dissolved oxygen values were significantly 
higher than in 1976, except perhaps after mid-August, when 
near-anoxic results were observed at the two deep locations. The 
largest change was observed at station 4 during June (compare Figure 
4a with Figure 3a). As aerator flow before mid-June was negligible, 
this is a strong indication of the effect of harbour-lake exchange. 
In addition, the highest dissolved oxygen recorded in the 
hypolimnion of station 258 occurred on July 5, just before optimum 
aeration was achieved. DO values at this location dropped below 1 
mg/L late in July but increased slightly above 1 mg/L for part of 
August. It cannot be said how much of this effect was due to the 
aerator. 


Figure 5 shows a series of dissolved oxygen-depth profiles taken 
July 13, 1977 on a roughly east-west line between the canal and 
station 258. This shows mixing of lake and harbour water due to 
lake water intrusion as far west as station 254, with apparently two 
lenses of higher DO water at station 254. Lower dissolved oxygen 
near the bottom at station 269 may have been due to short-term 
harbour instability; profiles taken at this location on other dates 


in 1975-77 have sometimes shown bottom DO values of greater than 8 
mg/L. The effect of lake water on the vicinity of station 4 was 
also shown by the low pollutant concentrations in the sediments and 
similarity to stations 269 and 1030 observed in 1977, as well as the 
mixed conditions at station 255 seen in several sediment sampling 
surveys (see section D, Sediment Chemistry). The effect of lake 
water intrusion and mixing on the harbour DO is to be studied 
further in 1978 using the three-dimensional numerical model. 


Experimental studies of sediment-oxygen demand and computation of 
the harbour DO stock has been discussed elsewhere (Section D, Oxygen 
Budget) and will not be repeated here. 


NITROGEN 


Data for 1976 and 1977 can be compared with results of previous 
years using techniques similar to those used in previous reports. 
Tables 1 and 2 show June-October surface averages for ammonia and 
nitrate at several stations as well as monthly surface means for the 
entire harbour. The June-October means were not computed in 1977 as 
a different set of sampling stations was used. Figures 6 to 16 show 
trends for anmonia, total Kjeldahl] nitrogen, nitrite-N and nitrate-N 
displayed for: (a) seasonal trends of harbour-wide surface means 
for 1972-77; (b) 1976 isopleths; and (c) 1977 isopleths. 


The most striking change in any nitrogen parameter is the dramatic 
increase in spring NH 4-N values in 1977, compared to 1975-76 

(Figure 6). Peak average values approached 5 mg/L, with station 4 
achieving a maximum of 6.3 mg/L on May 11. Observed results are 
comparable to 1972 figures, though the May 1972 peak was not sampled 
and may have been higher. A 10-station survey on June 9, 1977, 
yielded a result of 5.7 mg/L at station 20 in the southeast part of 
the harbour, 3.55 mg/L at station 4 and 3.55 mg/L at the other 
eight locations sampled. Two other surveys in 1977 also showed much 
higher concentrations at station 20. Although this was also true in 
1976, the 1977 mean (3.65 mg/L) was far higher than the 1976 mean 


(2.35 mg/L); it appears as though more NH is being discharged 
from the Hamilton sewage treatment plant or other location in this 
area. OQOutfall surveys in 1977 (Table 3) indicated the worst NH 
contamination to come from location 274 (Hamilton sewage treatment 
plant), with severe contamination originating also at other 
southeastern locations (220, Strathearne Street slip; 211, ditch 
carrying industrial effluent to Kenilworth Street slip; and 225, 
Parkdale Street). Values at all of these locations ranged upwards 
of 5 mg/L NH3-N to a maximum of 18 mg/L NH3-N (Table 3). 


Following the spring peak, the surface and bottom NH. 

concentrations decreased rapidly, such that July concentrations were 
comparable in each of 1975-77. In August, surface NH. 
concentrations were slightly higher in 1976 and 1977 compared to 
1975. During both years, epilimnetic ammonia concentrations were 
higher than hypolimnetic in the summer period (June-August), 
reflecting the surface nature of the discharges. Surface NHN 
concentrations are still considerably higher than the public surface 
water supply criterion of 0.5 mg/L for most of the season. 


The seasonal and depth variations of total Kjeldahl nitrogen 
(figures 9-11) parallel those of ammonia. This reflects the fact 
that the organic N content (TKN less NH.-N) of the harbour is 
roughly constant with time, exhibiting only irregular seasonal 
fluctuations. Most, but not all samples obtained in July-August 
1976 and June-early August 1977 exhibited organic N concentrations 
of around 1 mg/L while average concentrations before and after this 
period were 0.4 to 0.7 mg/L in 1976, and 0.6 to 0.8 mg/L in 1977. 
This tendency toward higher summer organic N values does not 
correlate with other parameters such as NO4-N, chlorophyll a, or 
total organic carbon. Organic N concentrations in both years were 
somewhat smaller than 1975 averages of 1-1.5 mg/L. 


The highest organic N concentration observed in 1977 at stormwater 
outfalls was about 10 mg/L. This value was obtained at locations 
211 (drainage ditch to Kenilworth Street slip) and 225 (Parkdale 
Street outfall). Both of these sewers are contaminated with 


industrial wastes. 
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Nitrite concentrations did not exhibit the extreme values picked up 
in 1975 (harbour-wide mean 1.1 mg/L NO.,-N on June 4, 1975). On 
June 14, 1976 and July 13, 1977, epilimnetic average concentrations 
of 0.35 mg/L NO.-N were recorded; concentrations of 0.25 to 0.3 
mg/L NO.-N were also observed in late August to early September of 
both years, while NO.-N concentrations of most other samples were 
below 0.2 mg/L. Epilimnetic concentrations tended to be slightly 
higher than hypolimnetic concentrations during most of the summer 
period. The variability of this parameter was too irregular to 
allow correlation with the operation of the aerator or any other 


chemical parameters. 


The general form of the seasonal surface nitrate curve (Figure 14) 
was similar in each of the years 1975-77, although the nitrification 
peak in 1976 occurred several weeks earlier than in 1975 or 1977, 
and the 1977 nitrification maximum was higher than any year since 
1970 (see Table 2). The nitrification peak in both 1976 and 1977 
coincided with the time of maximum BOD, on a harbour-wide average 
basis (Figure 24). In 1977, the nitrification peak is also related 
to the high spring ammonia values, although only a portion of the 
anmonia-N peak is nitrified. Total inorganic N, or (NH. + NO. + 

NO 3) - N gradually decreased from a mid-May maximum indicating a 
net biological uptake of N or denitrification. 


Combined seasonal-depth variations of NO 2-N in 1976 and 1977 
(figures 15-16) are less clear than they were in 1975, when a period 
of nitrification or sediment release during June was followed by 
denitrification under late summer hypolimnetic anoxia. In early 
June of 1976, nitrification or N release at the sediment-water 
interface allowed a build-up of NO..-N in the hypolimnion of the 

deep water stations. With the onset of anoxic conditions an 
irregular decrease in hypolimnetic NO 2-N occurred, with a minimum 
value being reached July 20. This was close to the time at which 
efficient ( 0.14 mos~l) aerator operation was established. 

However, this decline was interrupted by a period of NO 2-N 

increase between July 6 and 13. This was during the period in which 


dissolved oxygen was being supplied to the lower hypolimnion by 
artificial mixing and lake-harbour exchange as discussed earlier. 
During early August, a gradual release or build-up of bottom NO.3-N 
occurred, apparently connected with aerator operation; however, 
epilimnetic nitrate concentrations remained higher than hypolimnetic 
concentrations for most of the summer. After late August, 
denitrification occurred at the sediment-water interface, under 
near-anoxic conditions. The trend of NO2-N concentrations in the 
deepest waters, however, seemed to parallel the existence of 
oxidized conditions at the sediment-water interface. 


In 1977, however, the seasonal-depth variations of NO..-N are 

harder to explain. Throughout most of the summer, epilimnetic 
NO,-N concentrations were higher than hypolimnetic 

concentrations. Nitrification at the bottom appeared to occur only 
around the end of June, when rapid nitrification of ammonia was 
occurring throughout the water column. This occurred before 
efficient operation of the artificial mixing device had been 
achieved, but while considerable DO was still present in the 
hypolimnion at the east end of the harbour (figure 4a). Before the 
end of June, despite the continued presence of DO next to the 
sediment at least at station 4, irregular behavior occurred; whereas 
one might have expected a heavy NO. build-up during this period 
leading to a condition of hypolimnetic NO.2-N > epilimnetic 

NO 2-N. The reason for the lack of the expected trend is unclear, 
unless it is related to bacterial or phytoplankton populations. 
After mid-July, a tendency to decreased NO.2-N concentrations at 

the bottom occurred, despite aerator operation and DO values between 
1 and 2 mg/L. This is again in contrast to 1976 results. Such a 
decrease in NO 3-N concentration may be caused by nitrate 
assimilation by phytoplankton (Hutchinson, 1957). Alternately, 
denitrification within the sediments may be responsible (and thus 
indicative of anoxia within a few mm of the sediment-water 
interface). This is, however, at variance with 1976 observations; 
this difference cannot be explained at this time. 
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Nitrate values obtained at storm water outfalls (Table 3) were 
strongly affected by dissolved oxygen, owing to the high BOD and COD 
values generally exhibited. The highest NO concentration was 8 
mg/L in Red Hill Creek. Under DO values of less than 3 mg/L, 
nitrate concentrations were almost always below 1 mg/L. 


PHOSPHORUS 


Tables 4 and 5 give June to October surface means of total and 
soluble phosphorus for selected stations, and harbour-wide monthly 
surface means observed since 1967. Figures 17 and 18 present 
depth-time distributions of total phosphorus for 1976 and 1977, 
respectively, while Figures 19 and 20 present the same for filtered 


reactive phosphorus. 


Year-to-year trends are difficult to discern from the means of total 
and soluble phosphorus as presented in Tables 4 and 5. 
Consequently, the annual surface means were compared from one year 
to the next for each station using the Student's t-test. The 
results are given in Table 5a. A significant decrease in total 
phosphorus occurred at three of the four locations from 1975 to 
1976. Although the 1977 data suggested an increase from 1976, 
particularly in the 8- to 10- station survey data, the change was 
not statistically significant for total phosphorus. A significant 
increase in filtered reactive P did occur at stations 4 and 252 
(east end of the harbour) from 1976 to 1977. It must be remembered 
that soluble phosphorus is affected by large analytical variability 
as well as possible changes between sampling and analysis. In 
addition, only three surveys were performed in 1977 on the 
10-station grid, and one of these (August 10/77) coincided with a 
4-station survey which produced the highest average soluble 
phosphorus value of the entire season. These factors prevent one 
from drawing any definite conclusion about changes in phosphorus 
content in Hamilton Harbour from the available data. 


The worst areas of storm and sanitary sewer outfall contamination 
with regard to phosphorus are shown in Table 6. These areas are 
location 225 (Parkdale Street), 211 (sewer draining to Kenilworth 
Street slip carrying some industrial wastes), 273 (near Randle's 
Reef) and 274 (Red Hill Creek; Hamilton STP). Two of the five 
samples taken from location 274 exceeded the IJC guideline for 
maximum municipal sewage effluent total P of 1.0 mg/L. 


No seasonal trend is exhibited in the 1976 or 1977 total and soluble 
P data. A tendency for epilimnetic total phosphorus values to be 
Slightly higher than hypolimnetic values is observed at station 4 in 
both sunmers and stations 258 and 270 for parts of 1977. These 
changes are not noticeable in soluble reactive P. 


Both parameters showed occasional pockets of high results, which 
sometimes occurred on the same date (e.g. August 10/77) at several 
locations. This phenomenon has been observed previously (MOE, 1974, 
1977) but has not been explained adequately. 


Little, if any evidence is available for release of total or soluble 
phosphorus from the sediments in either year. Any slightly higher 
soluble P values that do occur are not correlated to either 
decreased DO or increased dissolved iron levels. Piccinin (1977) 
found a very high (but unstated) level of soluble reactive 
phosphorus at station 4 near the end of August 1976. Such a peak 
may easily have occurred between successive MOE sampling dates and 
have been missed; certainly there is no evidence here to confirm or 
deny this peak. Piccinin (1977) indicated the hypolimnion to be 
anoxic at the end of August; MOE data indicated about 3 mg/L DO at 
the bottom on August 23, but the DO monitor failed to function 
properly on August 30 or September 7; consequently, no proof of 
hypolimnetic anoxia at the end of August is available. 


Control of phosphorus concentration at the sediment-water interface, 
and proposed experiments to study this factor in more depth, are 
discussed in Section D, Sediment Chemistry. 


SILICA 

The seasonal variation of dissolved reactive silica (as Si) is shown 
as surface and bottom means in Figures 21 and 22 for 1976 and 1977, 
respectively, and in isopleth form for 1977 in Figure 23. Although 
a general similarity is apparent from 1975 (MOE, 1977) to 1977, 
differences are evident in each annual curve. 


In spring 1976, Si concentration was 0.6-0.8 mg/L, far lower than 
the 1.0-1.5 mg/L observed in 1975; consequently, the early summer 
minimum was reached earlier (mid-June) and with lower 
concentrations. With development of thermal stratification, 
epilimnetic Si concentrations become lower than hypolimnetic, as in 
1975; in July the concentrations again become homogeneous, 
presumably in response to artificial mixing, again as in 1975. 
However, the most homogeneous results occurred before optimum air 
flow was achieved, but at a time when hypolimnetic DO values were 
highest (see discussion on DO). In August 1976, irregular changes 
in Si concentration developed, along with a considerable increase in 
bottom Si. On August 17, large variations in Si content occurred 
across the harbour; the bottom Si concentration reached 1.1 mg/L at 
station 258 while surface Si concentration was 1.15 mg/L at station 
4 and 0.10-0.15 mg/L at other points. By contrast, the Si 
concentration at 3 m depth at station 4 was 0.25 mg/L; below this 
depth, the typical stratification (hypolimnetic Si  epilimnetic Si) 
had reappeared. This stratification persisted, with declining 
average values, until fall turnover. Although the bottom water at 
station 258 was near-anoxic at the appearance of high Si, the Fe 
concentration was negligible, suggesting that the high Si is not due 
to reduction of silicates at the sediment-water interface. Perhaps 
another metal is involved; or massive decomposition of diatoms could 
have released dissolved Si to the water column at this time. 
Hutchinson (1957, p. 794-5) suggests that silica can enter water 
from the sediments without complete anoxia. He suggests that 
diffusion of silica may be under control of temperature, but does 
not state why this should be so. 


The 1977 Si curves resemble the 1975 data more closely than the 1976 
period; however, the stratification did not develop until late June, 
Suggesting a delay in maximum diatom growth. In addition, the 
apparent homogeneous Si concentration at the end of July was not due 
to artificial mixing but was caused by a relatively high (0.45 mg/L) 
value at station 4. (See also Figure 23.) Although Si 
concentrations were somewhat lower in July-August 1977 then in the 
same part of 1975 and 1976, weak depth stratification did persist 
throughout the harbour in this period as Figure 23 indicates. 
Possibly a larger crop of diatoms may be responsible. These may 
have decomposed during autumn cooling and turnover, producing the 
relatively larger increase in dissolved Si concentration seen in the 
September-October 1977 period. 


HEAVY METALS 


In 1976 and 1977, iron and zinc were determined from all sampling 
points at every survey. The purpose of determining zinc was to 
assess the extent of toxicity to the harbour algae. Piccinin (1977) 
determined that zinc has a toxic effect on the alga Scenedesmus 
quadricauda at concentrations above 0.01 mg/L, although she states 
that algae cultured from Hamilton Harbour show an increased 
tolerance to zinc toxicity and that response of various species to 
zinc concentrations below 1 mg/L is inconsistent. 


As the zinc and iron concentrations showed no consistent spatial and 
depth variations and only irregular seasonal variations, 
harbour-wide monthly means were calculated. These data are given in 
Table 7. Mean zinc values range from 0.02 to 0.14 mg/L, with 
maximum individual results being 0.45 mg/L in 1976 (station 258, 
July 6) and 0.46 mg/L in 1977 (station 270, June 29). Mean iron 
values are about 0.3 mg/L, with occasional high values above 1.0 
mg/L observed in the epilimnion at station 4. No significant 
year-to-year trends are evident in either metal. No evidence of 
release of iron from the sediment-water interface was noted; the 
highest bottom Fe values observed in 1977 occurred at station 4 on 
June 15 and July 27 (0.7 mg/L) while the bottom dissolved oxygen was 
around 4-5 mg/L. By contrast, in August 1977 when the bottom DO was 
1 mg/L or less at station 4, the iron concentration was 0.3 to 0.5 


mg/L. 


A few samples were analyzed for manganese. Most results were below 
0.1 mg/L in the spring and fall; in September 1976 a hypolimnetic 
increase was observed at stations 4 and 258, where values of 0.2 to 
0.5 mg/L were found. This result is reasonable because manganese is 
released under less severe reducing conditions than iron. The redox 
potential to obtain 107° M (about 0.6 mg/L) concentrations of Mn 

and Fe corresponds to a pE of about 8 for Mn and 2 for Fe at a pH of 
near 7 (Stumm and Morgan, 1970, p. 533). 


The iron and zinc values may be compared to the Great Lakes water 
quality objectives (IJC, 1977), which state that iron should be 
below 0.3 mg/L and zinc below 0.03 mg/L. It is seen that zinc 
concentrations are greater than objective values most of the time, 
and that average iron concentrations just barely meet the IJC 
objective. That this is true is not surprising since the harbour 
sediments are known to be highly polluted with respect to these 
metals (Section D, sediment chemistry; MOE, 1977, Section C). 


Water samples were analyzed on one occasion only (May 11, 1977) for 
mercury. All samples were below 0.03 ug/L, indicating no problem in 
the water column for this metal. On the same date, aqueous copper 
was below 0.01 mg/L. 


Environmental aqueous heavy metal concentrations are strongly 
regulated by events occurring at the sediment-water interface (for 
example, see Hutchinson, 1957; Stunm and Morgan, 1970, and 
references therein). Relationships between iron and phosphorus 
within the water column and at the interface in Hamilton Harbour 
have already been explored (MOE, 1974, Section D). Studies of 
sediment samples obtained in Hamilton Harbour in the period 1975-77 
have suggested that control of heavy metal and phosphorus release 
from the sediment may be more complicated than the simple 
iron-phos phorus-dissolved oxygen system (see Section D, Sediment 
Chemistry, and MOE, 1978). 


Although most heavy metals become less soluble under oxidized 
conditions (Stumm and Morgan, 1970), zinc may become more soluble at 
higher redox potentials (Hem, 1972). Equilibria involving zinc 
carbonate or hydroxide in water of composition similar to Hamilton 
Harbour predict a zinc concentration of about 1076 M (0.65 mg/L) 
at a pH of between 8 and 8.5. Silicate equilibria predict a 
concentration of 1077 M (0.06 mg/L) at a pH of 8 and a silica 
concentration of 107" M (2.8 mg/L as Si). This zinc concentration 
is approximately that observed in the harbour water. Equilibrium 
concentrations decrease at higher pH. Under reducing conditions, 
ZnS precipitates over a very wide pH range. Examples given by Hem 
(1972) tend to support silicate control of zinc solubility in many 
fresh waters, and it is possible that observed zinc concentrations 
in Hamilton Harbour are controlled by this equilibrium. If humic 
Substances are present at a concentration of 5 mg/L or more, 
increased zinc solubility as humic complexes will occur (Wilson, 
1978). This may occur at the dissolved organic C concentrations 
present in Hamilton Harbour, but these interactions are highly 
dependent on the nature of the organic material present. 


Chemical speciation of toxic heavy metals in the water column is 
important in controlling their solubility and toxicity to biota. An 
initial study of heavy metal speciation in Hamilton Harbour was 
performed by Chau and Lum-Shue-Chen (1974), who analyzed "labile" 
(ionic and rapidly exchangeable complexed) and strongly bound 
concentrations of Zn, Cd, Pb and Cu in a Hamilton Harbour sample. 
They obtained concentrations in ug/L for strongly bound metals as 
follows: Zn, 27; Cd, 2.4; Pb, 1.3; and Cu, 49. The only detectable 
labile metal was Zn, at a concentration of 7 ug/L. Mancy and Allen 
(1977) reviewed the effect of organic complexing agents on heavy 
metal toxicity in detail, with particular reference to copper. 
Although some disagreement exists in the literature, they observed 
that in general organic chelating agents such as polypeptides or 
humic acids considerably reduce the toxic effect of heavy metals. 


An additional factor of importance in toxicity studies is the 
interaction between different heavy metals. At times, the full 
toxic effect of one metal may not be observed in the presence of a 
second metal (United Nations, 1976). This was indicated for the 
pairs cadmium-zinc, lead-calcium, and mercury-selenium. The effect 
of cadmium may thus be responsible for the increased tolerance of S. 
quadricauda to zinc toxicity and inconsistent responses mentioned 
above. On the other hand, increased toxic effects have been 
observed in the combination mercury-lead-zinc (United Nations, 
1976). A recent study of heavy metal toxicity to S. quadricauda 
(Wong, Chau and Luxon, 1978) revealed a synergistic toxic effect of 
a mixture of ten metals each present at the IJC (1977) recommended 
maximum concentration. The same study showed an additional toxic 
effect of Hamilton Harbour water on S. guadricauda compared to a 
test growth medium. Certainly the mechanisms involved in heavy 
metal algal toxicity are complex. 


Further studies on aqueous heavy metal concentrations in Hamilton 
Harbour would have to include determining metal speciation. This 
involves complicated electrochemical techniques such as anodic 
Stripping voltammetry (Chau and Lum-Shue-Chen, 1974; Mancy and 
Allen, 1977) which are at present beyond Ministry laboratory 
capabilities. However, studies in this field would correlate well - 
with proposed and existing studies on sediments and heavy metal 
toxicity. 


ORGANIC POLLUTION AND OXYGEN DEMAND INDICATORS 


In 1976, surface and bottom water samples from each location were 
analyzed for biochemical and chemical oxygen demand (BOD and COD), 
and total and filtered organic carbon (TOC and FOC). In 1977, 
samples from all depths and locations regularly used were analyzed 
for these parameters. In addition, samples from the stormwater 
outfall locations (Figure 1A) were collected hourly for 6-hour 
periods on each of two to four days during 1977, in an effort to 
measure the magnitude and time variation of oxygen demanding 
loadings to the harbour. 


To assess the extent of variability of these parameters and their 
relationships with one another, statistical methods of correlation 
analysis and analysis of variance were used, with emphasis on the 
1977 data due to the larger amount of information available. 
Three-way analysis of variance (depth-station-date) indicated that 
BOD, COD, and TOC all varied significantly at the 0.1% level with 
respect to sampling date. TOC exhibited a significant variation at 
the 1.0% level with respect to depth. Variation of other parameters 
with respect to depth, and of all parameters with respect to 
sampling station, was not significant. 


These parameters were next correlated with one another on the basis 
of harbour-wide means for each sampling date. As the results were 
highly variable and frequently indicated negative correlations, it 
was decided to correlate all 1977 data together. Turbidity and 
suspended solids were included in order to ascertain any 
relationships with particulate matter. FOC was not included at this 
stage aS analyses for it were not done on all 1977 sampling dates. 
The results, given in Table 8, indicate no significant correlations 
between any of these parameters, except for the expected turbidity- 
suspended solids relationship. The lack of relationships between 
TOC and turbidity or suspended solids is not surprising in view of 
the fact that FOC made up over 90% of TOC in most samples tested. 


To further explore the significant variations with sampling date and 
the lack of observed inter-correlations, and facilitate inter-year 
comparisons, harbour-wide means of each parameter were plotted 
against sampling dates for 1976 and 1977. These data are given in 
Figures 24 (BOD), 25 (COD), and 26 (TOC). Error bars are +1 
standard deviation. TOC data were plotted in 1977 as surface and 
bottom means due to the observed significant variation with depth; 
bottom samples were displaced slightly on the time axis for clarity. 


Figure 24 shows that in both years BOD exhibits a similar seasonal 
dependence, increasing to a June maximum followed by a decrease to a 
relatively constant value of 2-3 mg/L by the end of July. The 1977 
peak mean (6 mg/L) was slightly less than the 1976 peak mean (8 
mg/L). No significant changes occur for the rest of the season. It 
is interesting that the maximum BOD occurs at a time of the most 
rapid nitrification in the harbour, as well as when the harbour 
dissolved oxygen stock is decreasing rapidly. Normally, the 5-day 
BOD test measures largely carbonaceous biochemical oxygen demand; 
nitrogenous BOD is generally exerted at longer time intervals. In 
addition, nitrogenous BOD is not often well measured during the 
common BOD test as the nitrifying bacteria may not be well 
acclimatized to the conditions of the test (D. Weatherbe, personal 
comm.). Studies in the Grand River basin (MOE, 1976) showed that 
nitrogenous BOD made up a significant portion of the total BOD of 
sewage effluents. Additionally, an effect of nitrification on the 
measured BOD. values was observed by Chandler et al (1976) in a 
Massachusetts River. On the other hand, Hamilton Harbour 
Nitrosomonas populations in 1976 showed only an erratic seasonal 
dependence with no relationship to the observed BOD peak. The 
reason for this lack of correlation is unknown. 


Unlike BOD, the seasonal dependence and year-to-year consistency of 
COD and TOC is less evident. COD shows very little seasonal changes 
in either year, with the significant variation in the 1977 analysis 
of variance likely related to the high values’ observed on June 29 
and August 24. The former coincides with the BOD peak, while BOD 
was not measured on the latter date. The 1977 annual mean COD value 
(27 mg/L) was slightly higher than the 1976 mean (21 mg/L); unlike 
the BOD situation. TOC shows strange seasonal variations in both 
years, which do not relate to one another or to the other 
parameters. Surface TOC values in 1977 are slightly higher than 
bottom water values, as expected, considering the epilimnetic nature 
of the discharges to the harbour. No apparent relationship to the 
artificial mixing program is evident. 


The lack of correlation between BOD, COD and TOC in natural water 
samples (as opposed to effluents) has been generally noted in the 
literature (for example, Chandler et al, 1976; Jones, 1972). This 
is not surprising as some organic chemicals resist biochemical 
oxidation and others resist chemical oxidation, but none resist the 
catalytic oxidation procedure of the TOC analysis. In addition, 
some inorganic materials such as ferrous iron and sulphide are 
oxidized during the COD test. The BOD test is also highly dependent 
upon the nature of the seed organisms used in the test; if these are 
unlike the bacteria naturally found in the water, the result cannot 
meaningfully relate to oxygen uptake data. 


A weak but significant relationship (r=0.215) was obtained upon 
correlating BOD with FOC for all samples on the eight dates in 1977 
in which these parameters were both measured. This supports the 
Suggestion (G. Harris, personal communication) that dissolved 
organic carbon plays an important role in the water column oxygen 
demand. However, 5-day BOD undoubtedly represents only a small 
portion of available oxygen demand within the harbour, especially 
considering that the average water retention time (including 
lake-harbour exchange, industrial intakes and discharges) is about 
100 days. Nevertheless, the dissolved organic carbon content does 
present a problem, as it represents the ideal substrate for 
heterotrophic bacteria which are probably controlling much of the 
water column oxygen demand (G. Harris, personal communication). 
Some possible sources of this FOC are the Ottawa and Kenilworth 
Street slips, the Hamilton WPCP, and outfall locations 273, 211, 
220, and 225. More study is desirable in locating the sources of 
this FOC. 


In Table 9, the 1976 and 1977 means of all BOD, COD, TOC and FOC 
data obtained within the harbour are compared to mean values of the 
same parameters obtained at the various storm and sanitary sewage 
outfall locations. Additional details on the outfall results will 
be presented in a separate report under preparation; in this 
context, it should be mentioned that concentrations observed at the 
various outfalls varied widely from day to day and even within a 


given day. Standard deviations of most outfall data in Table 9 are 
close to or greater than the means. Table 9 indicates that the most 
concentrated effluents entered the harbour at locations 273, 211, 
225 and 274. Location 274 is Red Hill Creek, which contains the 
Hamilton WPCP effluent, while the other three locations all contain 
industrial wastes. Due to the large variability observed in the 
remainder of the data, none of the other means presented in Table 9 
are significantly greater than the 1977 mean values obtained within 
the harbour as indicated on the second line of Table 9. Except for 
locations 250 and 212 (Ottawa and Kenilworth Street slips), which 
have high industrial effluent flow rates and perhaps 220 (also 
contaminated with industrial effluent), these locations do not 
present a problem under dry conditions. Although very few results 
were obtained under conditions of stormwater runoff, results that 
did occur under these conditions indicated possible problems. For 
example, at location 108, maximum BOD, COD, and TOC concentrations 
observed shortly after commencement of rainfall were 32, 207 and 56 
mg/L; similar concentrations persisted one hour later, but 
concentrations had dropped considerably after two hours. 


MINERAL CHEMISTRY AND CONDUCTIVITY 


In the report on 1975 results (MOE, 1977), a model was described in 
which the conductivity of a water sample is calculated from the 
equivalent conductances of the component ions and ionic interaction 
effects (Rossum, 1975). This model was verified for single salt 
solutions and applied to Hamilton Harbour water samples. It showed 
that non-conservative (biodegradable) materials contributed 2-3% to 
the observed conductivity and that chloride and sulphate contribute 
an average of 36% to the conductivity, as compared to 29% for Lake 
Ontario and 22% for world mean river water. 


Surface and bottom water samples at the four major locations were 
analyzed once in 1976 and twice in 1977 for chloride, sulphate, 
sodium, potassium, calcium and magnesium. In addition, nitrate, 
anmonia, pH and alkalinity were measured during each survey; 
carbonate and bicarbonate concentrations were calculated from pH and 
alkalinity using the second ionization constant of carbonic acid 
(MOE, 1977). 


Observed concentrations and computed contributions to the 
conductivity of each sample are given in Table 10. Unlike the 1975 
Situation, where the calculated conductivity values were higher than 
the observed values by 6 to 12%, the calculated conductivities were 
within 5% of the observed values for all but one sample. These 
variations are better than the coefficients of variation (6-9%) 
observed in an interlaboratory comparison of Lake Ontario and 
Grindstone Creek conductivity samples (McGirr, 1974). It did not 
appear necessary to consider ion-pair formation as affecting 
conductivity, as was done in 1975 (MOE, 1977). Although this 
appears to indicate reasonable analytical precision, a deviation of 
the ionic balance is observed in the autumn samples of both years. 
Equivalent anionic concentrations were generally 5 to 8% higher than 
cationic concentrations. This indicates a systematic deviation in 
the analysis of one or more chemical parameters. That errors of 
this magnitude are not unusual was noted in an interlaboratory 
comparison study by Ekedah] and Rondell (1973), who noted that 
Systematic errors are more important than random errors. Deviations 
of the ionic balance in May 1977 were in random direction, and less 
than 6%. 


A larger deviation was indicated for the surface sample at station 
270 on October 13, 1977, where the ionic balance deviation was 13% 
and the conductivity deviation was 7%. Examination of the data 
indicates that the sulphate result is probably too high; a value of 
about 50 mg/L, similar to that observed with other samples taken on 
the same date, would produce better agreement. This shows the 
usefulness of these calculations in checking analytical data and 
Suggesting a parameter which may need re-analysis. 


Ranges and seasonal variations of major ion concentrations are 
similar in 1976-77 to those observed in 1975, except for the 
increased anmonia concentrations in spring 1977 (discussed under 
Nitrogen) which appear to be balanced by slightly increased 
bicarbonate. Calcium and bicarbonate were also very slightly higher 


in fall 1977 compared to fall 1976. The contributions of chloride 
and sulphate to the conductivity of Hamilton Harbour water averaged 
36% in November 1976 and May 1977, and 34% in October 1977. These 
figures are similar to 1975 results and are as expected, considering 
the major use of Hamilton Harbour as a receiving water for 
industrial and municipal wastes. If the non-conservatives are 
represented by the N parameters (ammonia and nitrate), the results 
in Table 10 indicate that non-conservative materials contributed 2 
to 5% of the total conductivity in 1977, the increased range being 
due to high spring ammonia N levels. 
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TABLE 1 
Ammonia Surface Data, 1966-1977 


(a) June-October averages for each station 
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Note: All figures are in mg/L as N 
Analytical precision + about 10% 
Standard deviations of above figures (including station 
and time variations) are generally 30--50% of the means. 
*4-station weekly surveys (biweekly in 1977) 
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TABLE 2 
Nitrate Surface Data, 1966-1977 


(a) June-October averages for each station 
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Analytical Precision + about 10% 

Standard deviations of above figures (including station 
and date variations) are generally 30-50% of station means 
and 10-50% of monthly means. 

* 4-station weekly surveys (biweekly in 1977) 
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TABLE 3 


Outfall Areas with Highest Nitrogen Concentration, 


NH3 and TKN 


Location 


274 
211 
225 
220 
2713 


NO 3 


Location 


274 
220 


212 (surface and bottom) 


Hamilton Harbour, 1977 


Typical Values, mg/l as N 
TKN 


NH3 

11-18 14-38 
4.7-12 12-32 
7.7-11.2 16-24 
2.1-15 6.3-19 
3.6-6.4 6.6-9.6 


Typical Values, mg/l_as N 


up to 8 
up to 3.4 
up to 3.4 


TABLE 4 
Total Phosphorus Surface Data, 1967-1977 


(a) June-October averages for each station 


Station 

20, G25Ge 94257 auleeSSty yha262 265 (270) 269 
1967 pile .091 .087 .067 077 . 103 126 
1968 .099 -060 .034 .035 .050 .092 .048 
1969 .093 084 044 043 049 .061 .067 
1970 TCL 5** pecUuet .036 044 .054 .046 .065 
1971 wah) «055 055 F031 .037 +037 .040 
1972 a ay -086 .072 .068 .065 .068 .078 
1974 -088 Oy O52 .057 .065 .066 - 
1975 viZ2 .072 -082 .054 -066 .058 .059 
1976 -073 S077 .061 .046 .048 .069 062 


(b) Monthly and annual averages for entire bay 


dun Jul Aug Sep Oct Period 

1967 BO As COG) <werL07 ay) 102 .095 
1968 048 OOUrn 00S .075 .068 .060 
1969 .058 048 .037 .087 .094 .063 
1970 sogot* 091 2065 .078 E059 .078 
1971 091 sQooF e072 .098 - .073 
1972 - Oy AES olbley 088 .093 .088 
1974 - C007, 2102 ~ s0517 .084 
1979* 082 086 .069 077 .063 .075 

# .084 074 = .082 .050 .069 abies 
1976* 064 2061" 5053 .060 044 .059 

# .077 054 - - .075 stig 
Lo S073 068 .063 5053 .060 .065 

# .098 -071 » 087 - - .085 
Note: All Figures are in mg/L as P 


**One questionably high value was onitted 

Analytical Precision + 0.003 mg/L or 10%, whichever is 
greater. 

Standard deviations of above figures (including station 
and date variations) are generally 50-100% of the means. 
* 4-station weekly surveys 

# 8-station monthly surveys (10 stations for 1977) 
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TABLE 5 
Soluble Phosphorus Surface Data, 1967-1977 


(a) June-October averages for each station 


Station 

20 2G * <25Ajeq 258 Wen 262 M3725 (270)... 269 
1967 .060 -060 .058 056 114 . 100 .098 
1968 .067 .023 .025 Uae .028 .047 6035 
1969 045 .041 ,021 027 .029 Oss -032 
1970 ~150** .007 .016 .020 .013 tULS 2013 
1971 UL .009 .004 .007 .004 .007 O12 
1972 .057 .020 50.17, ,015 .016 .015 .020 
1974 016 .005 -002 .004 .010 .003 - 
1975 .063 .016 .018 .016 017 .014 .019 
1976 .013 .009 .012 -O11 e017 .014 2015 


(b) Monthly and annual averages for entire bay 


gun dul Aug Sep Oct Period 

1967 .054 058.) .091 .080 .064 5075 
1968 .033 -034 .037 .040 .040 036 
1969 .018 sUz3°:- 014 .040 .071 .052 
1970 s000**% 5 =. 000") 4)", 033 i NN) 007 .029 
1971 .039 016 .026 .052 - .029 
1972 - 016 .024 .035 .021 024 
1974 - 002 .004 - 2087 .008 
1975* 014 O17.) .007 soll .010 .011 

# Ved 0S TOLy .021 .029 022 
1976* aQlZ <UL e007 .008 006 .010 

# .016 069 - - .016 .014 
LT “U2o 010 .022 006 .018 .016 

# .030 -OLbe> 025 - - O22 
Note: All figures are in mg/L as P. 


**One questionably high value was omitted. 

Analytical Precision + 0.003 mg/L or 10%, whichever is 
greater. 

Standard deviations of above figures (including station 
and date variation) are generally 50-150% of the means. 
* 4-station weekly surveys 

# 8-station monthly surveys (10 stations in 1977) 
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TABLE 5a 


t-Test of Annual Means of Phosphorus Data, 
Hamilton Harbour 


Station Total P Soluble Reactive P 


1975-76 1975-77 1976-77 1975-76 1975-77 1976-77 
TS ee ee EO 19/92/77 9/0 
é -2.04 NSD NSD NSD NSD 2.36 
252 NSD NSD NSD NSD NSD 2.64 
258 oe 09 NSD NSD NSD Zv62 NSD 
270 -2.64 NSD NSD NSD NSD NSD 


Note: The value of "t" is indicated where Significant at the 5% 
level. 
Negative "t" indicates a Significant decrease in the mean. 
"NSD" = no significant diff erence. 
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TABLE 6 


Outfall Areas with Highest Phosphorus Concentration, 
Hamilton Harbour, 1977 


Typical Values, mg/l] as P 
Location Total P Soluble Reactive P 


eee 


£29 3.1-54d 1.8-1.9 
211 2.0-3.9 1.1-168 
Cie 1.3-1.6 0.6-1.2 
274 0.3-4.7 0.14-0.55 
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Harbour-wide Monthly Means of Iron and Zinc Concentrations, 


1975 
0.31 
0.19 
0.22 
0.34 
0.26 
UeS7 
0.26 
0.28 


Iron 


1976 
0.33 
0.33 
0.19 
0.46 
0.27 
0. 38 
0.16 
0.30 


TABLE 7 


1975-1977 (mg/L) 


Le 
0.18 
0.40 
0.28 
OF35 
0.26 
0.20 
Oe0 
Oey 
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BOD 

CoD 

TOC 
Turbidity 
SS 


Correlation Matrix of 1977 Hamilton Harbour 


TABLE 8 


Organic Pollution Parameters 


BOD COD 
TOGO a0. 70 
1.000 


TOC 

0.156 
0.080 
1.000 
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Turbidity 
QO.1zZ 

-0.004 
0.063 
1.000 


TABLE 9 


Mean Values of Oxygen Demand Indicators, 
Hamilton Harbour, 1976-77 


Location BoD = COD TOC FOC 
Har bour-wi de 

(4 locations) 1976 3.6 Zr 5.5 Vey 
H arbour-wi de 

(4 locations) 1977 Bes 27 6.6 bol 
Outfall Areas 

82 6 26 8 7 
99 4 21 7 6 
103 5 5) 10 i 
108 8 54 1 8 
121 5 35 7 5 
273 30 130 90 ES 
250 (0.5 m depth) 4 32 6 4 
250 (2 m depth) 4 30 6 3 
212 (0.5 m depth) 5 29 7 5 
212 (2-5 m depth) 6 26 6 5 
217 60 130 70 30 
220 8 40 9 8 
225 140 600 1400 50 
274 14 50 20 10 
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ANOXIC OR NEAR-ANOXIC (<Img/ 4) 


FIGURE 3a : DISSOLVED OXYGEN PROFILES (mg/L) HAMILTON HARBOUR 1976, STATION 4. 
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Burlington 
Ship Canal 


0 ] 2 S 4 5 6 rg 8 9 10 i l2 
DISSOLVED OXYGEN (md/t) 


FIGURE 5: DISSOLVED OXYGEN - DEPTH PROFILES, HAMILTON HARBOUR AND 
BURLINGTON SHIP CANAL, JULY 13,1977 EAST TO CENTRAL HARBOUR : 
CANAL, 269, 254, 258 
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FIGURE 24: HARBOUR -WIDE MEANS OF BIOCHEMICAL OXYGEN DEMAND, 
HAMILTON HARBOUR, 1976 —1977 
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HARBOUR-— WIDE MEANS OF CHEMICAL OXYGEN DEMAND, 


FIGURE 2S< 
HAMILTON HARBOUR, 1976-1977 
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FIGURE 26: HARBOUR- WIDE MEANS OF TOTAL ORGANIC CARBON, 
HAMILTON HARBOUR, 1976-1977 
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SUMMAR Y 


The comprehensive sediment surveys conducted in 1975 and 1976 (MOE, 
1977, 1978a) indicated heavy enrichment of pollutants throughout 
most of Hamilton Harbour. Analyses of core samples indicated that 
surficial sediments were enriched by factors of 4 to 140 compared to 
sediments 50 or more centimeters below the sediment-water 

interface. The harbour sediments were found to be most severely 
polluted with respect to zinc, lead and chromium, and to a lesser 
extent, cadmium and mercury. The highest pollutant concentrations 
were found in deep water and adjacent to the outfalls. Enrichment 
next to the outfalls was observed to be a function of the pollutants 
discharged; these samples were classified into separate groupings as 
a result of factor analysis. 


INTRODUCTION 


In May and October 1977, single samples were taken by Ekman dredge 
from 25 locations in Hamilton Harbour. In addition, dredge samples 
were collected at the four major sampling stations (4, 252, 258 and 
270) monthly during the summer. The locations sampled are shown in 
Figure 1. The primary purpose of these surveys was the study of 
bottom fauna (Section J). Samples collected in May were analyzed 
for loss on ignition, free ammonia, total Kjeldahl nitrogen, total 
phosphorus, arsenic, iron, lead, zinc, ether extractables, mercury 
and percent water. In addition to most of the above parameters, 
October samples were analyzed for copper, nickel, cadmium, chromium, 
BOD and COD. Most of these parameters were also determined on the 
monthly summer samples. The results of these analyses are presented 
in this section and compared to the results from the 1976 core 
samples. 
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RESULTS 


Means and standard deviations of all results obtained in May and 
October are given in Table 1. No significant differences were found 
between the two sets of samples except for NH, which was higher 

in October and ether extractables, which was higher in May. 


The increase in anmonia-N could indicate some bacterial 
mineralization of organic N during the summer months. A 
qualitatively similar increase in sediment and interstitial water 
ammonia-N was found during anaerobic incubation of Wisconsin Lake 
sediments (Graetz, Keeney and Aspiras, 1973); unlike Hamilton 
Harbour, ammonia was released to the water column during their 
experiments. As reviewed by Keeney (1973), the interstitial water 
plays a vital role in sediment-water nitrogen relations; microbial 
reactions and redox relationships are all important in determining 
the seasonal nitrogen dependence. Therefore, a proper explanation 
of the observed increase in sediment ammonia concentrations cannot 
be offered without a detailed study of the pore water and exchange 
reactions. The sediment NH 2-N comprises only 1 to 5% of the TKN. 
That the sediment N is mainly organic in nature is not surprising 
for a water body such as Hamilton Harbour, which receives large 
amounts of organic waste inputs (Volume IT). 


Results of the monthly sediment samples (Table la) taken at the four 
major sampling stations indicated no significant changes from one 
month to the next, except for total phosphorus at locations 4 and 
258 in July, when values of 0.7 mg/g were obtained, in contrast to 
values of 4 to 5 mg/g obtained on other dates for these locations. 
This result does not relate to aerator operation or any changes in 
water chemistry (Section C, Water Chemistry), and cannot be 
explained. In addition, low values were obtained for all 
parameters at station 4 in May; results for samples taken during the 
summer at this location (except for July total P mentioned above) 
are similar to October results, implying that the May samples may 
not have been collected at quite the same location. 


HEAVY METALS 


Results for heavy metals in the 1977 dredge samples are similar to, 
or slightly higher than, those for the top portions of core samples 
obtained during 1976 (MOE, 1978a). Except for mercury, these 
results are also almost identical to those obtained for the dredge 
samples taken in May 1975 (MOE, 1977). 


As was done with the 1975 and 1976 data (MOE, 1977, 1978a), the 
concentrations of the 1977 sediment dredge samples were compared to 
values characteristic of "elevated" and "excessive" heavy metal 
concentrations of Great Lakes sediments (Fitchko and Hutchinson, 
1975). These authors considered "elevated" sediment heavy metal 
concentrations to be those of more than the mean + two standard 
deviations of the concentrations from uncontamined Great Lakes river 
mouth areas, and concentrations more than 5 times these values to 
indicate “excessive” values. Figures 2-8 indicate the areas of 
elevated concentrations for Hg, Pb, Zn, Cr, Ni, Cu and Cd for 1977. 
Averages of May and October data were used where available (Hg and 
Pb), otherwise data represent a single survey. Monthly data were 
not included as these represent only four sampling locations. A 
figure for cobalt was not presented as this metal indicated elevated 
concentrations only at the Ottawa Street slip (station 2) and was 
below the "elevated" limit elsewhere. 


Mercury, zinc, lead and chromium show "excessive" concentrations for 
the entire harbour, except for the shallowest areas and the zone 
immediately adjacent to the ship canal. Copper presents "excessive" 
concentrations in most of the deep zone and in the southeastern 
portion of the harbour. Cadmium has "excessive" concentrations in 
most of the deep zone. These metals are "elevated" in most of the 
remainder of the harbour. Nickel presents "excessive" 
concentrations only in the Ottawa Street slip and is "elevated" in 
most other areas except the northern portion of the harbour. 
Comparison of these results to 1976 data is facilitated by 
considering the station number and recorded water depths at the 


USS: 


north shore locations during the October 1977 cruise: 253, 4m; 25/7, 
10m; 260, 12.5 m; and 263, 14.0 m. This shows that only station 253 
can be considered as representative of the "shallow" pairs of north 
shore stations used in 1976 (MOE, 1978a) and indicates that the 
excessive concentrations obtained along the north shore for several 
parameters are not significantly different from 1976 results. When 
water depths at the various stations are considered, no significant 
differences between 1976 and 1977 samples exist, except for 

mercury. In the case of mercury, the 1977 dredge samples (Figure 2) 
present a far larger zone of excessive concentrations than the 
samples for 1976. Excessive concentrations occur along the 
industrial shoreline and in almost the entire deep portion of the 
harbour, except for an area in the eastern part of the harbour which 
was influenced by harbour-lake exchange. Concentrations of mercury 
in the shallow northern and western areas were below the "elevated" 
limit to "elevated" in 1977; this was similar to 1976 considering 
the distribution of sampling locations used in the two years. If 
the dredge sampler took a shallower cut than the 5 to 7 cm depth 
used for core sectioning in 1976, and the heaviest mercury 
contamination is located within the topmost few centimetres, this 
could explain why the 1977 concentrations are higher than those of 
1976. 


The high mercury concentrations observed in the deep water sediments 
are considerably less than those adjacent to the St. Clair River 
chlor-alkali industries, and comparable to downstream areas (Lakes 
St. Clair and Erie) affected by these industries. Values in central 
Lake St. Clair have ranged from 0.7 to 3 ug/g in the period 1973-76 
(MOE, unpublished data), and 1.0 + 0.5 ug/g in western Lake Erie for 
1976 (MOE, 1978b). The highest mercury concentration in Hamilton 
Harbour sediments in May 1977, 2.7 ug/g, occurred at station 262 
(Randles Reef) and appears to indicate that discharge of this metal 
is in that vicinity. This value is similar to the highest recent 
values observed in Lake St. Clair and the western basin of Lake Erie 
(MOE, 1978b; Walters, Wolery and Myser, 1974). The October value at 
the same location was only 0.93 ug/g; this difference presumably 
indicates a large spatial variability rather than an actual 

decline. These values are higher than those for most other Great 
Lakes areas, and indicate a need for further investigation. 
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RESULTS OF FACTOR ANALYSIS 


The results from the May and October 1977 dredge samples were 
processed by Q-mode factor analysis of the similarity coefficients 
obtained from comparison of the ratios of pollutant concentrations 
between the various samples obtained at one time, as was done with 
previous years’ data (MOE, 1977, 1978a). Q-mode factor analysis 
divides samples into different groups with samples of similar 
characteristics being grouped together. The results are shown in 
Figure 9 for May and Figure 10 for October. 


Although minor differences were present between the two series of 
samples, the factor patterns produced were similar to those observed 
in previous years. Both series of samples produced a principal 
factor relating all deep water areas. The stations affected most by 
waste discharges (268, Hamilton STP; 2, Ottawa Street slip; and 262, 
Randles Reef) each produced a unique factor, as did station 264, 
which is influenced by Cootes Paradise. Station 267 presented a 
unique factor in the May series only. Relationships observed in 
1977 along the north shore are less complex than those observed in 
1976, and relate at least partially to water depth at the sampling 
locations. In May, station 252, the shallowest location sampled, 
presented a unique factor, while the remaining north shore points 
were related to the principal factor. In October, one factor 
(factor 2) related stations 252, 253 and 257, while another factor 
(factor 8) related points 260, 263, 15, 270 and 265 in the west end 
of the harbour. Relationships in this part of the harbour were less 
distinct, however. Locations 15 and 270 were also related to the 
principal factor; locations 263 and 265 were related to factor 2 if 
the analysis was stopped at 7 factors. Factors 2 and 8 are similar 
to factors observed in 1976 at shallow and intermediate water depths 
along the north shore and west end of the harbour. Factor 2 
represents the shallowest locations while factor 8 represents 
locations of intermediate depth including the far west end of the 
industrial shoreline as observed in 1976. Factor 8 is sensitive to 
the chemical parameters used in the analysis, and does not appear 
when only the common parameters analyzed for in both months are 
included. Thus its presence appears to be related to the additional 
parameters determined in October. 
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In both months, the effect of lake-harbour exchange is again 
observed at stations 269 and 1030. An important difference, 
however, is the inclusion of station 4 in this factor in May only. 
The latter was completely unexpected, considering the deep water 
nature of this location, and the fact that the sample taken slightly 
to the east (at station 251, closer to the beach strip) was included 
in the principal factor. Pollutant concentrations for the May 
sample at station 4 were much lower than those for samples taken 
from other deep water areas; these were confirmed by repeat analyses 
for several parameters. This result (like the depth profiles in the 
1976 core taken at Station 4) suggests that some transport of sand 
through the Burlington Ship Canal, and southward, is occurring. 
Possibly the May sample was taken from a slightly different 
location; samples taken monthly during the summer at station 4 were 
apparently very similar to the October sample. The effect of Lake 
Ontario is not felt much farther west, as shown by the inclusion of 
Station 22 (not sampled in previous years) with the principal factor. 


It is possible that station 4, and the area immediately to the west 
(station 255) represents a zone of fluctuating transport processes. 
As in 1975 and 1976, no one single factor represented station 255 in 
May; however in October, this station was represented by the 
principal factor. 


As the zones close to the outfalls presented unique factors similar 
to those observed in previous years, they were again used to 
indicate areas which are most affected by waste discharges (Table 
2). The zones with highest pollutant concentrations are Randles 
Reef (Pb, Hg, Zn, Fe, Ether Extractables) and the Ottawa Street Slip 
(Pb, Ni, Cr, Fe, As, Cu, Ether Extractables). The Hamilton sewage 
treatment plant vicinity appeared to be less severely contaminated 
relative to other outfall locations in 1977. A possible source of 
cadmium in the southwestern part of the harbour may be indicated by 
the relatively high value at station 270. This fact, which was also 
observed in 1976, is worth further investigation. 
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Interparameter relationships in 1977 were also similar to previous 
years, as indicated by R-mode factor analysis for May (Table 3) and 
October (Table 4). R-mode factor analysis divides chemical 
parameters into groups with similar characteristics like geochemical 
behaviour or discharge patterns. Transition metals (Fe, Cu, Ni and 
Co) were associated in one factor and post-transition metals (Zn, 
Pb, Cd and Hg) were associated in a second factor. Nutrients (N and 
P) were related in a third factor in May, but were related to 
post-transition metals in October. The fact that both these 
relationships involving N and P have previously been observed does 
indicate year-to-year consistency. Consistency is also seen in the 
various measures of organic matter: COD (October) was related to 
transition metals as observed in 1975 and 1976, but LOI (both sets) 
was related to post-transition metals, as observed in 1975. In both 
months, ether extractables presented a separate factor which was 
somewhat related to Hg, and BOD (October) presented another separate 
factor. It is obvious that different organic components of the 
sediment react differently; a study of the chemical makeup of the 
organic compounds in the sediments would be interesting though 
probably difficult in practice. Somewhat surprisingly, ammonia was 
related most strongly to ether extractables in October. In May, 
ammonia was not included in the factor analysis as it was not 
determined on many of the samples. 


GRAIN SIZE ANALYSIS 


An initial at tempt at grain size analysis, by sieve and hydrometer, 
was performed on samples taken from stations 4, 252, 258 and 270 in 
September. The results are given in Table 5. The sample at station 
252 appeared to consist mainly of fine sand with some coarser 
fractions; unfortunately, the analyses are not consistent, with the 
total (sand and clay) percentage being greater than 100%. The 
results of the deeper locations are similar to each other, as might 
be expected considering other evidence. Considering the fact that 
the finest sediments generally tend to settle to the deepest points, 
the lower proportion of clay and higher proportion of silt at 
station 258 is somewhat surprising. 
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Grain size analysis is important as heavy metal content increases 
significantly as particle size decreases (Forstner, 1976 )t6 lf 
future grain size analysis is performed, chemical analysis should be 
done on various fractions to study the relationship between grain 
size and pollutant content, in particular as it relates to 
distributions throughout various parts of the harbour. 


DISCUSSION 
(a) COMPARISON WITH PREVIOUS YEARS 


In general, the 1977 Hamilton Harbour sediment chemistry results 
were similar to those observed in recent years, taking variations in 
sampling locations and method into account. The October results 
confirm the 1976 finding (MOE, 1978) that the industrial southern 
and southeastern areas produce a greater influence on the deep water 
sediments than do the southwestern part of the harbour and 
moderate-depth northern areas. Although 1977 sediment mercury 
concentrations are generally higher than 1976 values, they are still 
far lower than results obtained for areas in the St. Clair River 
directly affected by the chlor-alkali industries. The mercury 
results for Randles Reef (station 262) in May and October indicate a 
possible high spatial variability as already discussed. 


The body of general sediment chemistry data gathered in the period 
1975 to 1977 represents a good indication of general harbour 
sediment quality. It is suggested that general sediment surveys be 
discontinued for a few years and efforts be directed towards special 
Studies. In particular, the sediment-water interface deserves 
further study as chemical and biological reactions occurring at this 
location play a strong role in controlling aqueous heavy metal and 
nutrient concentrations. Interparameter relationships as depicted 
by R-mode factor analysis of the various sample series indicate that 
relationships between nutrients and heavy metals are more 
complicated than the simple iron-phosphorus-DO system. This topic 
has been discussed more thoroughly (MOE, 1978a, Section D), and 
future projects have been suggested for further study of the 
sediment-water interface. These include partial extraction schemes 


for nutrients and heavy metals in sediment samples, experimental 
measurements of nutrient and heavy metal uptake and release, and 
thorough study of the chemistry of the water from within a few 

centimeters of the interface with the aid of a special water sampler. 


(b) INTERSTITIAL WATER 


In addition to the above-mentioned studies, a study of the 
composition of interstitial water is desirable. This is important 
because the dynamics of nutrient and heavy metal availability to the 
water column is more dependent upon the composition of the 
interstitial water than of the sediments. Interstitial water makes 
up an average of over 50% of Hamilton Harbour surface sediments 
(Table 1), with the highest values being found in the deep-water and 
south shore areas. Interstitial water from reduced sediment cores 
in western Lake Ontario has been observed to contain 2 to 3 orders 
of magnitude more dissolved Fe and 3 to 4 orders of magnitude more 
dissolved Mn compared to lake water (Weiler, 1973). The same author 
also found around 10 mg/L NH. in interstitial water from one of 
these cores. It should be possible to integrate studies of the 
interstitial water with the partial extraction studies. 


(c) DREDGING 


The composition of the interstitial water also plays an important 
role in the environmental impact of dredging, as carried out for 
shipping channel maintenance. It is beyond the scope of the present 
discussion to present a detailed description of the complicated 
factors which control the uptake and/or release of heavy metals and 
nutrients during dredging. Patrick, Gambrell and Khalid (1977) have 
reviewed the factors regulating heavy metal availability during 
dredging. The most important contributor in these reactions is the 
pH - Eh relationship. Dissolved ferrous iron in the interstitial 
waters of the reduced sediment can be rapidly oxidized and 
precipitated as an amorphous gel, which is capable of adsorbing many 
heavy metals. Thus the heavy metal concentrations in the water 
column may remain low provided sufficient dissolved oxygen is 
present. On the other hand, if heavy metal sulfides are present, 
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oxidation can allow the release of heavy metals to the water column, 
although these may be adsorbed later on hydrated ferric oxide 
particles. These mechanisms have been used to explain transient 
pulses of dissolved mercury observed during simulated dredging 
experiments (Lindberg and Harriss, 1977). 


Organic complexation may mobilize many heavy metals which should 
otherwise remain insoluble. This is a complicated process which 
involves both dissolved and insoluble organic materials (Patrick, 
Gambrell and Khalid, 1977), and is bound to be important in the 
highly organic Hamilton Harbour environment. If metals are 
complexed to insoluble high molecular weight organic compounds such 
as humic acids, these may be released under oxidative degradation of 
the organic materials; however, this process is very slow, and the 
released metals would probably be adsorbed by hydrous Fe and Mn 
oxides. More important is the presence of low molecular weight 
soluble organic material in the interstitial water or the water 
column which is capable of complexing significant amounts of 
otherwise insoluble heavy metals. Despite the fact that Patrick, 
Gambrell and Khalid (1977) noted that most studies of water quality 
during dredging operations have shown little or no increase in 
dissolved metal concentrations, the effect of dissolved organic 
matter must not be overlooked. 


Maintenance dredging of the western portion of the Stelco shoreline 
(between stations 262 and 259) and near HHC pier 14 (just west of 
station 262) was performed in 1978. Studies in 1976 (MOE, 1978a) 
have shown that the central portion of the area does not contain any 
unusual contaminants; severe pollution occurs only in the region of 
Randles Reef. Considering the conflicting results concerning the 
effects of redox potential and organic matter reported above, it is 
difficult to predict the effect of such dredging on harbour water 
quality. The most important effect would be the increased DO demand 
which will be exerted by resuspended material. As the sulfide 
content of most harbour sediments is low (MOE, 1978a), the metals 
will not be present as sulfides; therefore, if dissolved oxygen 
levels are depressed, increased heavy metal concentrations may well 
be observed, if only temporarily, in the water column. 
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TABLE 1A 


MEAN RESULTS OF MONTHLY HAMILTON HARBOUR 
SEDIMENT SAMPLES, JULY-SEPTEMBER 1977 


July August September 
Loss on Ignition (%) 8.0 9.7 10.3 
Total Kjeldahl N (mg/g) Ea 357 BE5 
Total P (mg/g) 1.6 Bien) 3.7 
Cu 120. 130. 120. 
Ni 50. 50 Sk. 
Pb 390% 360 330. 
Zn 2200. 3500 3600 
Cd oS 9.2 9.1 
Se 240. 260 240 
Mn 2000. 1900. - 
Fe (mg/g) 80. 120: = 
As 24. 23% 24 
Co ie; ¥2: at 
Hg 0.7 0.7 0.6 


Note: Above results are means of samples taken at stations 4, 252, 258 
and 270. 
All results are in ug/g, unless otherwise noted. 


TABLE 2 


SHALLOW ZONES EXHIBITING HIGHEST POLLUTANT 
CONCENTRATIONS, HAMILTON HARBOUR, 1977 


May October 
Parameter Station Concentration Station Concentration 
Pb 262 1200 (2 860 
262 500 
Cd Not Analyzed 20 Sez 
270 10; 
Hg 262 (aed) 20 1.3 
20 TA 262 0.93 
256 ee 
Ni Not Analyzed 2 190 
Cr Not Analyzed 2 860 
Zn 262 7000 Not Analyzed 
Fe (mg/g) 2 340 Not Analyzed 
256 170 
262 17 
Ether Ext. 262 68 256 ies 
(mg/g) 256 40 267 1632 
268 25 20 14.0 
2 22.6 
As 256 50 2 57 
2 45 
Cu Not Analyzed 2 380 
Note: All results in ug/g unless otherwise stated. 


Variable 


Ether Ext. 


Hg 


R-MODE FACTOR ANALYSIS OF MAY 1977 
HAMILTON HARBOUR SEDIMENT SAMPLES 


Canmunality 


Hohl) 
~ 9758 
9556 
«90835 
59/30 
- 9854 
. 9909 
.9728 
poLOL 


Factor 
1 


S001 
.194 
Rev AS) 
544 
.006 
- 824 
904 
ES 
phy 


TABLE 3 
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Factor 
2 


.054 
002 
.098 
ooh i 
2900 
141 
207 
seo! 
ec 


Factor 
3 


896 


.061 
-296 


Factor 
4 


~405 
248 
.007 
.020 
241 
446 
. 168 
Jo50 
as As) 


TABLE 4 


R-MODE FACTOR ANALYSIS OF OCTOBER 1977 
HAMILTON HARBOUR SEDIMENT SAMPLES 


Variable Communality Factor Factor Factor Factor 
ih Z 3 4 
Cu 9656 902 245 2295 .069 
Ni - 9807 - 960 .180 - LOL .040 
Pb ALB regs: 415 .661 492 3285 
Cd 9187 a5 . 860 aya fs 286 
Cr .9555 896 So 340 HUNKS 
Co 8941 815 463 sO72 103 
As 2/939 ./66 240 .268 er | 
Hg 9233 . 306 Jee: 592 eco 
LOI .8468 494 643 BOOS 242 
BOD 9567 2178 . 388 209 7000 
COD 9306 2/08 2303 462 224 
NH3 .8778 -458 356 e135 .039 
TKN 9542 EEOO 840 344 ~220 
Total P 9126 .c4/ OUT, 081 .065 
Ether Ext. - 9530 400 #333 2153 w 309 
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HAMILTON HARBOUR STUDY 


1977 


VOLUME I 


SECTION E 


PHYTOPLANKTON AND BACTERIA 
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Phytoplankton communites were numerically dominated by the 
Chlorophyta throughout most of the year. This position was briefly 
relinquished to the Cryptophyta during June. The Bacillariophyceae 
were Sub-domi nant throughout the year with peak column 
concentrations appearing in early July. Total cell densities were 
greatest during the spring. A second peak was observed during the 
summer and densities declined into the autumn. Maximum cel] 
densities were invariably found in the upper few metres of the water 
column. These maxima coincided with the euphotic zone depth. 
Community structure was highly dynamic in the harbour, and 
interpretation of these fluctuations was confounded by the unstable 
physical regime and the diverse chemical loadings. 


Determinations of specific bacterial populations showed that the 
heterotrophs were numerically dominant throughout the year. These 
were followed in decreasing densities by the Sulphur oxidizers, 
nitrifiers, sulphur reducers and fecal coliforms. Popul ation 
densities were polymodally distributed over the year and were both 
spatially and temporally heterogeneous. Density fluctuations were 
not simple correlates of substrate concentration, temperature or 
dissolved oxygen concentration. Although neither processes nor 
reaction rates can be ascertained from standing stocks, the 
extremely high densities of the heterotrophs and Sulphur oxidizers 
Suggest that the metabolic processes associated with these bacteria 
were the major sources of oxygen depletion in the harbour. 
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INTRODUCTION 


Hamilton Harbour has been previously characterized as having a 
physically complex environment similar to other Great Lakes coastal 
zones. It has been further described as being hypereutrophic with 
respect to its nutrient concentrations (MOE, 1978). The vertically 
unstable environment combined with the high light attentuation of 
the harbour result in depressed algal productivity and biomass (MOE, 
1978). 


Previous studies have characterized the harbour as having very high 
bacterial standing crops (MOE, 1978). These studies have shown that 
the heterotrophic bacteria were present in the greatest densities 
followed in descending order by the sulphur oxidizers, Nitrosomonas, 
Nitrobacter, fecal coliforms and sulphur reducers. 


The determination of simple relationships between population 
dynamics and environmental variables has not been possible, not only 
due to the highly dynamic environment, but also because of the 
diverse nature of the chemical loadings to the harbour. A more 
detailed analysis based on several years of data will be required to 
gain a greater knowledge of the factors governing algal and 
bacterial community dynamics. 


PHYTOPLANKTON 


Phytoplankton samples were collected weekly at station 258 from 0.2, 
1.0, 3.0, 5.0, 10.0, 15.0, and 20.0 metre depths in order to examine 
vertical community, composition and dynamics (refer to Figure 1). 
Sampling was conducted from early May to mid-October. In addition 
to densities, cell volumes were determined according to the methods 
outlined in previous reports (MOE, 1978; 1977). 
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The phytoplankton communities were comprised of a three phyla: the 


Chlorophyta (greens), the Chrysophyta (diatoms) and the 

Cryptophyta. The phylum Cyanophyta (blue-greens) was also 
represented but only in very low densities at the end of the 

summer. In general, the number of cells per ml were greatest in the 
surf ace waters and declined with depth. Total cell densities were 
as great in the surface 3 m as they were in the remaining 1/m. This 
concentration of cells within the first few meters of the water 
column was expected since the euphctic zone depth (the depth at 
which the irradiance is 1% of the surface level) was previously 
ascertained to be approximately 4 m (MOE, 1978). There was however 
a surprisingly high density of cells at the 20 m depth throughout 
the year (refer to Table 1). This high bottom concentration was 
likely the result of the diffuse and transient thermal structure 
present in 1977 (refer to Figure 3, Section A). These temperature 
profiles show that the classical thermocline conditions of 19% 
change of temperature per metre depth were found very infrequently. 
Therefore exchanges of mass and energy between surface and bottom 
waters may have occurred repeatedly. 


The Chlorophyta comprised most of the biomass throughout the whole 
of the column during the spring. Tetrastrum sp.,.Shroederia sp. and 
Chlamydamonas spp. were the most abundant green algae. The 
cryptophytes were the next most prevalent group with Rhodomonas 
minuta and Cryptomonas ovata being the only members found. Diatoms 
Such as Stephanodiscus astrea, Synedra ulna and Asterionella formosa 
were also present but in low densities. During the early sunmer the 
cryptophytes attained numerical dominance, but by mid-July the green 
algae, such as Oocystis borgei, Pediastrum boryanum and 
Chiamydamonas sp. were once again present in the greatest densities 
(refer to Table 1). 


The diatoms, which were always numerically inferior, were trimodally 
distributed. Peak column densities occurred in July with smaller 
peaks in May and August. The general trends for the population 
densities throughout the water column were similar to those observed 
in the surface metre (see Figures 2, 3 and 4). The total column 
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diatom densities appear to be inversely correlated with the average 
water column contents of silica (Figure 5). All of the four major 
stations have approximately the same trend in their average column 
dissolved reactive silica contents (refer to Figure 5). 


Total cell volumes throughout the water column were greatest during 
the early and late summer. They also appear to be correlated with 
temperature structure since they declined rapidly below thermal 
structures (refer to Figures 6, 7 and 8). 


BACTERIA 


Bacteria samples were collected bi-weekly at the four major sampling 
stations (Figure 1). In addition to the health hazard indicator 
organisms (fecal coliforms), nutrient cycle bacteria associated with 
sulphur, nitrogen and carbon were also examined. Standing stocks of 
the various species were measured at several depths from late April 
until mid-October. 


In 1977, the heterotrophs, which utilize organic carbon, were the 
numerically dominant group. The sulphur oxidizing populations, 
which oxidize reduced forms of sulphur to sulphate, were the next 
most prevalent group. They were present at densities approximately 
one order of magnitude less than the heterotrophs. The Nitrosomonas 
and Nitrobacter populations, which convert ammonia to nitrite and 
nitrite to nitrate, were in total one order of magnitude less than 
the sulphur oxidizers. Their numbers in turn, were approximately an 
order of magnitude greater than the sulphur reducers and fecal 
coliforms, which were present in about equal densities (refer to 
Table 2). Although processes and reaction rates cannot be 
accurately assessed from standing stocks, these results do suggest 
that organic carbon breakdown and sulphur oxidation were the major 
sources of the large and sustained water column oxygen demand. 
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Heterotroph densities were greatest at the surface and declined with 
depth (frgures 95 10, IT) 23/714) 155. [6v 17, fo7 and tg, excep ae 
station 252, which showed vertically homogeneous distributions 
(figures 12 and 13). Station 4 had the largest annual geometric 
mean, followed in descending order by stations 252, 258 and 270. 

The populations were polymodally distributed at all depths over the 


year. 


The sulphur oxidizer populations, in general, declined with depth 
and were also polymodally distributed at all depths over the year 
(figures 9 through 19). The annual geometric means showed that the 
greatest population densities were at station 4, followed by 
stations 252, 258 and 270 (Table 2). 


The sulphur reducer populations, which were at least four orders of 
magnitude smaller than the sulphur oxidizers, showed slightly 
elevated surface concentrations (figures 9 through 19). These 
populations also showed polymodal density distributions. The 
population fluctuations were not directly correlated with those of 
the sulphur oxidizers. 


Nitrosomonas populations showed irregular depth distributions 
(figures 9 through 19). Population densities showed mid-summer and 
late summer or early autumn peaks. The greatest densities were 
found at stations 258 and 270. The Nitrobacter populations were 
greater than or equal to the Nitrosomonas densities. Population 
sizes were in general greater at the bottom. Densities were 
greatest at station 270, followed by stations 258, 252 and 4. 


The annual geometric means showed that the fecal coliform levels 
were slightly elevated over those of 1976 (Table 2). Densities 
frequently exceeded the 100 counts per 100 ml guideline for 
recreational waters. Coliform levels were generally greater at 
station 4, followed by station 252. Stations 258 and 270 had 
approximately equal population densities, which were also less than 
those of station 252 (figures 9 through 19). These populations were 
also polymodally distributed over the year. 
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Population size is unlikely to be uniquely correlated with a Single 
environmental parameter but rather a function of a multivariate 
environment. The major variables are most likely temperature, 
dissolved oxygen concentration, and substrate concentration. 
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Table 2: 


Fecal Coliforms 
Heterotrophs 
Sulphur Oxidizers 
Sulphur Reducers 
Nitrosomonas 
Nitrobacter 


Fecal Coliforms 
Heterotrophs 
Sulphur Oxidizers 
Sulphur Reducers 
Nitrosomonas 
Nitrobacter 


Fecal Coliforms 
Heterotrophs 
Sulphur Oxidizers 
Sulphur Reducers 
Nitrosomonas 
Nitrobacter 


Fecal Coliforms 
Heterotrophs 
Sulphur Oxidizers 
Sulphur Reducers 
Nitrosomonas 
Nitrobacter 


Surf ace 


Station #4 


179 
26,228, 500 
1,947,745 
217 

6,348 
2,467 


48 

14, 803, 900 
315 , 899 
212 

4,051 
3,100 


Station #258 


68 
10, 736, 600 
204 , 853 


Station #270 


49 
8,665,000 
CF5a01 

47 

2,700 

6, 608 


Station #252 
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Middle 


88 
7,783,400 
704,252 


Bacterial Populations of Hamilton Harbour - 1977 
(Annual Geometric Means), Counts/100 mL 


Bottom 


35 
6,709, 700 
41/5728 


96 

14, 984,100 
86 ,341 
274 

8,919 
8,028 


18 
4,110,600 
157,287 


18,473 
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FIGURE 4 : TEMPORAL ABUNDANCE OF CRYPTOPHYTA AT STATION 258 (Im). 
HAMILTON HARBOUR 1977 
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FIGURE 8: TOTAL CELL VOLUMES AND TEMPERATURE STRUCTURE AT STATION 258 
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FIGURE 10 : TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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FIGURE Il : TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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FIGURE 12: TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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FIGURE 13 #TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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FIGURE 14 : TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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FIGURE 15 : TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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FIGURE 16 : TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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FIGURE 17: TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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FIGURE 18 : TEMPORAL ABUNDANCE OF BACTERIA IN HAMILTON HARBOUR 1977. 
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SECTION F 


2-DIMENSTONAL NUMERICAL MODEL 


SUMMAR Y 


The two-dimensional depth-integrated numerical model of currents and 
water quality in Hamilton Harbour has been modified by including the 
real-time variation of currents and water quality at the major 
sources (the Ottawa St. slip and Hamilton STP) and the open boundary 
(exchange with Lake Ontario). This means that short-time variations 
of currents and water quality measured at these locations by 
recording current and water chemistry meters have been read directly 
into the model. Use of these data at intervals of minutes allows 
prediction of the effects on the harbour of large but short-term 
changes of water quality such as those observed at the Burlington 
ship canal. 


Following the procedures used in the recent Toronto Harbour 
modelling study, the model was calibrated to physical parameters 
such as wind stress and bottom roughness by comparing modelled 
output with currents and water quality measured at intervals of 
minutes at a fixed internal point not used as external data input. 
However, it was found that agreement of observed and modelled 
results was far from satisfactory. The location of the best fit was 
several grid squares south of the metered location and moreover 
varied in location with modelled time. In addition, the model is 
incapable of simulating short-term changes in current speed, 
direction or water quality. It was concluded that the major 
deficiency in this modelling effort is the use of a two-dimensional 
model to describe a three-dimensional (depth-stratified) case. 


Despite these shortcomings, computerized maps of water quality do 
indicate the control of water quality by advection (direct transport 
of pollutants by the water masses). This is shown by the shapes of 
the model concentration contours in the southeastern portion of the 
harbour. The contour maps also Show the importance of the 
Burlington STP and Grindstone Creek in affecting water quality near 
their outlets, and the time-varying zone of influence of exchange 
with Lake Ontario. 
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INTRODUCTION 


As part of the original Hamilton Harbour study (MOE, 1974), a 
two-dimensional depth-integrated numerical model was developed which 
is capable of assessing changes in water quality resulting from 
changes in harbour geometry, water withdrawals, and waste 
discharges, for varying meteorological conditions. The model was 
run for a number of conditions representing altered harbour geometry 
and industrial loadings, destratification, increased harbour-1lake 
exchange and varied meteorological conditions. The most effective 
methods of water quality enhancement were found to be artificial 
destratification and a 66% loading reduction of the major waste 
inputs. 


For more precise water quality predictions to be made, however, the 
model must be capable of considering phenomena which occur at very 
short time intervals (order of minutes). The importance of such 
phenomena has been demonstrated for Hamilton Harbour, and more 
particularly, the Burlington Ship Canal, by the recording meter 
results (MOE, 1975; Palmer and Poulton, 1976). This indicates that 
changes in conductivity of several hundred umhos, or in dissolved 
oxygen of over 50% saturation, can occur in a matter of minutes. 
However, the earlier model runs allowed for only the time variation 
of the boundary currents in the ship canal; flow data for the 
industrial and municipal sources and all concentration data were 
entered only as average values. In the present study, real-time 
data obtained from the current and water chemistry meter survey in 
1974 for the Windermere basin, Ottawa Street slip, and the 
Burlington ship canal, have all been used directly in the numerical 
model in an attempt to simulate better the rapid variations in 
conditions known to occur. 


Considering that the harbour is stratified in the summer months and 
the flow regime in the Burlington ship canal is complex with as many 
as three layers of different characteristics having been observed 
(MOE, 1975), a three-dimensional model is required for more complete 
prediction of harbour processes. Following the method used in 
Toronto Harbour (MOE, 1977), it is useful to operate the 
two-dimensional model in Hamilton Harbour and compare the results 
with the available real-time continuously recorded field data. 
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The version of the numerical model used in this study is identical 
to that used previously in Hamilton Harbour with the addition of 
real-time source and boundary inputs. The model has been described 
previously (MOE, 1974, 1977) and this will not be repeated. 


INPUT DATA 


Six pairs of recording current and water quality meters were 
operated in 1974. These installations occupied four distinct 
geographical locations (Figure 1). At the Burlington Ship Canal 
(locations Ml and M2), two Schneider robot monitors and two Plessey 
current meters were operated at depths representative of epilimnetic 
and hypolimnetic water. The offshore location (M6 and M7) comprised 
NERA self-contained water quality monitors as well as current meters 
(Plessey or Geodyne) also at epilimnetic and hypolimnetic depths. 
The other two locations each consisted of a single current meter and 
Plessey submersible water quality meter. The Ottawa Street slip 
installation (M5) was to monitor the combined effluents from this 
source, while the Windermere basin installation (M8) was to monitor 
the discharge from Red Hill Creek, including the Hamilton sewage 
treatment plant effluent. The offshore location (M6 and M7) was 
designed to provide an internal independent calibration point not 
used in the modelling runs. The details of the meter installations 
have been previously described (MOE, 1975). 


The proper functioning of all the instruments would have provided an 
adequate data base for calibration of a three-dimensional model. In 
practice, as already emphasized (MOE, 1975; Palmer and Poulton, 
1977), it is extremely difficult to maintain six sets of recording 
meters all operating simultaneously and in proper calibration. At 
no time were data available from the entire set of recording meters; 
in particular, chemistry meters M5 and M6 and current meter M7 
failed to produce any useful data except for chemistry at M6 after 
July 11,1974. The data fron current meters M5 and M8 were degraded 
by fouling of the rotors with suspended solids and algae; 
approximately 50% of the data record from current meter M8 indicated 
negligible speeds (less than 0.3 cm sul) and the remainder of the 
M8 results indicated slow currents which oscillated in either 
direction. For more details about meter operating periods, see MOE, 
1975, Section D. 
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In addition to these shortcomings caused by Operational problems, 
the data at M8 were unsuitable for direct input into the numerical 
model. As it was impossible to operate recording chemistry and 
current meters for any extended period of time next to the outfall 
of the Hamilton STP, interpolated mean daily flow results measured 
by the Hamilton STP were used with a time variation representative 
of the currents measured offshore of the outfall. Figure 2 shows 
the auto variance density spectra of the north component of M8 
current averaged at half-hourly intervals. The longest significant 
periodicity is 2.6h; this periodicity was superimposed upon the 
daily flows with an amplitude sufficient to give a coefficient of 
variation of about 20%. This represents the variability of the 
Hamilton STP flow as observed in 1972 (MOE, 1974, Section E). A 
constant total dissolved solids (TDS) concentration of 560 mg tl 


TWh 


was used. This represents a conductivity of 860 umho cm 
is an approximate average of measured results in Red Hill Creek, 
which receives the discharge of the Hamilton STP (Section I, 


Volume II). 


With these limitations, the data available for modelling in 
real-time included the Burlington Ship Canal (currents and 
chemistry), the Ottawa Street slip (currents), and the synthesized 
real-time current data set for the Hamilton STP described above. 
Daily flow records for Grindstone Creek, a small source relative to 
the other sources, were interpolated to the model time interval. 
All other flows, and all concentrations except at the Burlington 
Ship Canal were entered as averages. Except for the Burlington STP, 
these constant discharges are small in volume compared to the major 
discharges and lake exchange which were metered. The input flows 
and concentrations used are sunmarized in Table l. 


Two additional corrections in the flow data were required. First, 
at location M5, due to the fouling of the rotor as already 
mentioned, a gradual decrease in meter response occurred. It was 
assumed that 10 days were required for the response to drop 50% and 
an exponential correction factor was applied. Second, due to the 
multi-layered flow regime in the Burlington Ship Canal as described 
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previously (MOE, 1975, Section D), the current meter data set failed 
to describe completely the time variation of currents in the canal. 
In order to prevent apparent accumulation of water in the harbour 
during the model run, it was necessary to subtract 5 cm sol 
algebraically from the boundary current data. 


These corrections provided a consistent input data set to the model 
for the period June 29-July 3, 1974. The real-time input data also 
included hourly wind data as recorded by the Hamilton Harbour police 
at a location near their dock. The input wind record is shown in 
Figure: 3: 


SELECTION OF MODEL PARAMETERS 


As a component of the numerical model study in Toronto Harbour, 
which used real-time data inputs at intervals of minutes, a thorough 
investigation of model parameters was performed (MOE, 1977). This 
involved the checking of the hydrodynamic and mass balance time step 
length, effect of wind stress on a moving surface, and the bottom 
roughness coefficient (Manning's n), together with its interaction 
with water depth. 


The hydrodynamic time step criterion as suggested by Williams and 
Hinwood (1976) and applied to Toronto Harbour was used for the 1974 
Hamilton Harbour model run. The time step length At is expressed as 
a function of the space step Ax, as follows: 


1 < At _/gh mt eee (1) 
AX 


where: AX = Space step length 
g = acceleration due to gravity 
h = maximum depth 


An adequate time step length is obtained when further decreases in 
time step length produce no changes in calculated velocities or 
water levels, e.g. results are independent of time and space scales 
used. 
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Using this formula, and experience gained from modelling in Toronto 
Harbour (MOE, 1977), the required time step length for Hamilton 
Harbour would be of the order of 10 s using the previously employed 
Space step size of 152.4 m. As this would require excessive amounts 
of computation time, a new model grid was set up with a space step 
size of 240 m (see Figure 4). Tests with this grid indicated 
Stability was achieved with At=15 s. Following the methods used in 
Toronto Harbour (MOE, 1977), the time step length for mass balance 
calculations was set to five times this value, or 75 s. 


The other physical parameters affecting model results are the wind 
Stress and bottom roughness (Manning's n) coefficients. In common 
with the previous modelling in Toronto and Hamilton Harbours, the 
wind stress coefficient was set at 0.0032. However, the bottom 
roughness coefficient is not easy to define. Bottom friction is 
ordinarily determined in the field by comparing a measured velocity 
profile to a calculated profile based upon assumed wind stress and 
bottom friction. Although some velocity profiles were measured in 
Hamilton Harbour during 1976, the data were too variable to permit 
any meaningful calculation of Manning's n. 


For the Toronto Harbour modelling study, a literature review of 
typical bottom friction coefficients was conducted, and a grid using 
4 values of Manning's n was established on the basis that shallow 
waters tend to have higher n values than deeper or dredged areas. 
Accordingly, a similar grid was established for Hamilton Harbour 
(Figure 5). Several sets of n values were tested for 12 hours by 
comparing the fit of modelled data to observed data, and it was 
found that the values used in Toronto Harbour (n=0.040 to 0.060) 
were much too high to reproduce the observed currents at location 
M6. Test results suggested use of the values n,=0.020, 

no=0.025, n=0.030 and n,=0.040, and these were employed for a 
model run of longer duration, to be followed by comparison of 
modelled currents with observed currents and further refinement of 


the n values if necessary. 
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MODEL RESULTS 


With the factors affecting model performance as discussed above, a 
numerical model run for 48 hours was performed, using input data for 
the period June 29, 1974 at 12.00 h to July 1, 1974 at 12.00 h. 
After evaluation of the modelled results as discussed below and 
further alteration of Manning's n values, the modelling was extended 
to 96 hours (input data until July 3, 1974). Pertinent input 
parameters used are sunmarized as follows: 


1. Source and sink flow data according to text discussion and 
Table 1. 


2. Average wind field according to Figure 3. 

3. Wind stress coefficient = 0.032. — 

4. Manning's n values as follows: nj=0.020, ng=0.025, 
n3=0.030, and nq=0.040 used on grid shown in Figure 5. 
Hydrodynamic time step = 15 s. 

“Mass balance time step = 75 s. 

Minimum harbour water depth = 1.5 m. 
Mass balance dispersion coefficient = 1.7 més-l. 


On DMN 


The important fact of the volume of input data used in the real time 
modelling cannot be overemphasized. A model period of 48 to 96 h is 
long compared to the longest periodicities of significant magnitude, 
which are 12 h (lunar semi-diurnal tide) and 5 h (Lake Ontario 
uninodal seiche)(MOE, 1975, Section D). Such a model period 
therefore represents a reasonable test for the model. This 96 h of 
model run utilizes the following inputs: 96 wind vector readings, 
576 current meter readings at each of four different locations and 
288 water chemistry readings at each of 2 different locations. The 
validation test utilizes an additional 576 current meter readings 
and 192 water chemistry readings at one location. 

Modelling was also performed with Mannings n= 0.016, ny = 

02023; Ng = 030305 Ng = 0.040, and a= 0-025" Ny = 0.030, 

Leg 05035; Ng = 0.045. For the sake of simplicity, only the 

ny values will be referred to in the following discussion; but 
reference to a particular ny value implies use of the related 


"25 3 and Ng values of that particular group. 
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COMPARISON OF OBSERVED AND MODELLED VELOCITIES 


Observed and modelled velocities were compared by converting the 
calculated U and V components into a resultant Speed and compass 
direction, and plotting these together with the observed currents as 
a function of time for location M6. Different model grid locations 
were used in the vicinity of the current meter in testing the 
modelled and observed results. As will be seen, the location of 
best fit varied somewhat with Manning's n and even with time at 
constant n. The locations used for comparison are shown in Figure 6. 


Although the modelling was first performed for 48 h and later 
extended to 96 h, it is convenient to present data for the entire 
modelled period together. Comparisons of observed and modelled 
velocities are given in Figure 7 for location (14, 8) and n= 
0.016; Figure 8 for location (14, 8) and n, = 0.020; Figure 9 for 
location (14, 7) and ny, = 0.020; and Figure 10 for location (14, 6) 
and ny = 0.025. These locations represent the best fit to 

observed data obtained under the conditions used in each run, after 
summing the differences between observed and modelled velocities for 
several runs and columns near the validation point. 


The initial effort (n, = 0.020; Figure 8a) showed an under- 
estimation of speeds for the first 24 h. During the second 24 h, 
rapid changes in the observed direction were not simulated at all. 
In an attempt to improve this situation, modelling at n= 0.016 

was attempted. The result (Figure 7a) shows a much closer agreement 
to observed currents for the first 32 h. The model approximates the 
observed, rapid direction changes between 37 and 40 h, but does not 
Simulate the remainder of the compass oscillations in the 32 to 48 h 
period. 


In both instances, short period time variations of measured currents 
(order of a few hours or less) were not simulated at all. This was 
true for all modelled periods, and is also emphasized by the model 
autospectra presented later. 
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As the overall result for the first 48 h was somewhat better with 
ie 0.016, this value was used initially for modelling the next 

48 h. The result (Figure 7b) showed a rough simulation of observed 
speeds, but contained two sudden shifts in modelled direction which 
were not observed (the observed current direction was almost 
constant from southeast to south for most of this time period, 
shifting drastically only in the last 4h). Consequently the 
modelled current was at times almost opposite in direction to the 


observed. 


At this stage, the run with her 0.020 was extended to 83 h. The 
comparison at the same grid point (14, 8; Figure 8b) showed 
generally poor agreement, with a shift in modelled current direction 
to opposite the observed direction being found at about 72 h. For 
the same model run, however, the adjacent grid point (14, 7) showed 
a much better fit to the observed data between 50 and 72 h modelled 
time (Figure 9b); however, this grid point was unacceptable during 
the first 48 h (Figure 9a). These two sets of results (Figures 8 
and 9), derived from the same model run (ny = 0.020), showed that 
there are problems in the fit with observed data both with time and 
location. 


In an attempt to improve the output for direction during the last 48 
nh, the model was rerun for this period at ny = 0.025, starting 
from the result at 48 h for ny = 0.020. The best fit between 
observed and modelled results is showed in Figure 10. A good 
agreement between observed and modelled directions is evident for 
most of the modelled time, but speeds are too low for portions of 
the record, with the observed time variations of the speeds not 
being simulated. It is probable that this value of Manning's n 
would not have simulated the observed variations in current 
direction between 32 and 48 h, had the simulations included this 
time interval. 
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In general it is observed that the 2-dimensional model in its 
present form allows a qualitative simulation of currents as measured 
at location M6, but lacks the capability of modelling rapid changes 
in speed or direction. It was necessary to use modelled data from a 
location several grid points away from the actual meter location to 
achieve this limited degree of fit, and the grid location of best 
fit tended to shift from one point to an adjacent point as a 
function of modelled time. Use of the U component from 1 square and 
the V component from the next square (but still adjacent on the 
space - staggered grid) did not give any improvement. The result of 
that test is not detailed here. 


MODEL PREDICTIONS OF HARBOUR-WIDE CIRCULATION 


Predicted currents are shown in Figure 11 after 54 hours of modelled 
time with vs 0.020. This time immediately followed an interval 

of high winds, with the average wind speed for the previous hours 
being 7.5 m sol and the maximum hourly wind speed being 8.5 m 

Sark consistently fran the west. The strongest currents (up to 30 
cm svt) flow with the wind in an eastward direction parallel to 

the north shore, following the shore in shallow water around to the 
vicinity of the Burlington Ship Canal. South of the Canal, in the 
vicinity of the validation point, M6, the current regime is somewhat 
confused. A weaker flow with the wind is also evident next to the 
western part of the industrial shore, the slower speeds being 
related to the larger water depth in this area. In the deep central 
part of the harbour, an average return current is predicted, with 
the strongest return current being adjacent to the north shore 
current. The presence of a coastal jet along the north shore is in 
qualitative agreement with other modelling studies (e.g. Simons, 
1971), who showed a nearshore jet flowing with the wind in Lake 
Ontario, and a central return flow in the deeper waters. 


As the wind regime during the modelled period (Figure 3) was largely 
unidirectional, the overall features of the current pattern are 
similar at all modelled times. Consequently no additional plots of 
harbour-wide currents are presented. 
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MODEL PREDICTIONS OF WATER QUALITY DISTRIBUTION 


As in the previous models of Toronto and Hamilton Harbours (MOE, 
1974, 1977), water quality distribution was approximated by total 
dissolved solids (TDS), as measured by conductivity. This was used 
as it is largely conservative (non-reacting), and real-time data 
were available at 20-minute intervals for two depths in the 
Burlington Ship Canal. Although changes in concentration due to 
chemical and biological reactions are not simulated, the use of TDS 
does represent a reasonable approximation to the spatial 
distribution of water quality as created by advection and dispersion. 


Due to the large spatial gradient of observed conductivity in the 
southeastern part of the harbour close to the major discharges, it 
was decided to start the model with higher TDS concentrations in 
these areas. Results of a conductivity transect done June 14, 1977, 
between Windermere basin and a point west of the mouth of the Ottawa 
Street Slip were converted to total dissolved solids and entered 
into the model grid. This produced an initial total dissolved 
solids contour plot given in Figure 12. Although lake-harbour 
exchange and the other sources (Table 1) do exert a significant 
effect on the harbour TDS, no allowance was made for these sources 
on the initial TDS distribution. 


Computerized plotting representations of TDS contours predicted by 
the model are presented in Figures i3-16 for 24, 48, 60 and 96 h 
modelled time. These figures represent typical water quality 
distributions obtained; in general other time intervals did not 
produce greatly different results and are not reproduced here. 


The effects of the major discharges and their variations in real 
time on harbour water quality in the southeast portion of the 
harbour are evident in the figures. The control of water quality by 
advection is apparent when these figures are compared to the 
modelled currents. For example, for much of the modelled time, the 
water in the southeastern portion of the harbour (near station 20; 
model grid rows 6 to 12) tends to flow towards the southeast close 


to the Burlington beach strip (model column 6) and to the northwest 
closer to the industrial shore (model column 7) (see Figure 11). 
This flow pattern produces the typical shape of the 325 and 340 

mg L-! contours in Figures 14 and 15. As flow is weaker nearest 
the Windermere basin, the position of the 380 mg 7! contour is 
almost constant. 


The effect of lake-harbour exchange and of the other sources can be 
followed from the figures. The only other sources of sufficient 
magnitude to be observed are the Burlington STP and Grindstone 
Creek, which affect small areas in the northeast corner and far west 
end of the harbour, respectively. Although the model results do 
show increased TDS values near Randles Reef, the spatial extent of 
the increase was insufficient to appear on the present maps. Near 
the Burlington STP and Grindstone Creek, TDS concentrations are 
above 340 mg/L after 96 h modelled time. The extent. of influence of 
the Burlington STP is observed to include station 252 (one of the 
four major water chemistry and biology stations). Therefore, it is 
reasonable that effluent discharged by the Burlington STP will 
affect the results of the sampling program at this station and any 
unusual discharge from the STP should be observed within four days 
of discharge. As this source was programmed into the model as a 
constant, no inference about real-time variation of the plume area 
can be drawn from the results. 


The most interesting observed result in the figures is the variation 
of the area affected by lake-harbour exchange as depicted by the 310 
mg ty contour. This contour varies in location from within the 
canal (Figures 13 and 15) to an area encompassing roughly 0.4 km? 
around the harbour end of the canal, including station 269 (figures 
1 and 16). The zone of influence of lake-harbour exchange varies 
rapidly with time between these extremes due to the rapid reversals 
of water flow within the canal, as already well documented (MOE, 
1975; Palmer and Poulton, 1976; Kohli, 1979). Station 269 is 
included within the 310 mg 7! contour roughly 10% of the time. 
Profiling studies at this station indicate the existence of water of 
largely Lake Ontario composition (conductivity below 400 umho 

cm!) within the hypolimnion most of the time but rarely, if ever, 
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in the entire water column. Sediment studies (MOE, 1977b, 1978) 
indicate that this area is affected by lake-harbour exchange, but 
that an area slightly to the west (station 22) is not much affected 


by exchange. 


The results of the present model runs certainly indicate a zone of 
lake-harbour exchange effects within qualitative agreement with 
observed data. However, it should be expected that the three- 
dimensional model will define this zone far better, including an 
estimate of the proportion of time that the hypolimnion of station 
269 is occupied by lake water, as well as the spatial and temporal 


variation of parameters such as dissolved oxygen. 


In Figure 17, a comparison of modelled conductivity (total dissolved 
solids x 1.54) with observed conductivity at location M7 is 
presented. Large variations in the first 36 h of modelled time are 
not simulated at all; during the next 36 h there is closer agreement 
but the oscillations in the two time series do not coincide. 
Although this result was obtained with Manning's n= 0.020, the 
result with Manning's nga 0.016 was similar. A possible reason 

for the inability to simulate the observed time variation of 
conductivity is the observed variation of this parameter with depth, 
as seen in most of the water chemistry profiles at station 4 during 
the summer period. Figure 18 shows a depth profile of conductivity 
as measured July 26, 1977. As this behavior is typical during 
summer stratification, accurate simulation cannot likely be obtained 
using a 2-dimensional model. 


TIME SERIES ANALYSIS OF MODEL PERIODICITIES AND CROSS CORRELATION 
WITH OBSERVED DATA 


A set of model U and V velocity and concentration data at a 10 
minute interval from row 14, column 8 of the model grid (near 
current meter location M6, Figure 1) was analyzed by standard time 
series techniques (Jenkins and Watts, 1968). The results are shown 
in Figures 19 to 21. Unlike the modelled results for Toronto 
Harbour (MOE, 1977), which displayed several highly significant 
periodicities of the order of minutes, more than 99% of the variance 
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was found at near zero frequency. A minor peak appears in the 
V-velocity spectrum at 1.3 h (frequency of 0.75 hv), which is 
one-half the period of input flow data from the Hamilton STP 
discussed earlier. No apparent reason for this periodicity is 
evident. 


Cross-correlation of the modelled data at several grid locations 
with observed currents and chemistry at locations M6 and M7 was 
attempted. The only significant coherence in some cases was the V 
component at 2.6 h (the input periodicity of Hamilton STP flow). 
Few or no other significant coherences were obtained. These 
results, as well as the model-meter comparison presented earlier 
(Figures 7-10 and 17) emphasize the present inability of the model 
to simulate events occurring over short periods despite the use of 
real-time data for the canal and major discharges over intervals of 
minutes. 


DISCUSSION 


The comparisons between observed and modelled results as presented 
indicate that the present model is able to simulate the harbour 
currents and water quality on a qualitative basis, but that very 
little short-term variations were obtained. Despite the fact that 
currents were entered at 10 minute intervals for the major sources 
and open boundary, and chemistry values were entered at 20 minute 
intervals on the open boundary, the only significant 

cross-correl ation between the model and independent current meter 
data was for the 2.6 h periodicity which was Superimposed upon the 
input flows from the Hamilton STP. 


To obtain even a qualitative agreement with observed currents 
(Figures 7 - 10), it was necessary to select model grid points 
(Figure 6) which were somewhat closer to the industrial shoreline in 
the vicinity of the Ottawa Street slip. In addition, the location 
of the grid point with the best fit changed somewhat with time over 
the 96 h modelled period. Figure 11 shows that there is 
considerable local variation of currents in this part of the 


harbour, with the strongest currents near the eastern shore (beach 
strip) in a northwest to southeast direction, opposite to observed 
currents at meter location M6. These differences may be related to 
the actual spatial variations of water depth in this area, or more 
probably to the existence of thermal and chemical stratification. 
Three-dimensional model studies in the Great Lakes (for example, 
Simons, 1975) tend to show return currents against the wind in deep 
water areas, with frequent occurrences of current shear patterns 
(surface water transport with the wind and bottom water transport 
against the wind). The variations in fit between modelled and 
observed values at different times for separate grid points may be 
related to the use of a two-dimensional model in a stratified 
situation, and possibly to the limits on spatial resolution imposed 
by the present model grid. 


The validity of the present modelling effort is also highly 
dependent upon the values of input physical parameters such as 
surface and bottom stress, as already mentioned. Selection of these 
parameters is a difficult and expensive task as it involves repeated 
runs of the model with small variations in the parameters, followed 
by comparison of observed and modelled data for various grid 

points. In his independent work on numerical modelling of Hamilton 
Harbour using the same input data as used here, Rasmussen (1978) 
used two statistical methods to evaluate model fit. These were: 

(a) least-squares comparison between modelled and observed 
currents; (b) “minimax criterion", or minimizing the maximum model 
error (absolute value of difference between modelled and observed 
current). Jsing these criteria for 12 hours of modelled time, he 
obtained the best fit to observed currents with a wind-stress 
coefficient of 0.0032 and Manning's n,=0.0126 to 0.0130 (as in the 
present study, a grid of higher n values was used in rougher, 
Shallower areas). His wind-stress coefficient is the same as used 
here but the Manning's n is lower (corresponding to higher currents 
and more rapid changes in direction). 


During storm periods as presently modelled, the effect of wind 
stress is important during periods of strong wind whereas the bottom 
stress controls the rate of current decay after the event (Simons, 
1973). Additionally, bottom stress can be specified in terms of 
several different formulations (Cheng, Powell and Dillon, 1976); 
indeed, the formulation used has a bearing on the grid dispersion or 
smoothness parameter (local variation of the water level). Although 
more complicated in computational effort, the non-linear formulation 
used in the present Hamilton Harbour model emphasizes the depth 
dependence on bottom roughness with strong frictional effects for 
high currents and shallow depths. 


Spatial and temporal variations of surface stress, which were not 
invoked in the present study, could also be included. Donelan 
(1978) suggests that wind drag is highly dependent upon wave fields 
and thus upon fetch. Therefore, the wind stress coefficient should 
be a function of location and wind direction, increasing as the 
upwind distance from a grid point to the upwind shore increases. 


A significant physical parameter which is involved in the water 
quality model is dispersion. This includes the effect of sub-grid 
scale turbulence on both currents and water quality. Leendertse and 
Liu (1977) have suggested that vertical and horizontal turbulent 
exchanges play an important role in determining sub-grid scale 
dispersion, and they introduce modifications to their 
three-dimensional model to account for these effects. The vertical 
exchange is related to turbulent effects caused by flow moving over 
the bottom, air moving over water, and velocity gradients in the 
vertical direction. Horizontal exchanges are caused by large 
two-dimensional horizontal eddies. Non-linear horizontal eddy 
viscosity coefficients are introduced to accomplish the dissipation 
of vorticity which cascades from larger scales to smaller scales as 
a result of the two-dimensional horizontal turbulence. This 
coefficient is a function of the local velocity deformation 
calculated in the finite difference grid. It is possible that 


neglect of these turbulent exchanges in the present model could 
account at least partly for the failure of the model to simulate 
short-term exchanges observed at the validation location. 


Despite these limitations, the dispersion coefficient presently used 
does exert a smoothing effect on concentration contours and helps to 
counteract numerical instability effects (concentration waves) 
produced by the finite difference approximations near sources of 
high concentrations. This effect of dispersion was verified by 
performing a model run with the dispersion coefficient reduced by a 
factor of 10. The resulting output contained significant 
concentration waves in the area closest to the largest waste 


discharges. 


The basic shortcoming in the present modelling effort, and probably 
the fundamental reason for the limited success in simulating 
observed currents and water quality, is the use of a two-dimensional 
model to describe a three-dimensional problem. As already pointed 
out, other studies have shown a considerable effect of water depth 
and thermal stratification on water movements, and hence, water 
quality. Therefore, further effort be directed to the three- 
dimensional model, thus allowing for appropriate simulation of the 
Stratified case. In addition, dissolved oxygen concentrations may 
be modelled, using input data from the monitors which provide 
dissolved oxygen values at 20 minute intervals. This would include 
the effect of sediment and water column oxygen demand, as measured 
in the last few years (Section D; MOE, 1978, Section D). Given 
proper input data, additional parameters could be included. 
Furthermore, the model could be extended to multiple grid sizes, 
thus allowing a finer grid to be used in the southeastern part of 
the harbour (for example, Wickramartne et al, 1976). This would 
allow simulation of water quality gradients in better detail at 
minimum increased computer costs and provide for testing the effect 
of proposed construction projects such as landfilling along the 
beach strip, and alternatives for management of Windermere Basin. 


Although it has been shown that the 2-dimensional model fails to 
reproduce the detailed time-dependence of observed currents, the 
model is still useful for approximate predictions of water quality 
under non-stratified conditions. Concentration contours obtained 
using different meteorological conditions can be used to define the 
zones of greatest water quality impairment, as well as the extent of 
reduction upon reducing waste loads, and thus be useful in defining 
mixing zones or limited-use zones about waste outfalls as proposed 
by the Ontario Ministry of the Environment and the IJC. 
Furthermore, as recently done in Toronto Harbour (MOE, 1980), the 
model can be used to provide an estimate of the areas affected by 
stormwater runoffs in short time scales (minutes to hours) and the 
effect of diversion of such runoffs. 
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TABLE 1 


SOURCE AND SINK DATA USED FOR NUMERICAL MODEL 
HAMILTON HARBOUR, JUNE-JULY 1974 


Flow Conductivity 
Map I.D. Outfall or Intake (m3 s-!) (umho cm-1) 

A Hamilton WPCP 3.68* 860 

Slocts: 

2.98 

3.30 

3.26 
B Ottawa St. Slip** metered 592 
G Stelco Intake !2 -14.9 ~ 
D Stelco S. Open Cut Sewer Cat 677 
E Stelco N. Trunk Sewer 3.04 477 
r Stelco Intake !1 -2.40 - 
G Stelco !3 Open Hearth Sewer 2.40 569 
H Stelco 148" Plate Mill 0. 94 whet’ 
I Dof asco raw water intake -8.9 - 
J Dof asco (coke oven by-products, 2.74 628 

boiler house discharge, silicon 
steel) 

K Skyway WPCP 6.83 738 
L Dundas WPCP 0.066 738 
M Grindstone Creek 0.64* 625 

0.57 

0.49 

0.40 

0.34 

Note: 


Negative flows are intakes (concentrations not required). 

* Daily flows for 5-day period June 29-July 3, 1974 were 
interpolated. A variation of 20% was added to source A with a 
et of 2.6 hr (longest observed periodicity at location 


** Includes: Stelco (£. side lagoon, HC] regenerator and cold 
mill, hot strip mill, heavy gauge shear line), Dofasco (lagoon, 
eae sewer) and CIL (total plant less service water 
intake). 
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REFINED SPECTRAL DENSITY ESTIMATE 
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FIGURE 2 - VARIANCE DENSITY SPECTRA OF NORTH COMPONENT OF 
CURRENTS, LOCATION M8, HAMILTON HARBOUR 
JUNE-JULY, 1974. 
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MODELLED TIME, h. 
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FIGURE 7a - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES FOR LOCATION (14,8) AND 
MANNING’S N = .016,.023,.030,.040, 
STARTING AT 1200 h. ON 740629. 
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FIGURE 7b - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES FOR LOCATION (I4,8) AND 
MANNING’S N = .016,.023, .030,.040, 
STARTING AT 1200 h. ON 740701. 
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FIGURE 8a - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES FOR LOCATION (/4,8) AND 
MANNING’S N = .020,.025, .030, .040, 
STARTING AT 1200 h. ON 740629. 
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FIGURE 8b - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES FOR LOCATION (14,8) AND 
MANNING’S N = .020, .025, .030, .040, 
STARTING AT 1200 h. ON 740701. 
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FIGURE 9a - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES FOR LOCATION (14,7) AND 
MANNING’S N = .020, .025, .030, .040, 
STARTING AT 1200 h. ON 740629. 
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FIGURE 9b - COMPARISON OF OBSERVED AND MODELLED 


VELOCITIES FOR LOCATION (14,7) AND 
MANNING’S N = .020, .025, .030, .040, 


STARTING AT 1200 h. ON 74070). 
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FIGURE 10 - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES FOR LOCATION (14,6) AND 
MANNING S N = .025, .030, .035, .045, 


STARTING AT 1200 h. ON 74070I. 
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FIGURE |i - NUMERICAL MODEL PREDICTION OF CURRENTS AT 54 h. MODELLED TIME. 
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FIGURE 17 - COMPARISON OF OBSERVED (LOCATION M7) 
AND MODELLED CONDUCTIVITY FOR LOCATION 


(14,8) AND MANNING’S N 


STARTING AT 1200 h. ON 740629. 
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FIGURE 18 - DEPTH VARIATION OF CONDUCTIVITY, STATION 4, 
JULY 26,1977. 
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FIGURE 19 - AUTOSPECTRA OF MODELLED U VELOCITY AT LOCATION 
(14,8), HAMILTON HARBOUR. 
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FIGURE 20 - AUTOSPECTRA OF MODELLED V VELOCITY AT LOCATION 
(14,8), HAMILTON HARBOUR. 


F-43 


SMOOTHED SPECTRAL DENSITY ESTIMATE 


0.1 


0.01 


0.001 


0.0001 


0.0000! 


0 0.3 0.6 0.9 l.2 1.5 1.8 2.1 2.4 


FREQUENCY, h_! 


FIGURE 21 - AUTOSPECTRA OF MODELLED CONCENTRATION AT LOCATION 
(14,8), HAMILTON HARBOUR. 
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SUMMAR Y 


The three-dimensional numerical model of hydrodynamics, temperature 
and water quality as developed by Leendertse et al (1973) has been 
modified for use on lake systems and initial calibration runs of the 
hydrodynamic portion have been performed for Hamilton Harbour. The 
theory of the model is described in detail, with particular 
reference to the physical parameters for which calibration values 
must be specified. These parameters include coefficients for 
surface, interfacial and bottom stress, horizontal dispersion of 
momentum, and horizontal and vertical diffusion of heat. 


The model was operated using real-time input data for winds and 
currents in the Burlington Ship Canal obtained in June 1974 and June 
1976. Model performance was assessed by comparing modelled currents 
to those measured at an internal harbour location not used in model 
input, and attempts were made at improving the fit to observed 
currents by varying the values of the input physical parameters. 

The best results were achieved with a wind stress coefficient of 
0.020, a Chezy coefficient of bottom roughness = 700 cm’/271 | 
coefficient of interfacial friction = 1000 cm’, and dispersion 
coefficient for horizontal momentum = 10° cmés~!, However, 

modelled currents failed to follow the observed direction changes at 
the internal validation point; in particular, modelled middle layer 
currents were frequently opposite in direction to observed currents 
while bottom layer currents showed better agreement, despite the 
depth of the current meters used for validation, which corresponded 
to the middle layer. This suggests that because of harbour vertical 
instability, water flowing past the fixed current meter locations 
may be representative of different model layers at different times. 


Based on the results obtained so far, it is suggested that future 
modelling works should involve a larger number of model vertical 
layers, as well as vertical strings of current meters for 
calibration processes. Such a study is ambitious and expensive, 
however, and its merits should be carefully reviewed before 
undertaking. The simpler version of the model as outlined in the 
report may still be useful for a description of the effect of 
changing wind and lake conditions on average harbour currents, as 
illustrated by a series of computerized plots of modelled currents. 
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INTRODUCTION 


The development of a 2-dimensional model was included as a portion 
of the original Hamilton Harbour Study objectives (MOE, 1974). This 
model involved predictions of depth-integrated currents and total 
dissolved solids (TDS). The model has been extended (Section F) by 
adding input of currents and TDS values at the major waste sources 
(Hamilton STP and Ottawa Street slip) and at the open boundary 
(Burlington Ship Canal) at short-time scale intervals of 10 to 20 
minutes. These short intervals were required for adequate 
description of the rapidly varying conditions which were found to 
exist, particularly in the Burlington Ship Canal. 


Canparison of observed and modelled currents at an internal 
validation point showed inadequate agreement, with predicted 
currents opposite in direction to those observed for a portion of 
the modelled interval. It was concluded that errors were primarily 
caused by the use of a 2-dimensional depth-integrated model to 
describe a thermally stratified (June-July) condition, and that a 
3-dimensional model would provide an improved simulation. The 
present report summarizes initial efforts at applying a 
3-dimensional model to the problem of Hamilton Harbour circulation 
under conditions of summer stratification. 


THEORY 


Larger scale water motions in lakes are primarily caused by wind and 
modified by interfacial and bottom friction stresses. Motion due to 
inflows and outflows of water is normally observable only 
immediately adjacent to these flows. As natural lake basins are 
geometrically very complicated and wind motions highly variable both 
in space and time, an exact description of 3-dimensional time- 
dependant motions would require an infinite number of parameters. 


FORMULATION 


These fluid motions may be described by the classical hydrodynamic 
equations of momentum, continuity and mass transport (Leendertse, 
Alexander and Liu, 1973; Lick, 1976; Cheng, Powell and Dillon, 1976). 
These are conceptually similar to the 2-dimensional equations 
previously used in Hamilton Harbour (MOE, 1974), and are as follows: 


(a) Momentum equations 
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The symbols used will be defined below; the significance of the 
terms in equations (1) and (2) is as follows: 


Variation of velocity component with time. 


Convective transports of momentum in the x-, y-, and z- 
In larger water bodies like Lake Ontario, 


these terms are often neglected with resulting 


Simplification of computations (Csanady, 1978). 


This 


approximation is reasonable when the Rossby number is less 


For Hamilton Harbour, 


a Rossby number of 0.1 corresponds to a current speed of 
therefore these terms should be included 
for computational accuracy, though James and Eid (1978) 


[A] 
[B] 
directions. 
than about 0.1 (Cheng et al, 1976). 
about 8 cm sv}, 
chose to neglect them. 
ea 


Coriolis force due to earth's rotation. 


[D] Manentum change due to pressure variations in x- and y- 
directions. These pressure variations are caused by water 
level differences induced by wind set-up. 


[EJ Horizontal dispersive transport of momentum. These terms 
approximate the horizontal sub-grid scale turbulent 
stresses (Reynolds stresses). 


[F] Vertical dispersive transport of momentum. This term 
includes the effect of frictional stresses at the water 
surface and bottom, and between layers in the vertical. 


The vertical equation of momentum has been replaced by the 
hydrostatic approximation equation (3), as vertical accelerations 
may be neglected compared to horizontal accelerations. 


(b) Continuity equation (conservation of water volume). 


Fahtiay aes (4) 
The terms of this equation express the changes in water volume due 
to velocity changes in each direction. This assumes that water is 
incompressible, and that density variations are allowed to affect 
only the vertical acceleration (Boussinesq approximation). 


(c) Mass and heat transport equations 


agit aU ad Pp, $5 - 2 fe 8] = 0 (5) 
LB] Lc] 
[A] variation of mass or heat with time. 
[B] advective transport of mass or heat in the x-, y-, and z- 
directions. 
[Cc] diffusion of mass or heat in the x-, y-, and z- directions. 
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Equation (5) may be used for non-conservative substances (dissolved 
oxygen, BOD, bacteria, etc.) upon insertion of the approporiate 


decay terms. 
(d) Equation of state 
0 = pee et) (7) 


This equation expresses density changes as a function of temperature. 
In fresh water environments, small density changes due to dissolved 


solids may be neglected. 


The symbols in equations (1) - (7) are as follows: 


XV ez Cartesian coordinates positive 
eastward, northward, and upward, 
respectively 


u, V, W Respective components of velocity 

t Time 

a Coriolis parameter 

p Pressure 

S Dissolved solids concentration 

a Temperature 

re) Density 

re) Reference density, a constant 

rod Departure from p depending on 
temperature 

g acceleration due to gravity 

K Vertical mass diffusion coefficient 

KR? vertical thermal diffusion 


coefficient 


AH, AV Horizontal and vertical eddy 
viscosity coefficients 


Dix DY. Horizontal diffusion coefficients 


These equations as written above introduce no internal sources of 
momentum, fluid or solids; in the case of a harbour receiving waste 
discharges, appropriate source terms must be added. To complete the 
description of the model, appropriate initial and boundary 
conditions are added. The initial condition states that all 
velocities and water levels are zero; alternately, the final 
velocity pattern from a previous run (such as a spin-up) may be used 
as an initial condition. Boundary conditions require the 
specification of current, temperatures and salinities at all open 
boundaries (such as the Burlington Ship Canal); the input used for 
open boundaries is discussed in detail later. At closed (land) 
boundaries, the boundary condition is that all velocity components 
normal to the land boundary are zero. 


VERTICAL DISCRETIZATION 


In thermally stratified lakes and harbours, considerable changes can 
occur in the properties described by equations (1) - (6) within a 
short vertical distance. In order to achieve a reasonable 
approximation to this vertical variation within a finite amount of 
computer time and memory, the water body is divided into a number of 
discrete layers and the above equations are vertically integrated 
over the thickness of each layer. This vertical integration 
procedure is described by Leendertse, Alexander and Liu (1973); it 
suffices here to present the vertically integrated equations along 
with the assumptions used in their derivation together with the 
appropriate initial and boundary conditions. 


3 (hu) 3 (huu) 3 (huv) at a: i" Tear, 
pie ener eee igre 7 OU) eat (yD) apis yo pik ox 
PEP Wa a Stelse file the wert Ou: eas io = si = 0 (8) 
Dp wx Ay ox 0 Oy Ay Oy 0 ae 0 ets 
3 (hv) 3 (hvu) 0 (hvv) Z hi SEH ee hae) 
ot ui ox z oy (WV) ey (WV) as, Ce oy 
yz yz 
BHO ov = i ce ov Lat fs it = 0 
ery | AG 0 oy [ AG 5 E | k- |p k¥s (9) 
b 
mh, Ne tel d(hu) , d(hv) eset (10a) 
k-3 Ox dy 
=k 
Ig 
ge RSE TE oChu) , d¢hv)| 4 5 k=1 (10b) 
dt Q=1 OX oy 
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3 (hs) a 3 (hus) ‘. 3 (hvs) 
ot Ox oy 


0 (hT) 0 (huT ) + 


3 (hvT) 
ot ox ) 


In these equations, the quantities u, v, w, Ss, t and p are vertical 
averages over the layer k under consideration. A subscript k - ¥ or 
k + 45 represents a valve calculated at the interface between layer k 
and layer k -1 or k +1, respectively. See figure 1 for 
configuration of layers. 


In addition, 7 ** and7 9% = components of surface, interfacial or 
bottom stress. 


= any layer 

= layer being computed 

bottom layer 

= water level displacement from equilibrium level 
= thickness of layer 


Th (Saal Ope roe CS 
il 


= rate of water level increase due to addition of water at 
waste source point (= flow rate divided by area of grid 
Square into which discharge flows). 


The last term in equations (10b), (11) and (12) represents the 
addition of water, pollutant and heat, respectively, to the system 
from an outfall or its removal at an intake. An intake is 
represented by a negative F value. At present, all outfalls and 
intakes are assumed to be at the surface level (k = 1). 
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SPECIFICATION OF PRESSURE AND STRESS TERMS 


Equations (8) - (12) contain Stress terms which must be evaluated 
for every model layer as part of the boundary conditions for mode] 
solution. As the value of these terms is determined by physical 
constants whose value must be Specified for model calibration, their 
evaluation will be Presented in some detail. 


The pressure gradients in equations (8) and (9) are defined for the 
top layer (k = 1) as: 


is ee Ries 9P 4 
ax BR 9 FT Sashage 
OP} eka 1 OP) 
acy = —_— +k hy — 
Dy SP] y 2 8n} dy 


For the other layers (k>1) the pressure gradient is given by 


oP " oF et pre reese 
ax x ecket ay 
oP. ? ae — Oy ok 
dy ay 2 OES 


In the equations the Subscript k - } represents an average of the 
value between layers k - 1 and kK. The above equations are derived 
on the asumption that the Spatial variation of atmospheric pressure 
across the water body is negligible. 


XZ Wee yz Ws. he 
The stress terms Reiee Ek T ass and Ty In equations (8) and 


(9), which arise from the vertical dispersion terms [F] in equations 
(1) and (2), require careful and complete consideration as the mode] 
solution is highly sensitive to their values, which are determined 

in part by the choice of physical parameters contained therein. For 


the top layer (k = 1), the terms Te a and Ty fy) represent the effect 
2 2 

of wind stress at the air-water interface. At the bottom 

(k = b), the terms 74 and Ti represent the frictional effect of 


water moving over the sediment-water interface. At intermediate 
layers, these terms represent turbulent effects due to interlayer 
density and velocity differences. Each term will be considered in 


detail. 
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(13a) 


(13b) 


(14a) 


(14b) 


The surface stress terms are represented by quadratic functions of 


wind velocity as follows: 


tX> = Op We sin ¥ (15a) 
YS = Oy We cos ¥ (15b) 


where s (Superscript) = surface 


© = wind-stress coefficient 

ae = density of air 

W, = wind speed 

My = angle between wind direction and the model y - axis 


The bottom stress terms can also be represented by quadratic 
expressions: 


ee -. U u2 + v2) (16a) 
2 G 

aN e V (ue ty 2y%2 (16b) 
C 


where b (superscript) = bottom layer 
g = acceleration due to gravity 
C = Chezy's coefficient of bottom roughness 
(dimension L2 Tl ) 


Chezy's coefficient in equation (14) can in turn be expressed as a 
function of total water depth H by the Manning equation: 


Cie 1.9 r 1/6 


where n = Manning's constant (dimension 


(17) 
| -1/3 4) 


The corresponding stress terms in the constituent and heat balance 
equations (11) and (12) must also be considered. For a conservative 
Substance in equation (11), these terms are zero. For heat, the 
surface term (x/-2) for k = 1 represents heat flux at the 
Surface. An exact specification of this term requires that it 


consider heat transfer by mechanisms such as radiation, convection 
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and evaporation, which are all time varying. A simplified 
expression for heat flux q at the surface (Lick, 1976) is dependent 
only upon the air temperature ies and the model surface layer 
temperature qT): Integrated over the model surface layer, this 
expression becomes: 


Soe 
f =| = q°= Hhj e = 7,| (18) 
Ss 5 
where q = heat flux at the surface 
H = surface heat transfer coefficient 


h, = thickness of model surface layer 


At the bottom, a similar heat transfer coefficient may be used if 
the sediment temperature is known or assumed. Alternately, the heat 
flux can be considered as zero for the length of model operation 
used here. 


Several approaches were considered by Leendertse and Liu (1975) for 
the interfacial stress terms, which can be written simply for 
dissolved solids, temperature, x-, and y- momentum respectively as 
follows (see equations 8, 9, 11 and 12): 


ce] OT 
K Ss K! a» 7* and 


jel 
The stress terms 7%? and 74% arise from vertical integration of 

the vertical eddy viscosity terms [F] in equations (1) and (2). 
These terms are proportional to the velocity gradient in the 
vertical direction (see equations (1) and (2)); the proportionality 
constant is the eddy viscosity coefficient Ay. The sub-grid scale 
vertical turbulent fluctuations are then assumed to appear as 
increased viscous flow stresses (Reynolds stresses). Eddy 
conductivity coefficients are similarly introduced so that mass and 
heat fluxes are proportional to concentration and temperature 
gradients, respectively. In turbulent flow, these coefficients also 
depend upon the flow itself as this is the process generating the 
turbulence. 


In constant-density systems, the eddy viscosity coefficients may be 
represented by quadratic terms similar to those used for surface and 
bottom stress (Leendertse, Alexander and Liu, 1973): 


2 2) 4s 
AY, au if XZ, 9 AU AU AV 
ey ; y | ne, fee Ps iz (19a) 
2 2)% 
Seen Leeme yz 2 AV. J} jAU | AN (19b) 
Nv b7aa 5 \ het eA eee Hea a fz 
a) Rs 
dS 2 As y;AU AV 19c 
f 3 regen AG fra i rg ele 
ve oe ers a peace be 
OD INR i AZ (19d) 


Corresponding to the discrete depth formulation of Figure A, the 
partial differentials 0/dz have been replaced by the difference 
A/Az between the centers of the layers. The parameters v , K,,- 
and ki, are coefficients of interfacial friction. These 
coefficients (dimension L2) are the squares of the Prandtl's 

mixing lengths for momentum, salinity and temperature (Odd & Roger, 
1978, equation 2). Increased turbulence (increased Ay) manifests 


itself in an increased value of these coefficients. 


When vertical density gradients exist, the intensity of the 
turbulence is also affected by the buoyancy forces induced by the 
density gradient. Increased stratification suppresses turbulence; 
the onset of this Suppression may be determined from the value of 
the Richardson number, as follows: 


puateen€ 0p 02) 
Rate 9 (3U/dz)2 (20) 


where U is the velocity component under considerations for u and v 
velocity calculations or the resultant velocity for temperature and 
concentration calculations. Turbulence is maintained at Ri values 
below about 0.25 and suppressed at larger Ri values. 
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Leendertse and Liu (1975) investigated several formulas relating 
interfacial stress to the Richardson number, and selected an 
exponential form proposed by Mamayev (1958; cited by Leendertse and 
Liu, 1975): 
—mRi (21a) 


icCi=cac=are (21b) 
oO 


Values of 1.5 and 3.0 were suggested for the empirical coefficients 
m and r, respectively; vo is the vertical momentum transfer 
coefficient under neutral stability (Ri = 0). 


These two contributions to vertical eddy viscosity are combined in 
the present model to yield expressions for the 7 and Kk terms in . 


equations (8), (9), (11) and (12) as follows: 


1 = ha fe beh ap kek >| aoe a) | (i (22a) 
0 Clea fe exp(-m.Ri) + ( ~ = | rg (22b) 
; ce ¥ i. exp(-r.Ri) + [| + fe] | ! fy (32a) 
ee x I exp(-r.Ri) + ale) x cil | EB oan 


The differences are taken between layers k and k+l for subscript 
k+4s; for subscript k-s they are taken between layers k- and k . 
The coefficients Yo, Wy, Ky and Kj, occurring in the above 
equations are defined above, and their values must be specified as 
part of the physical constants affecting model calibration. 


The effect of stratification (Ri > 0) of the interfacial stress is 
thus seen as equivalent to a decrease in the value of the 
coefficient of interfacial friction. This is physically equivalent 
to a higher surface water velocity, and lower subsurface return 
velocities (but return velocities occurring over a larger portion of 


the water column). 


In addition to these vertical stress terms, there also exists the 
horizontal eddy viscosity terms [£] in equations (1) and (2). The 
eddy viscosity Au js another physical parameter to be specified 
during model calibration; an increased A, results in considerable 
smoothing of the horizontal velocities (Simons, 1971). The effect 
of this parameter, however, is generally less than that of Ay 
(Simons, 1974); it is often omitted from hydrodynamic models, 
including the 2-dimensional Hamilton Harbour model (MOE, 1974). 


FINITE-DIFFERENCE FORMULATIONS 


As with the 2-dimensional model (MOE, 1974), the coupled system of 
equations (8) - (12) is solved numerically by a finite-difference 
method based upon the space-staggered grid lattice shown in Figure 
1. In this grid the indices i, j and k represent single steps along 
the grid lattice in the directures x, y and z respectively. 
Temperatures and concentrations are calculated at integer values 

i, J and ks; u velocity components are calculated at half-integer 
values of i and integer values of j and k; v velocity components are 
calculated at half-integer values of k and integer values of i and j; 
and w velocity components are calculated at half-integer values of k 
and integer values of i and j. The resulting finite-difference 
expressions for equations (8) - (12) are described by Leendertse, 
Alexander and Liu (1973, 1975), and will not be repeated here. 


Several different integration methods are available for solution of 
the finite-difference equations. The integration schemes must meet 
three basic requirements (Cheng et al, 1976): (a) they must be 
numerically stable; (b) the wave motions must not decay in time when 
the dissipation terms are deleted from the governing equations; and 
(c) quadratic properties like total kinetic energy must be 
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conserved. Although the alternating-direction implicit (ADI) method 
used in the 2-dimensional model satisfies these criteria with a 
longer time-step than that allowed by explicit methods, the ADI 
method requires large computer memory and is subject to truncation 
errors when large pressure or momentum changes occur between nearby 
points (Cheng et al, 1976; Leendertse et al, 1973). 


Consequently, an explicit method of integration is used, in which a 
backward, central or forward difference approximation is substituted 
for the time derivative. Equations (8) - (12) are thus integrated 
by a three-time level method in which the forward time-level 
solution depends upon the solution at the two previous time levels. 
This method is known as the central-diff erence, mid-point rule, or 
leap-frog scheme (Cheng et al, 1976). In the first half of the 
computation, the finite-difference forms of equations (8), (9), (11) 
and (12) are used to calculate u, v, § and T at time level (n + 1) 
from the corresponding variables at level (n - 1). The dispersive 
terms [E] and [F] in the momentum equations and [C] in the 
mass-balance equations are evaluated from level (n - 1) while the 
other terms are evaluated from level (n) in order to retain 
computational stability. The vertical velocities, densities, 
pressure gradients and stress terms are next computed from the 
results just obtained for levels (n) and (n+ 1). This cycle is 
repeated as long as required. 


This leap-frog scheme has a rigid time-step limitation, beyond which 
the computation becomes unstable. This is as follows (Grammeltvedt, 
cited by Leendertse et al, 1973): 


Airis Ax 
7 25s, 8 (23) 


where Dass is the maximum water depth. This time step is a 
Shorter by a factor of about 3.5 than the maximum time step 
suggested for the ADI method (Williams and Himwood, 1976; Poulton, 
1977), but the computational advantages of the leap-frog method 
outweigh this disadvantage. 
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For the grid size (360 m) used in most of the Hamilton Harbour 
modelling described herein, equation (23) indicates a maximum time 
step of 20s. Camnputations at this time step, however, were found to 
be unstable, and a time step of 10s was required for stable 
computations. This time step was verified by comparing the results 
of a test computation with those of a similar run using At = 6s; 
the difference between these time step lengths were negligible. 

This method of time step validation is identical to that used 
previously in the 2-dimensional model studies (Poulton, 1977; MOE, 
1979). Time steps shorter than the minimum dictated by equation 
(23) were also used by Leendertse et al (1973, 1975) in their test 
computations. Several spatial smoothing schemes were also advocated 
by Harris and Jelesnianski (1964) in order to reduce the extent of 
numerical instabilities. 


MODEL INPUT TO HAMILTON HARBOUR 


The model as described in the previous section has been applied to 
Hamilton Harbour using a grid of space step length 360 m which is 
1.5 times the 240 m grid used in the previous 2-dimensional model 
studies (MOE, 1980). The axes of the two model grids are parallel 
and the grid points aligned closely so as to allow easy comparisons 
between the two models run with the same input data. This also 
allows comparison with the 3-dimensional model run on the 240 m 
grid; however, at present, the computer costs are too prohibitive to 
permit use of the 240 m grid on the 3-dimensional model. 


The 3-dimensional model grid, together with the real harbour 
shoreline and the locations of source and sink inputs is shown in 
Figure 2. The vertical structure of the model grid consists of 3 
layers of thicknesses 3, 4 and 9 m, respectively. Water depths are 
handled by using one, two or three layers per grid square according 
to the water depth at each point. The number of layers per grid 
point is shown in Figure 3. 


Interactions of the harbour with Lake Ontario (through the 
Burlington Ship Canal) and Cootes Paradise (through Desjardins 
Canal), as well as the mouth of Grindstone Creek, have been treated 
as open boundaries. The model has been programmed to accept current 
meter data as direct input to the model at real-time intervals of 
generally 10 min; the purpose of this input is to describe 
adequately the rapid changes in currents and water quality which 
occur within minutes at such boundaries as discussed in previous 
reports on this subject (MOE, 1974, 1977, 1980; Poulton, 1979; 
Palmer and Poulton 1976). Input from Lake Ontario is entered at the 
upper two model layers, and from the other open boundaries at the 
surface layer only. The real-time data input to the model is 
completed with the specification of hourly wind data as recorded by 
the Hamilton Harbour Police from a site located at the police docks. 


For the calibration attempts described in this section, real-time 
input data has been used fron two periods in which current and water 
quality meters have been operated: (a) June 29-30, 1974. This is 
the same period as used in the 2-dimensional model study (Section 
F); and (b) June 12-14, 1976. For both periods, constant average 
source and sink flow data were used (Table 1). Wind input for the 
1974 period is shown in Figure 4, and for the 1975 period is shown 
in Figure 5. For both periods, current meter data at the Burlington 
Ship Canal were read directly into the model at 10-minute intervals 
for two layers; in addition, for 1976 current meter data from 
Desjardins Canal were also read directly into the model at 10-minute 
intervals. Details of the current meter inputs have been described 
by Kohli (1980) and will not be repeated here. Daily flow records 
from Grindstone Creek and the average input value from Desjardins 
Canal (1974 only) were converted to current speeds by dividing by 
the model grid cross section area (Ax * top layer depth). 


MODEL CALIBRATION 


The model calibration is conceptually similar to that used with the 
2-dimensional model, except that there are more physical parameters 
to adjust with the 3-dimensional model. Modelled velocity 
components are compared at 10-minute intervals to current meter data 
obtained at the internal validation points shown in Figure 2. It 
should be noted that these metered results are used for calibration 
purposes only and are not entered into the model calculations. The 
calculated u and v components were converted into a resultant speed 
and compass direction, and plotted together with the observed 
currents as a function of time for the appropriate current meter. 
Different model grid locations were used in the vicinity of the 
current meter for the test procedure, and the locations selected 
represented the best fit to observed data as obtained after sunming 
the differences between observed and modelled velocities for several 
rows and columns near the validation point. 


The physical parameters and their meaning were described in the 
theoretical section and are summarized here: 


Parameter Meaning 
At | time step length 
AX Space step (model grid size in horizontal 
direction) 
C Chezy's coefficient of bottom roughness 
(eq. 16a and 16b) 
C* coefficient of wind stress (eq. 15a and 15b) 
Ax Ay horizontal momentum dispersion coefficients 


(eq. 8 and 9) 


Vo Vertical momentum transfer coefficient (eq. 
21a and 21b) 


Vy coefficient of interfacial friction 
(Prandtl's mixing length) 
These two coefficients (eq. 19a and 19b) 
determine the effect of density gradients 
(Richardson number, presently set to zero), 
and vertical velocity gradients, on the 
interfacial shear stress terms (7 x2) kh 
and CT Io) Riles 


Attempts at calibration so far, have included only the hydrodynamic 
portion of the model, at constant density. Despite the fact that 
thermal stratification existed during the calibration periods, its 
effect (reduction of vertical turbulence, corresponding to a 
decrease in the value of 44,) has been ignored at present. 


The model time step length was first verified by checking that 
modelled currents were not a function of time step length, as was 
previously done with the 2-dimensional model. The model was run for 
12 h with the 360 m grid. Input parameters included a wind speed of 
6.7 ms-+ from WSW, a Chezy coefficient of 200 cmé 572, pi 

5000 cm? and time step length At = 10 s. The results were 

compared with a similar run using At = 6 s, in order to test the 
effect of time step length. Figure 6 shows a plot of the modelled 
U, V, and w components for the first 10 h from point (17, 6,2), 
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which is the mid-depth of a point near harbour station 

# 270. The time series for the two points are displaced slightly 
for clarity. The results show that 10 s is an adequately short time 
for a model grid of this size; consequently this time step was used 
in further modelling work described here. Three other grid points 
plotted confirmed this conclusion. 


CALIBRATION WITH 1974 INPUT DATA 


Work on the 2-dimensional model (Poulton, 1980) and preliminary work 
on the 3-dimensional model not described here indicated that spin-up 
time is an important factor that must be accounted for in applying 
real-time data to the 3-dimensional model, particularly at low 
initial wind speeds. Under conditions used in the test run for 
verifying At, steady state was approached within the 12 h period, 
with a rough approximation being achieved within the first few 
hours. In order to reduce the computer costs of initial model 
calibration with 1974 real-time input data, it was decided to use 
the final state from this run as initial conditions for model 
calibration runs as the wind conditions used were similar to those 
at the beginning of the 1974 calibration period. 


The values of model physical parameters used in runs involving 1974 
input data are shown in Table 2, and the location of the current 
meter used for comparison purposes is given in Figure 2. The 
initial values of the physical parameters were similar to those used 
by Leendertse et al (1973) in their test basin. During the 
comparisons, two layers of grid points were tested. 


Although the depth (4 m from surface) of the current meter suggests 
the use of the top two model layers, simple observation showed that 
modelled surface currents generally followed the wind direction and 
not the observed current meter results. Layers 2 and 3 of the grid 
were therefore tested. 
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The best agreement found is shown in Figure 7. While the agreement 
in current direction is reasonable after 4 h, the observed speeds of 
5 to 13 cm/s were not simulated, the calculated speeds being fairly 
constant at 2 to 3 cm/s. Also, this comparison was found in model 
layer 3, which was at too great a water depth. 


For a second calibration attempt (run #2), it was decided that 
agreement should occur in model layer 2 if possible. Preliminary 
work had indicated that decreasing the coefficient of interfacial 
fieiction 4 would result in a large difference between layer 1 and 
2 currents, and a smaller difference between layer 2 and 3 
currents. It was thus decided to set 4 = 100 cm*; all other 
parameters were kept unchanged and the model was run for 12 h. The 
best comparison obtained is given in Figure 8. An initial increase 
in current speeds was now found for about 7 h. Speeds, however 
decreased to near zero in the last 5 h. The direction also became 
more variable, exhibiting a reversal in the last hour which did not 
occur in the field data. 


For another attempt at increasing speeds near the bottom, run #3 was 
done with the Chezy coefficient C increased to 700 cm? go 

Figure 9 shows the best agreement obtained. This result shows 
improved agreement compared to runs 1 and 2, and it was felt that a 
closer approach to model calibration was being achieved. 


As the physical parameters used in run #3 were considerably 
different from those in the spin-up run which output was used for 
initial conditions, it was decided to repeat the run from initial 
conditions using the same physical parameters as run #3. The 
real-time data were then applied after a reasonable spin-up period. 


However, the model was found to be unstable when run under these 
conditions for more that 19 h. For improved stability, the value 
of Vv, was increased to 500 cm? s7 and the model run for 24 h 

from rest. This run (run #4) was stable but failed to reach steady 
state at 24 h. As the effect of spin-up time of the order of 12 to 
24 h was found to be significant with the 2-dimensional model in 


Toronto Harbour (Poulton, 1980), it was decided to test this effect 
using two parallel runs (4A and 4B), done with spin-up times of 12 
and 20 h, respectively. Figures 10 and 11 show the best comparison 
with observed data that was obtained (using the same model grid 
point in each case). Run #4B (Figure 11) lost stability at 18.8 h, 
indicating complete stability still did not exist. In addition, the 
direct comparison between the two figures shows that the results are 
significantly different, and in only fair agreement with observed 
data. Clearly the model is not producing meaningful results under 
the conditions used in these runs; the agreement in the earlier run 
(#3; Figure 9) under similar conditions was somewhat fortuitous. 


In an attempt to improve stability, Vi was further increased to 
1000 cmé. This run (run #5) still became unstable; comparisons 
with observed data were similar to runs #4A and 4B and are not shown. 


In these unstable model runs, the modelled currents show large 
spatial variations in direction, frequently between neighboring grid 
points. An example is shown for layer 1, run #4B, in Figure 12. It 
was Suspected that increasingly large differences between adjacent 
velocity components contributed to instability and changes were 
needed to smooth out these variations. For this reason, the 
horizontal momentum dispersion coefficients A. and A. were 

increased to 1 x 10° cm* s~! and the model oe for 4 h with 
real-time input from the run #5 spin-up condition. This run (#5A) 
was stable for 24 h. However, no improvement was obtained in 
comparisions with observed results. At this point it was decided to 
attempt calibrations with a completely different set of input and 
validation data. 


CALIBRATION WITH 1976 INPUT DATA 


An input data set from June 1976 was selected, which was 
characterized by winds from the northeast for at least a portion of 
the record. The physical parameters affecting calibration were set 
to the same values as the last run with 1974 data (which produced 
Stability but not a good comparison with observed results). They 
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were: C = 700 cm? s*, », = 1000 cm, A, = 10 
cm¢s7+, The wind input for 12 hours preceding the model run 

plus 48 hours of modelled time is shown in Figure 5. The location 
of the current meter data used for comparison is shown in Figure 1. 
The model was spun up for 24 h with wind from the northeast at 5.3 
ms~ 2 and then run for 48 h with real-time data as already 


described. 


Comparisons with metered data are shown in Figure 13 for two 
locations in layer 2. This layer was selected for comparison as the 
current meter was located at 6 m depth. At both locations, the 
modelled direction was incorrect for the first 6 h. A fairly close 
correspondence between observed and modelled speeds and directions 
was observed for the remainder of the first 24 h, although the 
variation of model speeds at location (16, 11, 2) (Figure 13a) with 
time was inconsistent with the observed variation. For most of the 
second 24 h, modelled directions were almost opposite those 
observed, and modelled speeds were consistently much too low 
(Figures 13c-d). Only during the last 8 h at location (18, 12, 2) 
and the last 4 h at location (16, 11, 2) did the modelled directions 
approach the observed directions. Other locations within the same 
model layer yielded worse results; as with comparisons involving 
1974 input data, the illustrated locations represented the best 
agreements between observed and modelled data obtained by testing 
all grid points within several grid squares of the meter location 
(Figure 2). Obviously the model in the present condition was not 
able to simulate the observed current speeds and direction changes 
during the last 24 h. Further attempts at model calibration, by 
altering the values of physical parameters, were discontinued at 
this time. 


Despite the lack of agreement between observed and modelled currents 
over much of this time interval, it was felt useful to present a 
series of vector plots of modelled currents to show the predicted 
changes with time as the wind direction (Figure 5) shifts from NE to 
SW. These plots are shown in Figures 14a-e. At zero time (i.e. 
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after 24 h of spin-up with NE winds), surface currents were mainly 
from east to west, showing a deflection to the right of the wind 
‘direction as predicted by the Coriolis effect. Some return flow can 
be found in the deepest part of the west end of the harbour. Mid- 
depth currents in the deepest part of the harbour are deflected 
farther to the right of the wind direction, and form a counter- 
clockwise circulation at the west end of the harbour. As expected, 


bottom currents return against the wind. 


Their speed (see plotting scales on figures) is considerably lower 
than those of the upper two layers; surface currents are the 
fastest. Twelve hours later (Figure 14b), very little change in 
modelled currents has occurred, aS expected due to the constant wind 


direction (Figure 5). 


After 24 hours, the wind direction had begun to shift. Currents 
(Figure 14c) were much slower in the top two layers, although 
directions had not yet changed much. Bottom layer speeds and 
directions had not begun to change appreciably. 


After 36 hours, winds were from the SW (Figure 5) and surface 
currents (Figure 14d) had shifted around to follow the winds. 
Mid-layer currents in the deepest zone were again shifted to the 
right of surface currents, and bottom layer currents were in the 
opposite direction. These layers again both exhibited lower speeds 
than surface currents, especially in the middle layer at the west 
end and southeast corner of the harbour. Six hours later (Figure 
14e), under SW to W winds, the currents in all layers had shifted 
farther to the right; mid-depth and bottom currents were faster 
while surface currents were approximately of the same magnitude. 
Although not quantitatively verified, these results do indicate in a 
qualititative manner the response of the harbour currents to shifts 
in wind speed and direction. Similar observed variations of current 
Speed and direction between depths representative of the upper two 
model layers have been observed during drogue trackings carried out 
during the summer of 1979 in the east end of the harbour (Figure 15). 
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The fact that modelled mid-layer currents are frequently opposite in 
direction to observed currents suggests that bottom layer current 
directions may agree with the observed direction for portions of the 
time. Figure 16 shows a test of this hypothesis. It is seen that 
during the first 6 hours and from 32 to 40 h modelled time the 
current directions are indeed reasonably close to the observed 
direction, although modelled speeds are considerably lower than 
observed speeds. Between 40 and 48 h modelled time, results from 
either the mid-depth or bottom layer approach the observed 
conditions. This result suggests that because of the known harbour 
vertical instability, water flowing past the fixed current meter 
locations may be representative of different model layers at 
different times. Confirmation of this suggestion would require more 
vertical layers in the model and correspondingly increased computer 
costs. 
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DISCUSSION 


The results described in this section illustrate that calibration of 
a 3-dimensional model is a very complicated process. Some of the 
model coefficients were determined from literature surveys or from 
the prior work of Leendertse et al (1973, 1975) with the same 

model. These served to establish a range of values within which a 
realistic approach to calibration might be achieved. Additional 
refinement of the variables is exceedingly difficult and complicated; 
moreover, it requires an adequate field data base including currents 
from several depths in the harbour which were not available in this 
study. Rasmussen and Badr (1979) have applied mathematical 
optimization techniques including least-squares methods and mini-max 
criterion to the 2-dimensional Hamilton Harbour model. Although 
they report reasonable success, they restricted their effort to a 

12 h period in real-time; results could change as the time interval 
is extended. Furthermore, application of their systematic method to 
the 3-dimensional model with its larger number of coefficients, 
becomes a far more cumbersome and time-consuming process. 


The greatest shortcoming in the present work, as illustrated by the 
run with 1976 data, is the limited extent of vertical resolution in 
the model. Results of drogue trackings in the harbour during 1979 
(Figure 15) show that although currents very close to the surface 
are quite close to wind directions, currents within 6 m of the 
surface are considerably slower as well as in different directions. 
Aditional field data of this nature plus additional model layers are 
required to describe the vertical dependence of currents and provide 
a proper value for the coefficient of internal friction Vie This 
parameter is the one wnich has the greatest effect on determining 
the depth dependence of modelled currents and without a realistic 
value of v,, one cannot reasonably adjust coefficients like C and 9. 
It is thus seen that adequate vertical discretization of the model 
plus sufficient validation data as a function of depth are important 
for proper calibration of the model. 
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On the other hand, the method of model-meter comparison as used here 
Should be considered a very severe test of model validity. This is 
true despite the fact that the same method was used with better 
success in the 2-dimensional Toronto Harbour model. Simons (1974) 
used similar plots of observed vs. computed currents for his Lake 
Ontario model, but in addition, used vector plots of 3-day mean 
currents. Lam and Halfon (1978) used monthly mean currents measured 
at several locations in Lake Superior to calibrate their 
3-dimensional hydrodynamic model; results confirmed the general 
circulation pattern but frequently underestimated the speed of 
nearshore currents. They attributed this shortcoming to the lack of 
nearshore grid resolution. In the same paper, they were able to 
apply the hydrodynamic model successfully to a two-compartment 
phosphorus model. In other words, it may not be necessary to obtain 
an accurate prediction of currents at short time scales; estimates 
of the average harbour circulation pattern may be sufficient. 


FUTURE MODELLING WORK 


The purpose of numerical modelling in Ontario is to formulate 
discharge requirements which will ensure that the quality of 
receiving waters meets the criteria established by the Ministry of 
the Environment (Blue Book). Where local violations of the criteria 
occur adjacent to the discharges, models can be used to define the 
limits of limited use zones (mixing zones or zones of non- 
compliance) as defined by the International Joint Commission. Valid 
attempts at environmental assessment require a properly calibrated 
and verified model; moreover, they require an intelligent choice of 
the proper model to fit the problem at hand. Although the 
2-dimensional model of Hamilton Harbour (Section F) has only limited 
usefulness under conditions of thermal stratification, it proved 
considerably more successful when applied to Toronto Harbour 
(Poulton, 1977, 1979, 1980). Either type of model is useful for 
describing only short-term events such as storm runoffs or 
accidental discharges and is of no use in describing longer-term 
(e.g. seasonal) phenomena. 
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The present work has emphasized the need for adequate vertical 
discretization of the model plus sufficient validation data as a 
function of depth. Further model calibration should be done only 
with such data available. The following requirements are suggested: 


1. A minimum of 5 model layers are required, including 3 
layers at the Burlington Ship Canal. 


2. Vertical strings of current meters including temperature 
should be operated in the harbour for validation data, 
including at least three different depths. Open boundary 
data are also required for 3 depths in the Burlington Ship 
Canal. 


3. Wind data are required from a recording anemometer located 
over the water surface and operating at a 10-minute or 
shorter interval. 


4, Air and sediment temperature data are required plus a 
measure of solar radiation as a function of time for a 
realistic temperature calibration. 


In addition to these requirements, extension of the model to include 
turbulent energy computations might be helpful. Leendertse and Liu 
(1977) have produced an extension of the 3-dimensional model 
including sub-grid scale turbulent energy calculations from which 
vertical mass and momentum exchange coefficients are computed. 
Horizontal momentum exchange calculations are also improved by 
computing these as a function of local velocity deformations 
calculated in the finite difference grid. Conputational stability 
was also improved by introducing an implicit method for calculating 
the vertical exchange terms. Applied to Hamilton Harbour, this 
should ease the stability problems observed during conputations with 
runs #4 and #5 as earlier discussed. The computer program for these 
calculations is now available (Leendertse and Nelson, 1977) and need 
only be adopted to the present Hamilton Harbour situation. Before 
these calculations are implemented, however, the value of the extra 
effort which would be required in implementing these calculations 
Should be carefully assessed. 


Any continued modelling program as outlined above naturally 
represents an ambitious and expensive undertaking. Obviously, the 
benefits of such prediction capabilities must be carefully assessed 
before a commitment is undertaken. It would also be hoped that the 
deve loped model could be extended to other areas of concern in 
Ontario such as Toronto Harbour or the Bay of Quinte with minimum 
additional effort once the model has been validated in Hamilton 
Harbour. It is suggested that the model be let out to an external 
contractor such as a university where it could form a profitable 
Ph.D. or post-doctoral program as well as utilization of more 
economical computer facilities. In addition, valuable cooperation 
between the Ministry and the outside agency would be established or 
continued, particularly in the execution of the field program as 
elaborated above. 
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TABLE 1 


SOURCE AND SINK DATA USED FOR NUMERICAL MODEL 
HAMILTON HARBOUR, JUNE-JULY 1974 


MAP I.D. Outfall or Intake Flow 
(m3s-1) 
A Hamilton WPCP 3.68 
B Ottawa St. Slip* 14.06 
CD (Stelco Intake #2) less -12.2 
(Stelco S. Open Cut Sewer) 
a Stelco N. Trunk Sewer 3.04 
2 Stelco Intake #1 -2.40 
G Stelco #3 Open Hearth Sewer 2.40 
H Stelco 148" Plate Mill 0.94 
I Dof asco raw water intake -8.9 
J Dofasco (coke overn by-products), 2.74 
boiler house discharge, silicon steel) 
K Skyway WPCP 0.68 
L Desjardins Canal 0.066** 
M Grindstone Creek 0.64 ** 
Note: 


ex 


Negative flows are intakes (concentrations not required). 


Includes: Stelco (E. side lagoon, HCl regenerator and cold mill, 
hot strip mill, heavy gauge shear line), Dofasco (lagoon, Ottawa 
St. sewer) and CIL (total plant less service water intake). 


These points were treated as open boundaries. Flow data for 

Desjardins Canal (1974 only) and for Grindstone Creek were 

converted to the corresponding current speeds in cm/s by dividing 

by the cross-section area of the one-layer open boundary (3m x 

oS Ws In 1976, current meter data were used for Desjardins 
anal. 
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TABLE 2 
VALUES OF PHYSICAL PARAMETERS USED 


IN REAL TIME 3-d MODEL RUNS USING JUNE 1974 INPUT DATA 


Model Chezy Interfacial Momentum 
Run Spin-Up Run Coefficient Frigtion Dispersion 
# Time, h Time, h C hi Ay, A 
cm s-l cm2 eens 1 
bi i 12 300 5000 500 
2 12 WA 300 100 500 
3 2 12 700 100 500 
4A 12 24 700 500 500 
4B 20 18 .8* 700 500 500 
5 24 20* 700 1000 500 
5A 24 24 700 1000 10000 
Note: All runs used C* (wind stress coeff.) = 0.020 
me (vertical momentum transfer coeff. = 0.0148 cm? sl 


Run blew up at stated time. 


Runs 1, 2 and 3 were done from Run 1 spin-up conditions. 


Runs 4A and B were done from Run 4 spin-up conditions. 
Runs 5 and 5A were done from Run 5 spin-up conditions. 
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FIGURE !| - 3-D MODEL GRID REPRESENTATIONS 
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FIGURE 6a - MODEL U VELOCITIES AT LOCATION (I7, 6, 2) 
COMPARISON OF BESYLIS FOR AT = 6 AND (0 sec. 
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FIGURE 6b - MODEL V VELOCITIES AT LOCATION (I7, 6, 2) 
COMPARISON OF RESULTS FOR AT = 6 AND 10 sec. 
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FIGURE 6c - MODEL W VELOCITIES AT LOCATION (i7, 6, 2) 
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FIGURE 7 - COMPARISON OF OBSERVED AND MODELLED VELOCITIES 
AT LOCATION (II, 19, 3), JUNE 1974 RUN N21 
(C = 300 cm® s~!, »p, = 5000 cm?) 
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FIGURE 8 - COMPARISON OF OBSERVED AND MODELLED VELOCITIES 
AT LOCATION (II, 18, 2), JUNE I974 RUN N° 2 
(Cc = 300 cm® s~', v, = 100 cm®) 
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FIGURE 9 - COMPARISON OF OBSERVED AND MODELLED VELOCITIES 
AT LOCATION U(il, 18, 2), V (IO, 18, 2), JUNE I974 
RUN N° 3, (C = 700 cme s-!, p, = 100 cm®) 
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FIGURE 10 - COMPARISON OF OBSERVED AND MODELLED VELOCITIES 
AT LOCATION (10, 19, 2), JUNE 1974 RUN N°? 4A 
(Cc = 700 cm? s~', py, = 500 cm?) 
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FIGURE 13a - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES AT LOCATION (I6, II, 2), 
JUNE 1976, FIRST 24 hr. 
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FIGURE !3b - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES AT LOCATION u (18, Il, 2) 
v (18, 12, 2), JUNE 1976, FIRST 24 hr. 
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FIGURE I3c - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES AT LOCATION (I6, Ii, 2), 
JUNE 1976, SECOND 24 hr. 
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FIGURE 13d - COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES AT LOCATION U (18, II, 2) 
V (18, 12, 2), JUNE 1976, SECOND 24 hr. 
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FIGURE 14a - VECTOR PLOTS OF MODELLED 
CURRENTS AFTER 24 HOURS OF 
SPIN-UP TIME (AT THE BEGINNING 
OF THE APPLICATION OF REAL-TIME 
DATA, 1200 hr. ON JUNE 12, 1976) 
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FIGURE I4b - VECTOR PLOTS OF MODELLED 
CURRENTS IN THE THREE LAYERS 
AFTER J2 hr. MODELLED TIME 
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FIGURE 14c - VECTOR PLOTS OF MODELLED 
CURRENTS IN THE THREE LAYERS 
AFTER 24hr. MODELLED TIME 
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HAMILTON HARBOUR 3-D SIMULATION LAYER 2. CURRENT VECTOR PLOT 
PLOTTING SCALES: 
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FIGURE I4¢d - VECTOR PLOTS OF MODELLED 
CURRENTS IN THE THREE LAYERS 
AFTER 36 hr. MODELLED TIME 
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HAMILTON HARBOUR 3-0 SIMULATION LAYER 1. CURRENT VECTOR PLOT 
PLOTTING SCALES: Bay 


- a fh O, ~ 
ee Oe eg a WN aly of nga ne 
Bee sae RS RY GO a a MAE eee Ba ¢ 
Eb dele er 7 
a we we 


HAMILTON HARBGUR 3-0 SIMULATION LAYER 2. CURRENT VECTOR PLOT 
PLOTTING SCALES Basal 


HAMILTON HARBOUR 3-0 SIMULATION LAYER 3. CURRENT VECTOR PLOT 
PLOTTING SCALES t “Bass 


FIGURE I4e - VECTOR PLOTS OF MODELLED 
CURRENTS IN THE THREE LAYERS 
AFTER 42 hr. MODELLED TIME 
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FIGURE 15a - RESULTS OF DROGUE TRACKINGS IN HAMILTON HARBOUR ON 
SEPTEMBER 4, 1979, EAST WIND CONDITIONS 
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FIGURE I5b - RESULTS OF DROGUE TRACKINGS IN HAMILTON HARBOUR ON 
SEPTEMBER 27, 1979, SOUTHWEST WIND CONDITIONS 
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- COMPARISON OF OBSERVED AND MODELLED 
VELOCITIES AT LOCATION U (18, !2, 3) 
V (17, 12, 3), JUNE 1976, FIRST 24 hr. 
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FIGURE 


I6b - COMPARISON OF OSSERVED AND MODELLED 
VELOCITIES AT LOCATION (17, 12, 2), 
JUNE 1976, SECOND 24 hr. 


G-6] 


\ ~ 
. oo 


en 


A 


7 


| 


Ot SRAuUoIs & 


ia SMQo 


Th, SOIFRSO Wy TR. S330 59's 


wil 


qj 

= 
o 
ae 
» 
my 
Se 
> 

Cad 
ve 
< 

= 
7 

z 
7 
* 
ad 
° 


yale, 


eS GKOGIh .eFe! 


a , at 


HAMILTON HARBOUR STUDY 


1977 


VOLUME I 


SECTION H 


WATER MOVEMENTS 


SUMMAR Y 


Physical processes of Hamilton Harbour were studied during 1976 and 
1977 by operating self-recording current instruments at a mid- 
harbour location and in the Burlington and Desjardins Canals. In 
the harbour, the resultant currents varied from 0.7 to 1.8 cm.s7 1; 
in Burlington Canal, they ranged between 0.6 and 13.4 cm.s7 and 
in Desjardins Canal, they were from 3.3 to 8.1 cm.s7e, Both the 
connecting canals exhibited faster currents compared to the mid- 
harbour and lake coastal regions, probably due to the constriction 


of the canals. 


Flows through the two canals were estimated in both directions from 
measured currents. The exchange through Burlington Canal was an 
order of magnitude larger than through Desjardins Canal. Net flow 
through Burlington Canal was 0.2 to 1.4% of the harbour volume per — 
day, towards Lake Ontario. Net flow through Desjardins Canal varied 
from 0.07 to 0.15% of the harbour volume per day, towards the 
harbour. 


The spectral analysis of the currents at the mid-harbour showed 
diurnal and semi-diurnal periods and the effects of Lake Ontario. 
The analysis for the two canal locations indicated predominant tidal 
effects as well as the first two oscillation modes of Lake Ontario. 
At Burlington Canal, significant coherences were observed on 
occasions for periods from 2 to 30 h between different depths from 
June to August 1976, including vertical homogeneity of currents and 
water temperature inspite of the weak thermal stratification in the 
canal. Such a vertical homogeneity in other years may or may not 


exist. 


The mean water temperatures in the harbour ranged from 17.6 to 
19.2°C during summer 1976. Burlington Canal mean water 
temperatures ranged from 6.6 to 20.9°C, June to November 1976-77. 

A weak stratification of the Canal was observed during summer 1976, 
followed by no stratification from October to November 1976. The 
mean water temperatures in Desjardin Canal varied from 4.5 to 
20.8°C from May to November 1976. 
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INTRODUCTION 


Hamilton Harbour is a natural harbour located at the western corner 
of Lake Ontario. It contains 280 x 10° m° of fresh water and 

has a mean depth of 13 m. The harbour is connected to Lake Ontario 
by the Burlington Canal and to Cootes Paradise by the Desjardins 
Canal (see Figure 1). The port is very important to heavy industry 
on the south shore. The industry uses and recirculates 27 m>, 
sol (0.8% of harbour volume per day) of harbour water. Several 
creeks discharge 4.1 m°?, s7+ of water to the harbour while the 


municipal sewage outfalls add 3.2 m>. st of treated effluents. 


The Ontario Ministry of the Environment (MOE) has been studying the 
biological, chemical and physical processes of the harbour and how 
their complex interactions affect water quality (MOE, 1974, 1975 and 
1977). During summer of each year, hypolimnion waters experience 
oxygen depletion. Polak and Haffner (1978) concluded that the 
exchange of water between the harbour and the lake provides the much 
needed dissolved oxygen to the harbour. Lake-harbour exchange was 
estimated to be 1% of the harbour volume per day while the net flow 
towards the lake was 0.5% of the harbour volume per day during 
September 1975 (Kohli, 1979). 


This report discusses the results of the 1976-77 current meter 
operations in the harbour, Burlington and Desjardins Canals, 
including lake-harbour and harbour-Cootes Paradise exchange and the 
physical processes responsible for the currents and is a 
continuation of the series of reports on Hamilton Harbour physical 
data (MOE, 1974, 1975 and 1977; Kohli, 1978 and 1979). 


Table 1 presents the current meter operations in Hamilton Harbour, 
Burlington and Desjardins Canals during 1976-77. All current meter 
sites are shown in Figures 1 to 4 and 7. Locations 1117, 1116, 1101 
and 1102 were on the same tower in Burlington Canal at different 
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depths (see Table 1). All instruments were operated in the 
epilimnion, except 1102 which was in the hypolimnion. All data 
collected were numerically smoothed (Kohli, 1978) and then 
partitioned into monthly data sets (records). Statistical analyses 
of water movements and temperatures were performed to determine the 
spatial and temporal characteristics. Currents at the Bur lington 
Canal locations were examined to compute the lake harbour exchange 
using the excursion episode method (Kohli, 1979). Estimates of the 
flow through Desjardins Canal were made using average currents and 
the approximate cross-sectional area of the canal. 


RESULTS# 
CURRENTS 
HAMILTON HARBOUR 


Figure 1 shows the resultant southerly direction of currents during 
the stratified summer months of June to August 1976-77. Similar 
observations were made during August and September 1975 (Kohli, 
1978). The resultant currents at the mid-harbour location (1104) 
varied from 0.7 to 1.8 cm.s! while the arithmetic average speed 
ranged from 1.7 to 4.8 cm.s72 (see Table 2). The maximum speed 
recorded at this location during 1976-77 summers was 28 cms”, 

The persistence factor varied from 0.11 to 0.40 during this study. 
The results of the present investigation compared favourably with 
the previous study (Kohli, 1978). 


BURLINGTON CANAL 


Table 3 and Figures 2-4 present the sunmary of currents in 
Burlington Canal during 1976 and 1977. Figures 2 and 3 indicate 
that the resultant currents at the top two locations, 1117 and 1116 
(8.5 m and 7.5 m from bottom) were towards the lake for the entire 
period of study. However at 1101, 6.1 m from the bottom, the 
resultant currents were toward the lake from June to November 1976, 


Detailed current frequency tables for the various current meter 
locations have been published in Kohli, 1980 
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but they were toward the harbour from May to July 1977 (see Figure 
4). The direction of net currents at 2.1 m from the bottom (the 
lowest location 1102), was towards the canal wall during June and 
July 1976. At this site (1102), the recorded currents were 
integrated over 10 minutes and were less than 2 cm.s72 for 
approximately 96% of the time during June and July 1976. Therefore, 
most of the recorded speeds were actually less than the threshold 
speed of the Plessey Current Meters (3 cm.s7!), When the current 
speed is less than the instrument threshold speed, the direction 
vane can assume any random position and the direction recording 
becomes meaningless. Therefore, the computed resultant current 
direction towards the canal wall, in this case, was meaningless. 
Thus it may be concluded that the net currents in the canal are 


predominantly towards the lake. 


At location 1117 (8.5 m from bottom) the resultant currents ranged 
from 0.9 to 13.4 cm.s7+ while the arithmetic average speed varied 
Racoma7e. GAO 71627 cm.s71 (see Table 3). The maximum speed recorded 
during the study ranged from 25 to 83 cm.s7t. At location 1116 

(7.5 m from bottom), the resultant currents varied from 2.3 to 6.8 
cm.s7! while the arithmetic average speed ranged from 7.8 to 14.7 
cm.s7/ (see Table 3). The maximum speed recorded during the study 
varied from 48 to 129 cm.s7t, The resultant current at 6.1 m from 
the bottom (location 1101) varied from 0.6 to 7.6 cm.s~+, while 

the arithmetic average speed ranged from 6.9 to 14.1 cm.s72 (see 
Table 3). The maximum speed recorded at this level was between 44 
and 126 cm.s7+. The resultant speed at 2.1 m from the bottom 
(location 1102) was 0.5 cm.s7+ during June and July 1976. The 
arithmetic average speed during June 1976 was 0.6 cm.s72 and 

during July 1976, it was 0.7 cm.s7/, The level being closer to 

the botton registered the slowest currents in the canal. More 
specifically, the currents here were negligible for 49% of the time 
and less than 2 cm.s7! for 95% of the time during June 1976. In 
July 1976, the currents were negligible for 30% of the time and less 
than 1 cm.s7? for 70% of the time. The maximum speed recorded was 
25 cm.s~! during June 1976 and 3 cm.s + during July 1976. In 
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view of the very small currents (less than the instrument threshold 
speed of 3.0 cm. s!) persisting over longer periods of time, the 
results of the lowest location (1102) may be used but with due 


caution. 


The results of the top three current meters (1117, 1116 and 1101) 
are comparable to the earlier study (Kohli, 1978). The resultant 
currents of 0.6 to 13.4 cm.s ~ (present study) compare well with 
Aid toolse2témee™e during 1972 to 1975 (Kohli, 1978). The maximum 
speed recorded during the present study was 129 cm.s ~~ which 
compares wel] with the 144 cm.s ~ of the previous study. These 
faster currents in the canal may be attributed to the canal 
constriction (van de Kreeke, 1976) and the Helmholtz mode of the 
harbour (Freeman, Hamblin and Murthy, 1974). 


Careful examination of the current statistics of Burlington Canal 
(see Table 3) shows that the magnitude of the resultant current, as 
well as the arithmetic average speed, decreases with depth. Figures 
2, 3 and 4 show that the resultant currents in Burlington Canal were 
generally towards the lake at all three levels, except during May to 
July 1977 at location 1101 (see Figure 4), when the resultant 
currents were going towards the harbour. This shift in direction of 
the resultant currents at 1101 and the presence of a layered flow 
system in Burlington Canal during summer 1977, in contrast to the 
unidirectional flow during the summer 1976, was confirmed by the 
mean water temperatures (see Figure 5). A sharp temperature drop of 
8-9°C was observed between locations 1116 and 1101 during summer 
1977 while the temperature drop between the same locations and 
during the same period of 1976 was 19°C 


Flows through Burlington Canal were computed by the excursion- 
episode method (Kohli, 1979) and the results for the period June to 
November 1976, are presented in Table 4. As expected, the net 
exchange was towards the lake. The total exchange varied from 1.1 
to 1.6% of the harbour volume per day, while the net exchange 
towards the lake ranged from 0.2 to 1.4% of the harbour volume per 
day (see Figure 6). These results compare well with the previous 
estimates when the total exchange was 1.0% of the harbour volume per 
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day, while the net exchange towards the lake was 0.5% of the harbour 
volume per day (Kohli, 1979) during September 1975. Harris et al 
(1979) computed the average mass exchange through Burlington Canal, 
based on the mass balance of total dissolved solids in the harbour 
over a year. Their results are compared with the present study and 
Kohli (1979) in Table 5. Harris' calculations assume the harbour to 
be in steady state over a year - an unreasonable assumption in an 
otherwise dynamic harbour. Table 5 shows that Harris' flow 
estimates towards the lake and the harbour are 2 and 7 times larger 
than the excursion-episode estimates, but their net flow estimates 
toward the lake are 40% of those by the excursion-episode method. 

As the excursion-episode method computes the flows from direct 
measurements in the canal, taking account of the periodic transport, 
it is considered a better method. 


DESJARDINS CANAL 


Desjardins Canal is located at the western end of Hamilton Harbour 
(see Figures 1 and 7), connecting the harbour to Cootes Paradise. A 
Single current meter at mid-depth (Tocation 1109) operated under the 
railway bridge, from 27 May to 24 November 1979. Only one 
instrument could be installed due to the shallow depth of 3.7 m. A 
summary of current statistics is presented in Table 6. The 
resultant currents varied from 3.3. to 8.1 cm.s7 2, while the 
arithmetic average speed ranged from 9.8 to 24.2 cm.s72. The 
maximum speed recorded was 88 cm.s7 2. The faster currents in the 
canal may be attributed to the constriction (van de Kreeke, 1976). 


Based on the average current speed towards the Harbour and Cootes 
Paradise, with the corresponding percentages of occurence, the flow 
in both directions was computed for the period May to November, 1976. 
As such estimates in Burlington Canal were 2 to 3 times larger than 
the excursion-episode method (Kohli, 1978), a correction factor of 
2.5 has been applied to obtain the more realistic estimates of the 
flow presented in Table 7. Although the correction factor was 
transposed from Burlington Canal to Desjardins Canal, it produced 
results that were compatible with the water balance of Cootes 


Paradise (see discussions below). The total exchange through the 
canal varied from 0.07 to 0.15% of the harbour volume per day (see 
Figure 8), while the net exchange was 0.02 to 0.05% of the harbour 
volume per day (average net flow of 1.11 m says towards the 

harbour. The exchange through the Burlington Canal is an order o 


magnitude larger than through the Desjardins Canal. 


Spencer Creek discharges an average of 0.81 ms? from June to 
November (1961-75), while the Dundas STP discharges an average of . 
Oxte m>.s7! into Cootes Paradise (Semkin et al, 1976). Thus, 
Spencer Creek and Dundas STP, together add 0.93 m*.s + to Cootes 
Paradise while the net average outflow from Cootes Paradise to 
Hamilton Harbour was computed as 1.11 mesa The small 

difference of 0.18 m’s,~2 may be attributed to other smaller 

creeks flowing into Cootes Paradise, ground water, precipitation or 
to the choice of the specific correction factor used for computing 


flows in Desjardins Canal. 
AUTOSPECTRA 
HAMILTON HARBOUR 


Table 8 presents a sunmary of the major spectral periods observed at 
the mid-harbour location (1104). Diurnal (24.0 h) and semi-diurnal 
(12.0 h) periodic motions were observed in the east-west directions 
during July 1976. Semi-diurnal oscillations were also observed 
along the north-south direction during June 1976 and July 1977. 
These are associated with the lunar tidal motions. The 17.1 h 
(August 1976) and 20.0 h (June 1977) periods along the north-south 
direction were due to the lake-wide seiche. In June 1977, along the 
east-west direction, 3.6 and 4.6 h periods were significant and 
these may be due to the lower modes of free oscillations of Lake 
Ontario (Rockwell, 1966). The effects of Lake Ontario are therefore 
observed in the harbour. 


BURLINGTON CANAL 
A summary of significant spectral periods in Burlington Canal is 


presented in Table 9. The 12.0, 5.2 and 3.2 h were the most 
predominant periodicities observed at locations 1117, 1116 and 1101 
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for all monthly data sets and in both major directions. The semi- 
diurnal periodicities may be due to the tidal motions. The 5.2 h 
and 3.2 h are the first and second modes of free oscillations of 
Lake Ontario (Rao and Schwab, 1974 and Palmer and Poulton, 1976). 
At the bottom location 1102 in Burlington Canal, longer periods of 
24.0, 20.0 and 10.9 h were observed. The periodicities of the 
Burlington Canal indicate that the currents are largely affected by 
tidal motions and first and second modes of free oscillations of 
Lake Ontario. 


Summary of 95% significant coherences is presented in Table 10 for 
Burlington Canal water movements and temperatures at 3 levels 
(locations 1101, 1116 and 1117) during June to August 1976. 
Significant coherences between current speed and water temperatures 
at the same location were generally present between 2.0 and 17.1 h. 
Current speed between two locations exhibited coherences from 12.0 
to 30.0 h. Coherence for water temperature between two levels 
ranged from 2.3 to 24.0 h. These significant coherences between 
water currents and temperatures at any two levels indicate the 
vertical cross-correlations of currents and temperature in 
Burlington Canal during June to August 1976. It is not known if 
such vertical homogeneity of currents and water temperature existed 
during the same period of other years. 


DESJARDINS CANAL 


Table 11 summarizes the major spectral periods in Desjardins Canal 
from June to November 1976. The most commonly occurring periodi- 
cities in this canal are 12.0, 5.2 and 3.2 h. The 12.0 h period is 
due to the tidal motions. The other two periods are the first and 
second modes of free oscillations of Lake Ontario. In November 
1976, a 13.3 h period was observed along the channel direction; this 
may be due to the lake-wide seiches. 


It is interesting to note that the Burlington and Desjardins Canals 
have similar periodic motions, namely the tidal motions and the 
first two modes of free oscillations of Lake Ontario. Therefore, 
the water movements in the two canals were caused by the same 


phenomena. 
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WATER TEMPERATURE 


HAMILTON HARBOUR 


At the mid-harbour location (1104), water temperatures were measured 
from June to August 1976. The mean water temperatures varied from 
1/265 0 1 vse progressively increasing from June to August (see 
Table 12). The standard deviation of the temperature decreased from 
1.1°C in June to 0.6°C in August 1976. The minimum recorded 
temperature during a month increased from 14.3 to 16.7°C, while 

the maximum recorded temperature during a month remained almost 
constant at 20%. 


BURLINGTON CANAL 


Table 13 presents the frequency of occurrence of water temperatures 
in Burlington Canal during 1976-77. At the top location 1117, mean 
temperatures varied from 6.6 to 20.9% with standard deviations 
ranging between 1.4 and 2.4°C from June to November 1976. At the 
next location down (1116), from June to August 1976 and May to July, 
1977, the mean temperature ranged between 15.2 and es with 

the standard deviations from 1.4 to 2.6°C. The mean temperature 
varied from 6.6 to 19.3%, with the standard deviation ranging 

from 1.2 to 2.7°C at the mooring 1101, from June to November 1976 
and May to July, 1977. During June and July 1976, at the lowest 
site 1101, the mean water temperatures were 11.0 and 927°: with 
Standard deviations of 2.4 and 2.9% respectively. Figure 5 
presents the mean temperature graphs for June and July of 1976 and 
1977, at 3 levels. The figure shows a very sharp decline in mean 
temperature with depth during 1977, compared to the 1976 curve. 


DESJARDINS CANAL 


Mean water temperature in Desjardins Canal from May to November 76 
varied from 4.5 to 20.8°C (see Table 14), while the standard 
deviations ranged from 1.0 to 3.5%. The maximum temperature 
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recorded in the canal was 26.9°C during the study period. As the 
canal is generally shallow, with a maximum depth of 3.7 m, 
isothermal conditions are expected to prevail. No bathythermographs 
were taken to support the assumption of isothermal regime. 


CONCLUSIONS 


The currents at the mid-harbour location during 1976-77 were similar 
to those observed during 1975. The resultant currents varied from 
O¥7*+to-I's cna while the arithmetic average speed ranged from 
PP? 201 4t8 anes oy. In ge Canal, the resultant currents 
varied from 0.6 to 13.4 cm.s -, while the arithmetic average speed 
ranged between 6.9 and 16.7 cmas7*. The maximum speed recorded in 
Burlington Canal during 1976-77 was 129 émgst ty During 1976 the 
resultant currents in Desjardins Canal varied from 3.3 to 8.1 
cm.s7t, while the arithmetic average speed ranged between 9.8 and 
24.12 cm.s7?, Both canals exhibited relatively faster currents 
compared to mid-harbour and lake coastal currents, probably due to 
the constrictions. 


Net flow through Burlington Canal was estimated as 0.2 to 1.4% of 
the harbour volume per day, while the total exchange through the 
canal ranged from 1.2 to 1.6% of the harbour volume per day. Net 
flow through Desjardins Canal to Hamilton Harbour was estimated as 
0.02 to 0.05% of the harbour volume per day, while the total 
exchange through the canal varied from 0.07 to 0.15% of the harbour 
volume per day. The exchange through the Burlington Canal was an 
order of magnitude larger than through Desjardins Canal. 


Diurnal and semi-diurnal motions were observed at the mid-harbour 
location. Lake-wide seiches and modes of free oscillations of Lake 
Ontario were also present, indicating the influence of Lake Ontario 
at the harbour location. Both Burlington and Desjardins Canals 
exhibited the semi-diurnal periodicities, as well as the first two 
modes of free oscillations of Lake Ontario. Thus, the tidal motions 
and lake effects appear to influence the water movements in the two 


canals. Significant coherences were observed from 2 to 30 h between 
different levels in Burlington Canal during June to August 1976, 
indicating vertical cross-correlations of currents and water 


temperatures in the canal. 


Mean water temperatures in Hamilton Harbour varied from 17.6 to 
Peat during June to August 1976. The maximum temperature 
recorded at this location was 20°C. In Burlington Canal, the mean 
water temperatures varied from 6.6 to 20.9°C from June to November 
1976, and May to July 1977, with the standard deviations of 1.2 to 
2.7°C. The Canal had a weak stratification during summer 1976 and 
a strong stratification during summer 1977. No stratification was 
observed in the canal from October to November 1976. In Desjardins 
Canal, from May to November 1976, the mean temperature varied from 
1560 20.8°C, while the standard deviations ranged from 1.0 to 
3.5°C. Isothermal regime was assumed in the canal. 
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TABLE 8: Summary of Major Spectral Periods (Hours) 
Hamilton Harbour, Lake Ontario, 1976-77 


(95 Percent Confidence Level) 


LOCATION | PERIOD NORTH-SOUTH EAST-WEST 


12.0* None 


None 76.0.242.0, 6.0% 
17.1%, 6.0*, 4.8*, 4.3*| None 

200" 6 dU. Ues need 4.6*, 3.0% 

12.0 None 

None Me gs 


* 80 Percent confidence level 


TABLE 9: Summary of Major Spectral Periods (Hours) 
Burlington Canal, Lake Ontario 1976-77 
(95% Confidence Level) 

LOCATION PERIOD ALONG CHANNEL ACROSS CHANNEL 
Jun 76 12 005/28 oe2 | 1220, bee eaee 
Jul 76 12708 5 oe ance fee eee Parke sre 
Aug 76 | 5.2, 3.2* 12.0, 5.24 3.2 
Oct’ 76.-| “12.0%, 5.2, 1.2 1250), 5ezeeoee 
Nov 76 | 12.0, 5.2, 3.4 | | 1220, 5.2) 3.5% 
Uunip7 ee lee. Onnbe cemeeee 12505 5320e ec 
JWT. Jale le nkes 0 eibeoearced 4) 12,0, Seeeeiee 
dun 76" ||" 1200%). 5.0n13t2e 12,0*S 5. 2.0c0es 
dutevome | le Ome sac 12705-5722 3.2% 
Augu76 © | 412.0s,5 62.8 | 420, 5.24 352% 
Oct 76 | 12.0, 5.2, 3.2 | 12.0, Secenegee 
Nov 76 TORO. 502, 302 1250576. 7s cemeeas 
Jun 77 20.0*, 12.0, 8.0*, |  12.0*, 8.0%, 5.5%, 3.2% 

a 52k ee : 

Jul 77 WS BORD |  20.0*, 12.0*, 5.5% 
un 7b) tis oNOee 15 42g 28 12.0*, 12.0*, 5.5% 
JuMM760y | 2c 0s eee sn 2 i) 2.0 5e2 
Aug 76 12.05.5,25 302% 1250952 3c 
Oct 76 1250 5 ocaee* 1280 5 2s aes 
Nov 76 12200525 go texk 1000, 5225. 352% 
Jun 77 12.0*, 6.0%, 5.2 1. Oe5e2 
Jul 77 12%0% 12.0*, 5.2% 
Jun 76 5.2*, 4.8 5.2% 
Jul 76 20.0% 24.0*, 10.9% 


* 80 Percent Confidence Level 
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TABLES eUL: Summary of Major Spectral Periods (Hours) 
Desjardin Canal, 
Hamilton Harbour, Lake Ontario 1976 


(95% Confidence Level) 


LOCATION PERIOD ALONG CHANNEL ACROSS CHANNEL 


1109 . -UUna/Oa "1S 12205 2620 | 270; *540 
| Jul 76 12806502 | 12.0, 6.0*, 5.2, 3.2* | 
Aug 76 1260805 23.2% 12.08 5.2... 3s5.enoee 
Sep 76 12.0*, 6.3*, 4.8 1200", 25142 
Octa/6. |p 1280s Sev sos L270, 522, 3Q* 
Nov 76 133%, 572 12-0, “bee 


* 80 Percent Confidence Level 


TABLE 12: Temperature Frequency, 
Hamilton Harbour, Lake Ontario 1976 
Percentage of Occurrence 
Location Code 1104 


Temper ature 
Range °C 


Mean OC 

| Std. Dev.°C 

Minimum OC 

Maximum OC 

Series Length 
(h) 
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